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THE  GCARS  SYSTEM 
FORTRAN  IV  PROGRAMMERS  MANUAL 
PART  A  -  PROGRAMS   FOR  DATA  COLLECTION  AND 
EVALUATION 

by 

A.  Keith  Turner 

INTRODUCTION 
The  Generalized  Computer-Aided  Route  Selection  (GCARS )  System  has  been 

developed  to  assist  engineers  in  defining  optimal  highway  alternatives  in 

i  \* 

terms  of  many  conflicting  highway  location  factors  (5,  10)  .  The  use  of 

the  computer  dictates  that  information,  or  "measures",  of  the  various  factors 

which  are  considered  pertinent  to  the  particular  analysis  be  reduced  to 

numerical  models.  This  generally  involves  sampling  of  areal  data.  The 

entire  subsequent  analysis  depends  on  the  quality  of  the  information 

collected;  thus  the  development  of  simple,  logical,  and  reproducible  sampling 

procedures  is  an  important  prerequisite  to  the  successful  application  of 

the  GCARS  System.   In  addition,  procedures  for  testing  the  sampled  data  for 

accuracy  and  for  lack  of  bias  in  depicting  the  actual  conditions  form  a 

necessary  part  of  the  GCARS  System. 

This  part  of  the  Programmers  Manual  describes  the  computer  routines 

developed  to  (l)  assist  in  the  design  of  sampling  procedures  for  topographic 

data  (Program  VECTOR),  (2)  convert  card  data  to  standard  GCARS  System  tapes 

(Programs  LOADS  and  LOADG),  and  (3)  test  sampled  data  for  accuracy  and  lack 

of  bias  (Programs  TREND  and  RESMAP).  Parts  B  and  C  of  this  manual  document 

the  programs  for  other  phases  in  the  GCARS  System  (7,  8).  Overall  System 

design  concepts  and  standards  are  summarized  in  the  GCARS  Users  Manual  (6). 


Numbers  in  parenthesis  refer  to  references. 
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METHOD  OF  OPERATION 

During  development  and  testing  of  the  prototype  GCARS  System,  large 
volumes  of  data  describing  the  various  highway  location  factors  had  to  be 
reduced  to  numerical  form.  Many  factors  were  shown  in  graphical  form  on 
various  types  of  maps.  Figure  1  shows  the  IAKR-V  digital  coordinatograph 
used  to  "digitize"  the  map  data.  This  machine  punches  the  data  onto  cards 
through  the  attached  card  punch. 

Figure  2  shows  the  interconnections  and  flow  of  operations  through  the 
various  programs  and  subroutines  described  in  this  manual.   Program  LOADS 
was  developed  to  allow  the  conversion  of  large  unwieldy  card  data  decks 
produced  by  the  LARR-V  to  more  convenient  tapes.   Program  LOADG  performs 
similar  functions  for  girdded  data  arrays  produced  later  in  the  analysis. 

The  remaining  Programs  perform  various  statistical  and  display  functions 
which  are  designed  to  assist  the  engineer  in  evaluating  his  data. 

PRESENT  CAPACITIES,  OPTIONS,  AND  RESTRICTIONS 

a)  Size  of  Matrices 

The  GCARS  System,  as  presently  developed,  has  been  standardized  to 
work  with  50  x  50  matrices.  Neither  dimension  can  exceed  these  limits 
without  correcting  the  DIMETEION  statements  in  the  programs.   Changes  in 
shape,  for  instance  to  100  x  25,  can  be  accomplished  without  increasing  the 
overall  storage  requirements . 

Program  Vector,  involved  with  the  quantification  of  small  land  form 
units,  is  dimensioned  to  accept  only  a  20  x  20  array  of  elevations;  larger 
arrays  were  not  believed  necessary. 

b)  Irregularly -Spaced  Data 

Program  TREND  performs  regression  analysis  on  large  numbers  of  irregularly- 
spaced  data-points.  The  maximum  number  of  such  observation  that  can  be  processed 
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LEGEND 


A      THE   READING    TABLE 

1  X-axis  arm 

2  Y-oxi»  arm 

3  keyboard  arm 

4  keypack 

5  cursor  carrioge 

6  readout  button 

7  X-axis  encoder 

8  Y-oxis  encoder 

9  magnetic  hold-down  strips 

10  roller  for  long  charts 


B      TELECORDEX    UNIT 

11  X-axis  telepack 

12  Y-axis  telepack 

13  event  counter 

14  memory  and  control    unit 

15  power  switches 

16  constant  data  switches 

17  patch  panel 

C      IBM   MODEL    029    CARD  PUNCH 


FIGURE 


THE   LARR-V    DIGITAL    COORDINATOGRAPH 
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at  once  is  very  large;  over  10,000  observations  have  been  successfully 
analyzed  by  this  program.  In  other  programs,  1000  irregularly  spaced 
observations  is  generally  the  maximum  number  that  can  be  handled. 

c)  Input  Options 

Both  Programs  TREND  and  RESMA.P  can  accept  data  on  either  cards  or  tape. 
For  purposes  of  consistency  the  various  types  of  tapes  are  given  different 
but  consistent  names  (and  numbers)  with  all  these  programs.  They  are  as 
follows:       TAPE1  -  any  temporary  "scratch"  tape 

TAPE2  -  a  tape  containing  irregularly  spaced  observations. 
These  tapes  are  blocked  into  1000  observation 
(3000  word)  logical  records. 
TAPE3  -  a  tape  containing  gridded  data  matrices. 
TAPEU  -  a  tape  containing  original  values,  computed  values, 
and  residuals  for  all  observation  points.   It  may 
be  generated  on  request  by  Program  TREND. 
TAPE5  (SYSIN)  -  is  the  computer  system  input  file  from  the 

card  recorder. 
TAPE6  (SYSOT)  -  is  the  computer  system  output  file  to  the  printer. 
PUNCH  -  is  the  output  file  to  the  card  punch. 

d)  Output  Options 

Programs  TREND  and  RESMAP  allow  many  output  options  including  different 
forms  of  printer  display,  punched  cards,  and  tape. 
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METHOD  OF  COMPUTATION 
a)  Program  VECTOR 

The  terrain  through  which  a  highway  may  pass  is  often  an  important 
location  factor.  In  extremely  rugged  areas  it  may  be  the  dominant  factor. 
Topographic  maps  are  one  type  of  graphical  measure  of  the  earth's  surface. 
This  surface  is  not  uniformly  easy  to  sample  however.  Closely  adjacent 
points  can  have  a  great  variety  of  spatial  relationships  which  collectively 
form  the  "topography"  of  an  area.  Topography  is  classifiable  into  any  number 
of  categories  which  are  frequently  given  names  such  as  "smooth",  "rolling", 
"hilly",  "crested",  etc.  Topographic  maps  show  these  differences  by  contour 
patterns  of  different  types.  It  is  possible  to  study  topographic  maps  and 
subdivide  them  into  areas  containing  different  types  of  topography.  Now 
each  subarea  has  an  associated  degree  of  difficulty  of  sampling;  obviously 
the  more  difficult  areas  should  have  a  greater  sampling  density  than  the 
easier  areas  so  that  all  portions  of  the  map  are  sampled  with  a  similar 
reliability.  Some  quantitative  measure  of  this  sampling  difficulty  is 
required. 

The  development  of  a  suitable  quantitative  measure  of  terrain  roughness 
utilized  vector  analysis  procedures.   Program  VECTOR,  a  modified  version  of 
a  computer  program  of  the  same  name  developed  by  Hobs on,  performs  the 
calculations  (3).  The  computations  are  quite  simple.  Within  a  rectangular 
array  of  elevation  readings,  a  set  of  intersecting  triangular  planar  surfaces 
is  defined  by  groups  of  three  adjacent  values,  as  shown  in  Figure  3-A  and  3-B. 
Normals  to  these  planar  elements  are  represented  by  unit  vectors.  A  statistic, 
known  as  Fisher's  dispersion  factor  (K)  indicates  the  spread  of  these  unit 
vectors  (2,  3,  9).  It  takes  on  small  values  for  highly  dispersed 
distributions  characteristic  of  "rough"  areas  (Figure  3-D)  and  large  values 
for  low  dispersions  characteristic  of  "smooth"  areas  (Figure  3-C). 
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DEFINITION     OF     TRIANGULAR     PLANES     FROM     EliVATlON     MATRIX 
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TYPE  2    TRIANGLES 


"SMOOTH"    TERRAIN 


LOW    VECTOR    DlSPFRSlON 


ROUGH      TfcRRAIN 


HIGH    vrcTOR    DISPERSION 


(adapted   from   Hobson) 


FIGURE    3.     PROGRAM    VECTOR   ANALYSIS    PROCEDURES 
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The  technique  was  tested  on  25  terrain  types  (9).  As  a  result, 
Fisher's  dispersion  factor  (k)  was  found  to  be  an  excellent  descriptor  of 
terrain  roughness.  The  log  K  transformation  was  found  useful  and  could  be 
correlated  with  other  topographic  measures  such  as  log  (variance  of  dips) 
and  log  (maximum  dip). 

However  the  roughness  does  not  completely  measure  the  sampling  difficulty. 
A  perfectly  smooth,  planar,  horizontal  surface  can  be  described  by  a  single 
elevation  value;  any  planar  inclined  surface  requires  three  elevation  values, 
and  associated  X  and  Y  coordinates,  to  describe  it.  Gently  curved  surfaces 
require  still  more  data  points  to  satisfactorily  describe  their  overall 
orientation.  The  overall  orientation  of  subareas  might  be  described  in  a 
number  of  ways.  The  range  of  the  elevations  is  the  simplest  value  to  obtain, 
requiring  no  prior  sampling  procedure.  Accordingly  a  terrain  variability 
factor,  which  includes  both  roughness  and  range  of  elevations,  is  defined 
as  follows : 

i   -U--1-J.  j>     u.  Range  in  elevation 

terrain  variability  factor  =  ~= — ^ r— \ 

Log  (roughness  ; 

Values  of  this  factor  are  large  in  "rough"  areas  and  small  in  "smooth"  areas. 
The  terrain  variability  factor  is  used  in  designing  a  sampling  procedure  for 
a  topographic  map.  The  technique  is  described  fully  by  Turner  (5). 

Program  VECTOR  computes  the  terrain  roughness  and  the  terrain  variability 
factor  for  any  elevation  matrix  submitted  to  it.   Its  associated  subroutines 
TILTS,  WULFF  and  HISTO  are  designed  to  give  additional  statistical  descriptors 
of  the  terrain  which  may  be  used  in  more  carefully  analyzing  terrain  (9). 
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b)  Program  LOADS 

Program  LOADS  reads  X,  Y,  and  Z-coordinate  data  from  punched  cards 
and  writes  the  data  onto  tape.  The  tape  is  developed  according  to  GCARS 
standards  -  one  record  containing  title  information  followed  by  a  series 
of  logical  records  each  containing  1000  sets  of  observation  data.  Since 
each  set  of  data  contains  X,  Y,  and  Z  coordinates,  the  records  are  3000 
words  long.  Since  the  number  of  observations  is  unlikely  to  be  an  exact 
multiple  of  1000,  the  last  record  contains  any  remainder  and  is  therefore 
of  variable  length.  All  data  is  in  relocatable  binary.  An  end-of-file 
mark  is  placed  at  the  end. 

Program  LOADS  can  accept  two  types  of  data,  integer  or  floating  point 
Z-coordinates.   Integer  Z-values  occur  when  data  such  as  soil  conditions 
have  been  sampled  with  an  integer  code  identifying  each  soil  type.  LOADS 
converts  such  integer  codes  to  the  corresponding  floating  point  value  (FAC) 
supplied  as  control  data.  These  may  represent  one  of  several  evaluation 
schemes.   Thus  Program  LOADS  allows  for  the  simultaneous  conversion  of 
measures  to  values  during  the  transfer  of  data  from  cards  to  tape.  When 
floating  point  2-values  are  indicated  no  conversion  is  attempted.   Card 
formats  can  be  varied,  and  more  than  one  set  of  observations  can  be 
included  on  a  card  (the  maximum  is  ten). 

Program  LOADS  can  also  accept  any  number  of  data  subsets,  separated 
from  each  other  by  end-of-record  cards,  and  combine  them  into  one  large 
data  set.  The  identities  of  the  individual  subsets  are  lost  in  this  process, 
A  particular  1000-set  record  may  contain  values  from  several  subsets. 

Data  in  each  subset  may  be  referenced  to  its  own  origin,  with  the 
following  restrictions:       (l)  all  X-coordinates  are  assumed  to  increase 

toward  the  right 
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(2)  all  Y- coordinates  are  assumed  to  increase 
upward . 

(3)  all  coordinate  systems  are  assumed  to  have 
a  common  orientation. 

Each  subset  includes  one  control  card  giving  the  X  and  Y  coordinates 
of  its  origin  in  the  master  coordinate  system  and  a  scale  factor  used  to 
convert  the  subset  coordinates  to  a  consistant  scale  (Thus  different  subsets 
may  have  different  scales).  Figure  k   shows  how  the  coordinate  conversions 
are  specified.   This  option  allows  data  to  be  collected  from  a  number  of 
adjacent  maps  then  combined  to  produce  a  single  data  set  covering  a  larger 
area, 
c)  Program  TREND 

Trend  surface  analysis,  a  variety  of  regression  analysis,  is  used  in 
the  GCARS  System  to  test  value  surfaces  for  adequacy  and  lack  of  bias.  Trend 
surface  analysis  is  widely  used  in  the  earth  sciences  to  separate  regional 
from  local  effects  for  any  type  of  spatial  phenomenon.  The  technique 
involves  the  fitting  of  comparatively  simple  trend  surfaces  to  data,  according 
to  the  least-squares  criterion.  Statistical  measures,  such  as  partial  and 
total  F-tests  (l,  5)  and  the  simple  correlation  coefficient,  are  used  to 
estimate  the  significance  of  the  surfaces.  One  or  more  of  these  surfaces 
are  selected  to  represent  the  regional  effects.   The  residuals  to  the  surfaces, 
representing  the  lack  of  fit  of  the  surface  to  the  data,  ere  plotted  as  contour 
maps  and  are  studied  for  information  concerning  local  anomalies.   This  technique 
is  described  in  detail  by  Turner  (5),  and  in  the  geological  literature. 

Trend  surface  analysis  Is  also  used  to  convert  topographic  data  to  earth- 
work cost  values.  Estimation  of  probable  earthwork  costs  throughout  an  area 
requires : 
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FIGURE  4.      PROGRAM   LOADS  COORDINATE  DEFINITIONS. 
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1)  development  of  topographic  data  by  sampling  topographic  maps, 

2)  approximation  of  grade  lines  for  all  possible  highway  routes. 
Trend  surface  analysis  develops  a  series  of  smoothed  surfaces  of  varying 
complexity.  These  trend  surfaces  are  assumed  to  approximate  all  possible 
grade  lines  in  the  area  while  the  residuals,  measuring  the  lack  of  fit  of 
these  surfaces  to  the  elevation  data,  represent  the  magnitudes  of  cut  or 
fill  required  to  build  a  highway  having  a  grade-line  following  the  trend 
surfaces.  The  residuals  define  those  areas  lying  far  above  or  far  below 
the  local  mean  elevation.  Thus,  topographic  residuals  are  a  suitable  value 
of  the  earthwork  cost  factor  (5,  10). 

In  cases  where  the  topography  is  too  complex  to  be  satisfactorily 
explained  by  a  single  trend  surface,  it  is  possible  to  fit  trend  surfaces 
to  overlapping  sections  of  the  area.  After  such  "piece-wise"  fitting  of 
trend  surfaces,  a  trend  mosaic  can  be  assembled  from  the  various  parts,  much 
as  a  mosaic  is  assembled  from  individual  aerial  photographs. 

Program  TREND  is  used  in  the  GGARS  System  to  fit  trend  surfaces  for 
the  above  purposes.  Trends  up  to  the  fifth  degree  can  be  fitted.   The 
current  program  represents  the  latest  version  in  an  evolving  series  of  programs 
that  were  previously  called  "POLYFIT"  (h) . 

The  computation  of  the  trend  surface  coefficients  is  accomplished  by 
the  matrix  algebra  procedures,  as  described  in  Chapter  2  of  Draper  and 
Smith  (l).  The  procedure  involves  the  construction  of  several  matrices, 
resulting  in  the  construction  of: 
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(1)  an  (XY)  (XY)  matrix. 

i 

(2)  an  (XY)  Z  matrix. 

(3)  the  inverse  of  the  (XY)  (XY)  matrix  =  ((XY)  (XY))"1 
The  matrix  of  coefficients  =  A  is  computed  by  matrix  multiplication  as 

follows : 

A  =  ((XY)  (XY))_1(XY)  Z 

(XY)  is  a  matrix  of  independent  variables 

i 
(XY)  is  the  transpose  of  the  (XY)  matrix. 

Z  is  a  matrix  of  the  dependent  variable 

(XY)  (XY)  matrix  contains  the  sums  of  all 

the  powers  and  products  of  the  independent  variables 

(i.e.  X,  Y,  X2,  Y2,  XY,  etc). 

Figure  5  is  a  generalized  flow  chart  of  Program  TREND  and  its  associated 
subroutines  to  show  the  sequence  of  operations.  The  matrix  inversion  process 
performed  by  Subroutine  INVER  is  crucial  to  the  accuracy  of  the  method.   One 
way  of  checking  on  the  accuracy  of  the  computed  inverse  is  to  multiply  the 
inverse  by  the  original  matrix.  The  result  should  be  an  identity  matrix.   If 
it  is  not  the  computed  coefficients  should  be  suspected  as  being  erroreous . 
This  check  may  be  requested  as  an  option. 

Program  TREND  has  been  developed  so  that  all  input  functions  are 
contained  within  Subroutine  IREAD  while  most  of  the  output  functions  are 
within  Subroutine  RITER.  Subroutine  POLY  evaluates  the  coefficients  of  any 
surface  at  any  X  and  Y  coordinates  to  compute  the  value  of  the  surface.   It 
is  used  to  compute  values  of  the  surfaces  on  a  grid  network  to  prepare 
for  the  production  of  contour  maps  of  the  surfaces. 

Subroutine  ANOVA  computes  the  various  analysis  of  variance  statistics 
and  displays  them.  Methods  of  utilizing  these  statistics  are  described 
by  Turner  (5). 
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FIGURE    5  .    FLOWCHART  FOR  PROGRAM  TREND. 
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Program  TREND  utilizes  Subroutines  PRTCONU  and  STATS  to  produce  contour 
maps  of  the  surfaces  on  the  printer.  These  subroutines  are  described  in 
Part  C  of  this  manual  (8). 

A  large  number  of  input  and  output  options  are  available  during  the 
use  of  Program  TREND.  These  may  be  requested  through  the  use  of  appropriate 
control  cards  as  described  in  the  documentation  of  Subroutine  IREAD  which 
utilizes  them  (see  pages  A-59  to  A-67).  An  example  of  a  job  submission  deck 
is  shown  in  Figure  8,  sample  input  data  is  shown  on  page  A-85,  and  the 
resulting  sample  output  on  pages  A-86  to  A-96. 
d)  Program  RESMAP 

Program  RESMAP  constructs  residual  maps  showing  the  lack  of  fit  of  the 
various  degrees  of  trend  surface  to  the  original  data.  Gridded  value  arrays 
for  the  original  data  and  for  each  of  the  trend  surfaces  are  supplied. 
Since  the  original  data  normally  consists  of  a  series  of  irregularly  spaced 
X,  Y,  and  Z-coordinate  observations,  these  must  be  converted  to  a  set  of 
grid  values.  This  is  accomplished  by  the  use  of  Program  MAPPER  and  Subroutine 
GRID  as  described  in  Part  C  of  this  Manual  (8).  Gridded  value  arrays  for 
each  trend  surface  are  supplied  from  Program  TREND. 

RESMAP  takes  each  trend  surface  array  in  turn  and  computes  the  residual 
value  for  each  point  in  the  array  as  follows : 

Residual  (l,J)  =  Interpolated  "Original"  Value  (K,j)  -  Trend 

Surface  Value  (l,J) 
where  I  and  J  are  the  row  and  column  coordinates. 

These  residual  values  are  then  contoured  by  PRTCONU.  An  option  allows 
these  residual  arrays  to  be  punched  on  cards  or  stored  on  tape. 
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e)  Program  LOADG 

This  program  will  load  gridded  data  arrays  onto  tape  in  a  standard  GGARS 
format. 

Each  array  will  appear  on  the  tape  as  a  series  of  binary  records.  The 
first  record  contains  identification  information,  as  follows : 

FACTOR  -  Uo  characters  identifying  the  factor. 
AREA    -  30  characters  identifying  the  area. 
NCOL    -  an  integer  giving  the  number  of  columns 
ITYPE   -  an  integer  identifying  the  type  of  array  data. 
Each  row  in  the  array  forms  a  subsequent  record.  The  array  is  thus  stored 
row  wise  from  the  top  down. 

A  single  array,  or  a  group  of  several  arrays,  separated  from  subsequent 
arrays  by  an  End-of -Record  card,  will  form  a  single  binary  file  on  the  tape 
(i.e.  when  an  End-of -Record  card  is  encountered  an  End-of -File  is  written 
on  the  tape). 

Such  binary  files  can  contain  factor  combinations  which  are  often 
analyzed  together  by  Program  NETSUM  (see  Part  B  of  this  Manual). 

LCADG  can  load  any  number  of  binary  files  onto  tape  at  any  one  time. 

PROGRAM  OPERATING  CHARACTERISTICS 
a)  Program  VECTOR 

Program  VECTOR  will  load,  compile,  and  execute  in  approximately  60K 
(octal)  central  memory  locations.   Each  analysis  requires  about  7  seconds. 
Figure  6  shows  the  sequence  of  control  cards  required  to  perform  an 
analysis .  Examples  of  typical  input  and  output  are  shown  on  pages  A-  39 
to  A- 1*1. 


*0n  the  CDC  65OO. 
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b)  Program  LOADS 

Program  LOADS  requires  UOK  (octal)  central  memory  locations  to  load, 
compile,  and  execute.  Execution  requires  23. 5K  locations.   It  loaded 
approximately  2500  observations  onto  tape  in  just  under  13  seconds.  Figure 
7  shows  a  typical  deck  setup. 

c)  Program  TREND 

Program  TREND  requires  60K  (octal)  central  memory  locations  to  load, 
compile  and  execute.  Execution  requires  U7-5K  locations.  Running  times 
depend  on  the  number  of  observations,  the  number  of  trend  surfaces  filled, 
and  the  requested  options.  Program  TREND  computed  all  five  degrees  of 
trend  surfaces  for  just  over  2500  observations  in  U5  seconds  (including  load 
and  compile  times).  This  included  the  production  of  printer  maps  for  all 
surfaces. 

Figure  8  shows  a  typical  deck  setup  for  card  data.   If  data  is  supplied 
from  tape,  the  setup  is  similar  except  that  the  data  cards  are  omitted. 
There  are  many  options  which  are  described  in  detail  in  the  documentation 
of  Subroutine  IREAD,  pages  A-59  to  A-66.  A  sample  of  input  data  and  the 
resulting  output  is  shown  on  pages  A-85  to  A-96. 

d)  Program  RESMAP 

Program  RESMAP  will  load,  compile,  and  execute  in  50K  (octal)  central 
memory  locations.  Running  times  are  dependent  on  the  size  of  the  arrays 
and  the  number  of  maps  produced.  The  time  can  be  estimated  by  allowing 
approximately  3  seconds  for  each  copy  of  each  map  produced  and  adding  another 
3  seconds  for  load  and  compile  times.   Figure  9  shows  a  typical  deck  setup. 
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FIGURE  7.   PROGRAM  LOADS-   TYPICAL   DECK   SETUP. 
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FIGURE   8.   PROGRAM    TREND -TYPICAL    DECK    SETUP. 
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FIGURE  9.     PROGRAM  RESMAP-   TYPICAL    DECK    SETUP 


A-22 


e)  Program  LQADG 

Program  LQADG  requires  UOK  (octal)  central  memory  locations  to  load, 
compile  and  execute.  Running  times  are  dependent  on  the  number  of  arrays 
being  loaded  and  their  sizes.  An  allowance  of  5  seconds  for  each  array 
should  be  adequate.  Figure  10  shows  a  typical  deck  setup. 

PROGRAM  LISTINGS  AMD  DOCUMENTATION 
The  following  pages  contain  listings  of  the  various  routines  in 
order.  All  are  written  in  basic  FORTRAN  IV.  Each  routine  is  documented 
through  the  use  of  appropriate  comment  cards  at  the  beginning  and  throughout 
the  program. 

Each  program  has  been  analyzed  by  Program  TIDY,  available  as  a 
library  program  at  Purdue,  which  resequences  all  statement  numbers,  removes 
extraneous  numbers,  and  generally  makes  the  program  neater.  For  this 
reason  the  statement  numbers  in  each  routine  increase  progressively  in 
increments  of  five,  and  all  FORMATS  (if  any)  are  collected  at  the  end. 
The  cards  are  each  given  a  sequence  number  which  appears  at  the  far  right. 
a)  The  Use  of  FORSTAT  Data 

Each  routine  was  next  analyzed  by  another  program,  FORSTAT.  This 
checks  the  syntax  of  source  (FORTRAN)  programs  and  produces  cross-reference 
listings  of  all  variables  and  statement  numbers  found  within  the  routine. 
The  cross-reference  listings  utilize  a  series  of  "FORSTAT  Numbers"  which 
are  found  at  the  far  left  of  the  listing.  A  vertical  ink  margin  line  has 
been  added  to  separate  these  numbers  from  the  rest  of  the  listing.   The 
cross-reference  listings  are  found  at  the  end  of  each  routine. 
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FIGURE   10.  PROGRAM    LOADG  -TYPICAL   DECK   SETUP. 
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All  symbols  defined  in  the  routine  are  listed  alphabetically  and  are 
identified  according  to  type,   definition  and  references    (where  referenced) 
The  following  symbols  are  utilized  in  the  cross-reference  listings :- 
Types  I  =  INTEGER 

R  =  REAL 

D  =  DOUBLE  PRECISION 
C  =  COMPLEX 
L  =  LOGICAL 
N  =  NAME  LIST 
ASF  =  ARITHMETIC 
Definition  C  -  COMMON 

D  =  DIMENSION 
E  =  EQUIVALENCE 

All  other  terms  are  believed  self-explanatory.      These  cross- 
reference  listings  will  be  most  useful  should  any  reader  desire  to  make 
changes  to  these  routines. 


A-25 


PROGRAM 
LISTINGS 

FOR 
CDC  6500 


A-26 


b)    LISTING     OF    PROGRAM     VECTOR 

PKO<tk4M    VECTOR  (  ThiPill  »  OUTPUT  »PtlF'CF4,1  APF.«i  =  INPUT  tTAPEfsrllTPUT) 

C 

C  SO'JKOL  r  T  "i£     v/tCTnw 

c 


To  vmALY/F  TFPRAlN  Tfi  OF IfcRMlME  SOMFACF  ROUGHNESS  STaTISITC 
tuO  TFCHNTOUFS  APF  USUI— ]  )  FlcVATlONS  /I  Hfl-  A  Nil  L  Y  7FO  FOP  Mt"AN  HANG' 
VAPIAhCF  ,  STANDARD  OFVIAT1CN.  tie.  2>  GROUPS  (IF  rnP.EE  POINTS  APF 
Tn  OM  T''F  PLAM-S  WHOSp  OR.TF  NT  a  T  [On  IS  nF  TfRmT  MF.O.  t-'F  aS'.iPF_<=  OF  VF  C 
STRFNUTH  AND  nlspFPSlOU  A.'F  USt|)  TO  KsTTMuTF  RUUKMF'FSS.  TvPF  1 
AMD  TYPF  ?  1pIanO|FS  i'«i«  Rf  FTlTFn.  a  vapIRIIITY  faCIOp,  V,  IS 
r«irui  iTF0   from    Thf    Raimgf:    of    flf  vaT  iok  «;/Ti  if-    Rouo,hnfss   "facijbf,    k. 


riy    - 
HOU(?HmFSS    In 


THF    PRnQRA*     IS    PaPTLY    riASEO    ON    A    PPor.wAM    iiFSOMhrr 

h,      n,      UOHSONt      lVr,/, 

FOPlpAn      IV     PHOfiP.lMS     TO     DFTFOH'F      SIIHFAOF 

IopocRaPhY  m)P  l  hK  cue  340n  COnPUTpR, 

COmPUTfH  CUuTRThlTiun  |4i  kaNSas  GFi>l_Oo,TOAl  SUl'VFY 

l.A'vPFNrF.  KANSAS. 


Thf    G(  APS    STSTr.M 
PUHUUS     UNIVERSITY 


«.     K,     TupivFq* 
CIVTI      E.nO,|nfFR  InR, 
Iiji-JF     lVr-w. 


CAPO  OROFR  IS 

A)   CA>'OS  JHICH  APPtAR  ONCF.  — - 
1  )  Ei. 11  OF  pFOopO  (7-H-9) 
■>)     (  1  n  F  Capo  (IiT|)   -(COLS  1-nn) 

A  SFN'FRaI.  IT  TIF.  fhlrn  WILL  APPEAR  A1- 
OF  FACH  OAl a  SET . 


A  HKAOTF'P  AT  ThF  HF 


R) 


CA"US  uH|(-h 
I  b]7F/sCALF 
"-{001. Si 
II-  (001  S4 
KSIFP- (0 
YSTFP-IO 
hnn-  ( CO 
I  0"TPUT  OPT 
TPlr\T-(c 
IPLOl- (r 
I S<lP- (0 

iHjsl- (r 
ITYPF.-IO 


aRc- 
CAp 

-3)  . 

-ft)  . 

01  S 
"US 
LS    i 

IOM 

01  1 
oi  1 
01    s 

OLS 

"I-     1 


pe'jui 

n    (  m  <  n 

IS     Tmf 

is  Tne 
1  1-20) 
21-30) 

1  -411  )   • 
OAK!) 
1  .IF     I 
)  .  If     I 
I, IF     1 

ld-1 1 ) 
?)     . 


RFC 


'  S  1  F  P 


,,H|* 

.IS 

.  ts 

1  s 

(id 
pi  0 1 
Pi  di 
Sf  JR 

I 


(IF    This 

uMSLFSiH 
STFOFOGRA 

■>)      Fill*  MAT     CAP 
SHFOTF IF 

<,  )     PHOnLFM    CA 

a   hhf 

c)      SFPJES     OF 


CAP 
LIST 

M     Op 

0     (F 
S    Fo 

»r    ( 

WHlO 
DA  I  A 


0  lb    H 

1  NO,     (IF 
Hlbl  0 

mT)        - 

PF-A  1     0 

ppUp) 

"ILL 

CAPOS 


IF 

IF  I 

IF  1 

IF  I 

IF  I 

Lank 

ALL 

3UAW 

(0(iL 

■      I'A 

-(C 
SF> 
O  1  V 


ms 
NTS 
TyP 
1  YP 
Typ 
S  I 
T* 
S.) 

S     i 
Ta 

01  s 

VF 


pap 

YST 
OF 
OF 
ST 
ST 

i»UT 

IS* 

STp 

STf 

A      I 

IS 

1(1 

1 1  ■> 
f  =  i 

F  -■? 
F  =1 
AND 

1  «w 

-Rnl 


T    OF    Faoh    oaIa    «fT.    --- 
F  P  »  n  I B  N ) 

ROWS. 

C0|  HF'NC. 
/F     nF    C o (_  1 1  ►.■) p- .    SPACING     IN    F|.Ev 

n    of  Pnv  spacing  in  flfvaTi 

H      [N     pFr.UFFS     OF     Y     A»  I  c  . 
IPl!HlSni),lHlST(?),[TYPF) 


i*HK.1(H  OF  M.iPMa(  VFCl 
]  SI  yigC  iif  TPi  ANo,|  F  Sip 
Supopf  s<~FO. 

)  =)   ,  A     HI  Si  0(.*4M     OF      0  IF 

)  =  1  ,  .1  HISTOGRAM  (,f  STp 
1  TyPF  )  1  PI  anG|  F'S  j,wf 
»  TYPF  ?  TP1  AMfH  F  S  APF 
,TYPF«  1  AND  ?  TPl«Mi| 
ARO  OUTPUT  IS  PROOUCFIr 
'.I  F     SlPi  KF  S    afjO    0IPS»  A 


0"S     IS 
IkFS    A 

S      IS     F' 

I<FS  I 
FlTTf 
FlTTf 
Fs  HPF 
.-  — Ty 
MO    NO 


SFT    WHICH    Fni  LO»S 

1-HO) 
TO     lOF'-TIFY     THF     OATa     SF  T , 

FLFVATICN    "aTpIX. 


Vt'C 
v/FO 

S.         VFO 

F  VFO 

hSt>'FO 

TOR  VFC 
VFC 
VFC 
VFC 
VF.C 

VFO 

VI  C 
Vf  r 
>/FC 
VKC 
VFC 
VFC 
VFC 
"FC 

VFC 

VFC 
VFO 
VFC 
VFC 
VFC 
vFC 
VFC 
VFC 

---V/FC 

VFC 
VFC 
VFC 
vf  C 
VFC 
OlFJMMOVFC 
VFC 

\^r 

VFC 

VFC 
VFC 
V^ 
ATTON  UVFC 
ON  HNlTvFC 
VFC 
VF  C 
VFC 

PRoriUc«Fc 
MD    OIPSVEC 

Vt  C 

LOTTED. VCC 

S    P|  OTTVFC 
0.  Vfc  C 

11.  V  F  C 

FITTFOvFC 
PF     1     TRVF.C 
VF  C 

\Jl  r 
VFC 

«FC 
vF  r 
VFC 
tff  f 


111 
1  ? 
U 
1  6 
1« 


?ft 
?« 
in 

1« 


S.I 

6  1 

IS4 
ftft 

^^ 

70 
7? 

7^ 


RO 


PFOUIPFO  SUrPoUTxmfS  — 

iiLis  —  which  cai.culaTfS  strjkfs  amo  oips  fpoh  vfctor  cosp 

SlATS    —    WHICH    CALdHAlFb    STATTsTJCai     mFasUpfS. 

rfliLFF"    —    WHICH    OFNKpaIF.S    WU|.FF    MFT    STfRFuCRamS    fpcm    vfciors. 

HISTO    --    WHTCH    CFNFuAlFS    HlSlOOPAF-S    OF     SlRlKF     AM)    0TP    0IS1RT 

OPTIMAL    SURMnoTlNfS    -- 
MONF. 


VFO 

vFC 
VFC 

FS.VFC 
VF  0 
VFC 

RUT VFC 
VEC 
VFO 
VFC 
VFO 


Rft 

QR 

100 

102 

I  n* 

I  Oft 
1  Ofl 

1  in 
1  1  ? 
1  1  <> 
1  ift 
If 
l?n 
i  ?2 

I  ?4 
l?ft 
1?H 
1  31 
1  3? 
1  34 
13ft 
13R 
1  wl 
14? 

I   4<- 

1  Aft 

1   U  - 

I  SO 

is? 

154 

'  Sft 


A-27 


in. 
1  1. 
1?. 
|3. 

14. 
IS. 
1  ft, 


1  7, 
Ifl. 
19, 
20. 
21  . 
22. 
?3, 
?4, 
?s. 
?fi. 
?7, 


i?9, 
30, 
31. 
32. 
33, 


14. 
IS. 
36. 
37. 
3«. 
39. 

40. 
41  . 
4?. 
41. 
44. 

45. 
46. 

47. 
4«. 

49. 

so. 
51. 
52. 

ST. 

S4. 

55. 

H  7. 
ULt. 


10 


. \i>  r  l  sh 

v'Ff  IhO 

niMFNSTON  F(?n,?n)»  xnoh*  moo) ,  yndm^  (nnni .  ZNnPM(soo),  nPS(HOn).  vec  162 

lSTl((flOn)t     711L(«)«     MTlhli     PW0M(8)                                                                                               VFr  164 

nlMFN^TON      TH1STI?)                                                                                                                                                                     VFC  1*6 

I/KC  lh» 

QcMi    r  t tlf,*/\<;tf  p. *imD  daTa  foh^.at   cflMD1;                                                              vfer  1 70 

HE..D     (S,t.n)      T  1  T I                                                                                                                                                    Vfr  174 

RE/W1     n,m     M.M.XSTFH.YSlEPtDlFtN                                                                                                    i/F  C  176 

READ     [c, .  1  (n  )     ]pi  nl  .  ISKTH,  (  IhlST  (  I  )  .  T=  1  ,?]  .  T  I  Yf>f                                                     \)fC  I7A 

ptiMi     <«i»90)     "'T                                                                                                                                            VET  lhi 

— "                                                                                   yEC  1P^ 

WPTT6     ni  i  V    HFAMInHc;                                                                                                                                                     VET  1R4 

ver  inf. 

w H  I  T I      (  6  •  '  o s )                                                                                                                                                  ijFf  1»« 

"Kl  TF     IfiillOl                                                                                                                                                            VEC  190 

WHITF     (fit  115)     -r  T  TL                                                                                                                                    VEC  I"? 

WkITF     (6,110)                                                                                                                                                  vFC  144 

WKJTF     (fi,]?0)     F»,M»XS7EP.YSTEP.nIRN,FMT                                                                              Vtf  Ifs 

Rt'An    <  s  .  ]  l  ^  )    Pnnn                                                                                                                              vFC  I'Ju 
w «  I  I  F     (  h  .  U  n  >                                                                                                                                 vi-r^nn 

whitf    (b«i?m    ppob                                                                                                   vfcr  ?o? 

WHITE     (6,110)                                                                                                                                                            UFr  204 

^J^C  ?06 

V  (■  (  ?  0  H 

i/Ef  ?in 

yec  212 

vFc  ?i* 

vft  ?]6 


QFM1     ANn     WHITF     FLFV'aTJON     '-"AlHT* 


?o 
c 

r 


tiO     I  o     1  =  1,1. 

HE  ail     C>.F'in      (F  ( I  .J)  •  Jul  »N) 

w H I T F      (6.1301 

no   15    i=l,w 

•IN  r  IF       (ft.l3<=)       (F    I  I  ,11  ,  J=l  ,N) 
WH1TF"      (6(140) 

TAIL    STaIS     (M»M.?0.20»F.EHAH,EVAHtESIG»EMAX.EMTM) 
WHITE    (  e. ,  i  i  n  ) 

UK  I TF  (6.145) 

WK 1 7F  (1,140) 

WHITF  16, 175)      FHAH,fV4H,EStlJ.FMAXtFMlM 

hm  I  TF  (6,1  |0) 

RFMFHeTF    MOHMflLS     m    ThIaNHLES 
TTYPf     nFTFF.»TMFS     TYF'F     OF     THIaNGLF    HElNfi    FTTTFO 

MM=M-1 

NN=N-1 

IF     (TTYPF.I  T.1.0P.ITYPF.GT.3)     TlYPF=l 

fiO    Tn    (20 ,30, Hn) »     1 TTHF 

K«] 

FIT     1  YPF     1     Th  T  ANfll  F  S 

(FIHST     rulAMPLF     IS    UtFINFii    HY    P051TIOMS     (  1  ,  1  ),(  1  ,?).(?,  1  )  ) 


NNK=1 

WHITF  (6.150)  NNK 
wHl TF  (6,1101 
01)  ?5  |:l,"u 

T  T  =  T  ♦  1 
no  25   1=1  , KM 

1.1=  I*  I 

A  =  (F  (  I  »J  I)  -F  (  I  ,  Jl  )  "YSlF.P 
H=(F(1 (J)-E(IT(J) )*XS1EP 

r=x«;  (fp«ysTf° 

OIVaSORT  (A»»?*B»*<?*r.«*2) 

XWOKM  (K  I  =A/i'TV 

YmODM ( K ) =H/nTV 

7M0KM ( K ) =r/p  1  V 

K  =  K*  1 

"=  (E(  !  T  ,.IJ)  -F  <  IT  ,J>  I »YSTfcP 

5=(F(!(JJ)-FUl(JJ)  )«XST£P 

r  r  X  q  T  F  P  »  v  s  T  F  p 

nitfsSOHT (A»»?*H»«2*r«»2) 

XMORM (K) sA/OIV 

YMOPM (K) =M/nTv 

/MOUM(K)  =  (Vol  V 


v£( 

?1H 

VFT 

??n 

VE( 

^?j 

vFr 

2?  4 

VFC 

?2ft 

VEC 

??H 

VFC 

?10 

vk: 

?  i> 

VFC 

234 

vFr 

2  ih 

vfc 

234 

WE( 

?40 

vEr 

?«2 

VFC 

?44 

VFC 

?4., 

VfcC 

".« 

VFC 

^so 

i/FC 

?s^ 

VFf 

?S4 

>/Fr 

?Sft 

VFC 

P'iH 

yEC 

?in 

VEC 

?62 

VFC 

264 

>;c  ( 

-»'-'. 

>'EC 

?*,« 

wft 

?  rn 

vfr 

p  7^ 

"FT 

^74 

VEC 

?7h 

«M 

2  r« 

>.F( 

'no 

wF  f 

^>ap 

VF( 

?«4 

VF( 

?flft 

VFT 

?«H 

v/F( 

?«0 

wf  r 

?y^ 

vf  r 

?q* 

wt-r 

pwr, 

wf  r 

^u 

>lFf 

ion 

VF.C 

(0? 

vtr 

?04 

vFr 

10h 

l/Ff 

1(M 

A-28 


60. 

m  . 

<•?. 

63. 
64, 
65, 
66, 
fi7. 
fifl. 
69. 
70. 

71  . 

72. 
73. 

75, 
76. 

77. 
78. 

74. 
BO, 
HI. 

R2, 

B  1. 

04. 


P5. 

4S 

P6. 

R7. 

SO 

PH. 

PP. 

90. 

91. 

«?. 

"3. 

94. 

«=;. 

ss 

P6, 

"57. 

60 

PH. 

Inn. 
ldl . 
In?. 
In3. 
In4, 
Ins. 

1  Oh, 

ln7. 


InH, 
Inf. 
110. 
111. 
11?. 

m. 


114. 
IIS. 

1  is, 


P  II      I  YPF    ?     TUTaNbI.eS 

ifjw-.T    TWlANGLF     TS    OtHufi'    HV    POSITIONS     1 1  .1  )•  (2.1  )»<?•?)  ) 


is 


NN*  =  ? 

w«iTF  (Sil^n)    ^:^iK 

tertl  TF  ( s. 1 1 0) 

00     TS  T=  1  ."'' 

I  T=T*1 

00    3S     1=1  ,  rM 
JJ=J*1 

a=  if  i  ii  ,.D  -f  ( 1 1 ,  jj)  )»ystfp 

«= (F i I , j) -E ( I  I » J) ) »XS(Ep 

r=XsIEP«YS IFP 

DItfsSOPT (A**?*b*»2*r»«2) 

tUDO*'  (K  )  =  A/HT  v 
YNOHM(K) =H/nI V 

7mo«m (k ) =r/ni v 

k=k*1 

»  = (F ( 1 , j) -F ( I. JJ) ) OYSIE" 

H  =  (Fill  ..IJ)  -F  (  I  .  JJ)  )«<STtP 

C=xsTFP»YSTf  P 

OTVrSOPT ( A*«?*B«»2  +  C»«2) 

XMO»M (K) =A/nT V 

YN09M  (K)  =H/OT V 

/now"  ( k) =r/pxv 


MK: 


IF      (TSk-JP-l)      4S.S0.4S 


C  rcl CJLATF  strike  ami  dip  OF  ail  triangles 

C       anT     !»■     1S*IP=1     00    (JOT    LIST    VA|UfS 


*H  ITF      (S,  ]SS) 

wPI  TF     (6« 1 60) 

CONTINUE 

00     SS     Kb  1  ,NK 

A|   PHA  =  XM1P(<  (K) 

HFTAsYMOSM (K ) 

OAMMA=7K0Ph  (  K) 

CALI      IIITS     (Al  PHA.BE  TA  .C.AMMa  . ST P 1 K E , 0  I P . 0  I PN ) 

nps(K) =  Mp 

STq(K) bsThIkF 
COimTTNOF 

IE     ([SKlP-1)     6n,6S,*-0 

rtKlTF     (6.16S)      (KiSTf'(K)  .DPS  (K)  iKsl.NK) 
WHI  TF     (fulini 
rO'iTTMHF 


C       rri^PUfF    viFAN    amo    sTaNoaHO    i)EVt«TionS    OF    STRIKFS    amo    OIR« 


TALI.    STaIS     MtM<.l.H0n.SlH,blp.nAR,eTPV<\w,ST*SIO.STpMAX,STRMtM 

CAl  L     STAlS     (].NK,l.HU0,OPS.uPSHAH,OPSV.iR,0PbSlfi.PPSMAX,npSMlN) 

WHITE     (6,170) 

WHITF     (M17S)     <;TpbAH«STi<VARiSTRSI<>,STR'-aX,STPMf' 

WKTTF     (S.laO) 

WH1TF     <  *> ,  1  7S  )     pPSbflH,0CSVAO,0PSSTf>,0PSi-'AX.nPSMI(>' 

WH  1  TF      (M,  1  1  r\) 

1  =  1 


<-  A I  L  HISTO  FOP  HTSTOOHnMS  CF  STPJKES  ANO  OIPS 

IF  (  IHIS1  ( I )  ,Nf  .1 )  GO  ro  7n 

f'ALI   HISlO  (MK.OPSMAX.n.OinPS.PHOH,  T  .N>'K) 

T  =  1  »  1 

IF   (IHIST (I) ,NF.l)  00  Td  7=. 

C<M.L  HISIn  (NK.STHMAX.SThMINtSTK.PPOH , T ,MMK) 

COOT  T  Nl« '  F 


r    COMPOTE  VECTOR  STRFnPTh 


Torx=o.o 

TO  I  YsO.il 

TOT/aO.O 


VEC  3'4 

»FT  316 

VFr  319 

vtc  3?n 

«<F-r  \?? 

VEC  3?4 

"EC  3?h 

vEC  3?« 

VEC  330 

VfC  33? 

VFC  ^34 

VEC  336 

VEC  33" 

vEC  340 

VEC  34? 

VfC  344 

V?C  346 

VEC  34R 

\JFC  3SH 

VFC  3S? 

vEC  3S4 

VEC  3S6 

VEC  3SM 

VEC  3*0 

VEC  36? 

vEC  164 

VEC  366 

VFC  3S« 

VFC  370 

VEC  37? 

vFC  374 

VEC  376 

VFC  37H 

\lfC  3R0 

VEC  3H? 

VEC  3«4 

VFC  3«6 

VEC  38H 

VEC  3P0 

VEC  39? 

VEC  3P4 

VFC  IPS 

vEC  3<J" 

VEC  400 

VFf  40? 

VFC  404 

\J^c  4  06 

VEC  40H 

VFC  410 

VE(  41? 

VFC  41 4 

VFC  41 S 

VEC  41 H 

vFC  4?u 

VEC  4?? 

VEC  4?4 

\I(C  4?6 

VFC  4?" 

VFC  430 

vFC  <»3? 

VfC  43* 

VEC  436 

VFf  43* 

vFr  44 n 

vEC  44? 

VfC  4  44 

VFC  44S 

yF.C  44" 

VEC  4S0 

\/FC  4S? 

VFC  4S4 

VSC  4S6 

vfc  4S^ 

VFC  4  fc  n 

vFC  4h? 


A-29 


117, 

\\». 
119. 
120. 

121. 


122. 
123. 
12*. 
125, 
126. 
127. 
12H. 
129. 
130, 
131. 


132. 
133. 
134. 
135. 
136. 
137. 
13A. 


Ill, 
HO, 
141. 
142, 


143, 
144, 
145. 
146, 
147. 

14P. 
149. 


150. 
151. 
152. 
153. 
154. 
155. 
156, 
157. 
]Sflt 
159. 
160. 

161. 
162. 
163. 

164. 

16b. 

166. 

167. 

16B, 


no  nn  k  =  1 .nk 

TnlXs I OTX+XmOPMIK) 

TnTY  =  TOTY»YK'ORM  (K) 
TUT  7=TOTZ*7N0R" (K ) 
/'H  =  ToTX»lnTX*TnTY*Ti)Tr*T0T7»ToTZ 

VrcTOR  STHEHGTH  SmoWN  HY  FIlHFR  Ri  nH  R, DISPERSION  PY  F« 


«1  =  c;ort  </R) 

FNKsN* 

RsRl /FNK 

FSTK=(FNk-1.0)/(FNk-R1) 

WK1TF   (6, 1  AS) 

WHITF  (fi.UO! 

WH  ITF  (6,  iq(i)  Pl  ,H 

WHITF  (h.105)  F5TK 

Vs(fmax-FmTN)/aI  OGIO(ESTK) 

wrJTF  (6,2no)  y 


C  COM 

r  m 
c 


100 
105 

no 

]  1? 
120 


125 

1  30 
1  35 
140 
145 
1  50 
155 
1  h(i 
1  f-u 
1  70 
175 

1  po 
1  85 

1  tip 

i  t»5 

?nn 

205 

210 


PI)  r  E  VECTOR  "FAN  AND  MEAN  Pi  ANE  ORIENTATION 
RFCTTON  COSTMFS  OF  MEAN  VECTOR 

AI_Rha  =  toTX/r1 

BETlaTOTY/R] 

RAmM«sTOI7/R1 

WHITF  (6,215)  A|  PHA  .HETa,GA»iMA 

CALL  TIL'S  (ALPHA  ,RE  T  A  .(jAMm  A  ,  STR  I  KF  ,D  I  P  ,0  I  RN  ) 

WHtTF  (6,210)  STPIKE.DTP 

wRl  TF  (6,lln) 

HFN  TPLOl  =  l,  CAI  L  WULFF  TO  RLOl  A  STfKE^GRAm  OF  THE  IP 


IK  (  I  PI  OT.NE.l  )  (iO  TO  M5 

CALL  mil  It-f-  F  (NK,  XNORM,  YhORM, ?NOPM , PROP  .NNK.  ) 

IF   (  1  TYPE  .EQ.3.AND.NMK..E0.1  )  GO  10  30 

CO    TO  5 


FOR" 
FOR" 
FOR" 
FORM 
FOR 

1  ) 
FOR 

F  OHM 
1,13, 

2  U«T 
FORM 
FOR" 
FOhM 
FOR 
FORM 
FORM 
FO 

FORM 
F  ORm 
FORM 
FOR" 

1  I  AT  T 
FOR" 
FOR" 
FORM 

1  ,F 
FORM 

1  cTor 
FORM 

IK)  , 
FORM 
I  ,  F  I  '1 

F  OHM 

I  ST  it 

FM' 


AT  (EA1 0) 

»T  (2T3,4X,3FlO.0) 

A!  (  1 1  .3X . T 1 .4>  .21  )  .  1  1) 

AT  (lHl  »]  fix  .30HCOMPIITAT1ON5  HY  PROGRAM  VECTOR) 

at  (1ho,4  9h---- — ..-.-...._..- 


AT  (1H 

AT  ( 1H0 

7H  XSTF 

A  FORMA 

AT   ( 1H0 

(  1H0 

(1H 

(1HP 

(1H0 

(1H1 

(  1H0 

(  1H0 

(1H 

(  1H0 

(IN 

ONJ=iF  1  0 

AT  ( 1H0 

AT  ( 1 H 1 

a  1  (inn 

6,34H, 
AT  (  1H 
.  T5.F12 
AT  ( I H  0 

i<;  ,fi 

a  t  (  1  HO 
,S,7H 

AT   (  t  HO 

tkf  =  ,f>* 


HAl  0) 

]4hVAL0ES  PF.AO  1n//,12h  KO  OF  ROWSs, 13,12 

»,FH.2,7M  YSTEP=.E».2»  1«R  Y-AXIS  OTPFCUr 

=.12Ab) 

32PCOMPU1 AT  ION    NFSlll  TS    FOR    SllHAKEA     ,PA10) 

23HORlplNAI      ELFVaTJON    iiaTA//) 

12F10.2) 

3BHSTAUSTTCAL  MEASURES  OF  FLFVATIOm  OaTa 
32HTHE  PRDr-HAM  15  MOW  FITTING  TYPF   ,11,11 
3lF'SThIKE  AC I)  01P  OF  A|  L  T  R  T  ANGLES// ) 
lX,35HTRIAil<iLF  NUMrFR        STRTKF       DIP 

3  <  7  X  ,  1 3  «  1 1 X .»  E  H  .  4 ,  ?  X  » F  H  .  * )  ) 

44HS1  Al  ISTTOAL    NEASURES    OF     STRIKES    OF      Ml 
5HmEAN=,K10.4, 12H  VARlANrF=,F10.4,22H 

4>llH  MAX  1MJM=,F1  0.4,  1  ]  H  MINI  M(|M=  »F  1  0 

41H5TAT  I  ST  ICAL     MEASURES    OF    I1IPS    OF     TRlAf^ 

3bH"FASURFS    OF     VICTOR    STRE'fiTH    ARE     Pl  =  ,F'l 

HEpF    H=h1/M)vmpn    OF     TH1ANRI.FS.     ) 

5MHVKCTOR    nlSPFRSlON.MFASURFO    RY    FTSHEHS 

4,40H,       LOW    VA|  UKS     lNDICATF     VARUHLF     TF.HR 

60HTRE    VARl»bI|lTY    FACTOR,     V    =    ELEVATION 

.3) 

51HTHE    UIRFCtlniM    COSINFS    OF     T  HF    mfam    VtCT 

F  I  A  =  ,F  10. 5, MR       IjAMmAs.F1  0, 5) 

59HTHE    ORIFNfATlON    OF     THE     MFAN    PLAMF     IS    A 

4  »6H,     DIPs  ,F  M ,a  ) 


vFC 

*6  '• 

«K 

466 

^fc 

46fl 

vE  r 

4  r  0 

VF  ( 

472 

VEC 

4  74 

Tk                                  VFC 

4  7h 

VFC 

47H 

VEC 

4M0 

VEC 

4«2 

VEC 

!•«« 

VEC 

4H6 

VFC 

4M 

vK 

b<ia 

vF  r 

4Q2 

vFr 

4Q4 

vFC 

i.'J'i 

(/EC 

4-<»fl 

vi  r 

BOO 

ve  ( 

5^? 

Vf.  r 

Sl, 4 

VEC 

B06 

Vf  r 

so  9 

VEC 

si  II 

VEC 

s)  ? 

VFC 

S14 

VEC 

516 

VEC 

S1H 

VFC 

520 

VFC 

522 

IaNC-LES                   VFC 

524 

VFC 

S^h 

Vf  ( 

SPH 

VFC 

53n 

VEC 

si? 

Vf  f 

s  14 

VFC 

536 

VEC 

SIM 

Vt  C 

S4  n 

VEC 

Si? 

vec 

544 

VEC 

Mh 

54H 

sso 

VlC 

VEC 

SS? 

H    NO    OF    COLS  =  VF'C 

SS4 

Ns,FH.?,//13HVLC 

556 

VFC 

SSM 

VEC 

Sftn 

Vf  r 

562 

VEC 

SS4 

VEC 

566 

)                                   VEC 

SfcH 

H    TPTaNKLfS.) VFC 

s  70 

VEC 

S7? 

)                                     VFC 

574 

VFC 

s  76 

AmGLFS)                 VEC 

s  r« 

STANDARD    DF.WEC 

SM0 

.4)                               VFC 

SM? 

LFS)                                 VFC 

5  •<  ,x 

VF  C 

SHA 

2.4,7m,            Rsvfc 

S«M 

VEC 

SRO 

DISPERSION    FAVEC 

S^? 

ATN.)                         VEC 

S    <4 

RANGE/l  OGlFSTVFC 

sgf, 

vEC 

59M 

OR    ARE     ALPHAsvEC 

600 

VEC 

602 

S    FOLI  OttS. —    VF  C 

604 

'iv  r 

60CS 

vFf 

*-0« 

A-30 


P 

o 

UJ 
> 


/ 

i/ 

U 

UJ 

U. 

LJ 

<_> 

u 

- 

2 

2 

! 

U; 

jj 

.' 

r 

I 

r 

U 

U. 

LI 

u 

U- 

U 

U 

u 

u_ 

U 

- 

Q 

QC 

cr 

</ 

u. 

b. 

u. 

-1 

c 

C_ 

c 

u 

L 

CI 

Qt 

1 

1 

UJ 

UJ 

UJ 

a 

I 

5 

S 

i 

5 

1 

y 

*' 

Z 

* 

UJ 

I 

.u 

-. 

vl 

u 

: 

-J 

<i 

*- 

il 

<J 

¥- 

U 

1- 

<i 

c 

> 

1/ 

i 

h- 

C 

r\ 

<v 

(Y 

1/ 

-j 

■"' 

■i. 

~ 

< 

O 
O 

Q. 


O 


CO 

ID 


Ul 

o 

2 
UJ 

E 

U. 
UJ 

cc 

I 

CO 
CO 

o 
or 
o 


iiictiro'CtiiittJtiriairiiarctaiiiiifri 
oooooooooooooooooocoooooo 

U.U.U.U.U.UU.U.U.U.U.U.UU,    U.U.U-U.U.U.U.U-U-U.U. 


CO 

cc 
o 


uicfcif'Cifcircirci/circifcircifcifcir'Cif'Circiroifcifcifcif'circ 

HHHHHHHHiHHHHHHHHHitHr<4|\||V(\| 


A-31 


• 

• 

• 

■ 

• 

• 

• 

• 

• 

o 

,-. 

<r 

in 

i   iT 

a- 

pj 

i— « 

\S 

j 

f- 

j. 

h- 

a 

j. 

- 

-» 

.T 

i 

IV 

a 

o- 

If 

-3 

a 

a 

c* 

< 

r 

- 

-t 

<- 

P- 

L 

a 

c 

•3 

c\j 

r^ 

„• 

r^ 

r\j 

P» 

m 

m 

a- 

o  ^ 

«• 

<J 

m 

- 

>i 

•3 

■J. 

u 

a. 

cv 

f»i    C_" 

,""1 

r- 

^ 

— . 

o 

4 

~ 

4 

i 

a 

r*j 

r*-    - 

T 

■3 

m 

<\j 

- 

.£ 

* 

<C 

■» 

A 

fV 

rvj  o 

, 

. 

m 

a 

r 

o 

T 

cr 

C 

c: 

r» 

o 

^-.  r~ 

t\ 

m 

(\ 

C 

- 

lT 

<t 

J 

IT 

o 

a 

rv 

(\  a 

■f 

• 

C    c 

c 

rv 

ir 

J 

I 

a 

or 

cr 

>r  <v  <r 

CT 

>t 

c 

c 

r- 

r 

a  a 

tr 

c 

(\ 

n-  f 

p 

r\ 

tr 

C 

c 

*j 

o 

r* 

PI 

J 

j 

a 

a 

a. 

—  j- 

^ 

■I 

■*  ~~ 

rH 

~ 

rH 

"""' 

■"* 

^ 

' 

■*    r>  a  <   a   ■»  a1       -3  cccccr        (v  rv  rv  rv  r\j  (V       t\j  trc  r»  i--»n" 


a  iinaaair 


a.aaxaaaa.a.a.a.aaxa.act'-      __,-.-.___.- 


151        «  z 

c  a      t-      a  a  > 

_J  _l  Ul  HHH 


a  x  z  ir  a 

<  <  t-<  ►-  «i 

a  s  7  in  > 

V   U,  if  U    V   V, 

a  a  a  a  a  a 

a  c  c  o  c  o  u. 


■a  C 

q   x  z  e  ir  cr  ih 

«       «»->-(-«      l~     *       3     IT, 

u;uju.'UiUju.ik(SXH 


t-  i-  a  u.' 

m  c  «-a. 

•-  _j  sc  >- 

i  i-i  a  m>-       -5 

»-.-<»-..-.-,  -J  ->  * 


A-32 


~2  pj- 

C  (- 

—  ._> 

I-  C 

c  a 

z  x 

3  D 


a  a  a  a  ara  aaaaiacrcracicrcta'aa.cicrcicr 


it  in  k  z  o  a 
i/i      <  Si  4  —i—<  <  in  u.  x  o.  z  a.  z 

a  ki-      cr>-s3ul>i-_ixvr-j      u.  a  lu  a  tu  a 

c  a<oaoo:aaQ_i-i-(-i-       -j  o  ►-  c  i-  o 

or       i-*o>->->-»-»-i-»->-»-t-iCoc~       nyir^tn^a 
laai^i/iif.  i/iwi^uiunf,  ht-i-i-»->j«x>->-Ns 


1/ 

w 

1/ 

UJ 

UJ 

UJ 

u 

l_> 

u 

z 

^ 

2 

UJ 

UJ 

u. 

I 

X 

I 

u 

u. 

u 

u. 

u. 

U. 

u. 

u 

u 

a 

or 

a 

u 

u 

u 

c 

c 

t 

L 

LU 

J. 

a 

a 

s 

UJ 

u 

Ui 

a 

a 

a 

u 

Z 

z 

z 

u 

=> 

- 

c 

1 

Z 

r 

l/J 

1 

UJ 

n 

C 

I 

^ 

< 

.1 

a 

5 

x 

H 

U 

>- 

<1 

C 

J> 

tr 

1 

h- 

«I 

^ 

c 

C\J 

>f 

-£ 

m 

c 

[»1 

^> 

A-33 


d)    LISTING    OF    SUBROUTINE    TILTS 


SU-.POUTINF     TILTS     (  ALPHA  t  riK  Trt  .  GAMMA  .  ST  R  T  KF  » ll  I  P  <  n  T  RN) 
SUHROUTINF    TILTS 


PIJRPOSE-- 

Tn    lit  I  F  PMI  fJF    STRIKE    AMI    (J  I  P    OF    PLANES 


ThE     bOHS    SYSTfM 
PiiRHUE     UN!  VE  wc;  |  r  v 

»,      k  .       IU-I|\E'M 

CJVTI   KNMINFFKrl: 
HIIJF  |^f,1'. 


CALL  TTLTS(ALRHAiHFlA,GAMMA,STRIKE.iJIP,nIKM 

WHERE  -  ALPHA  =  DIRECTION  COSINF  OF  NORMAL  VFC TOR  To  X-AXIS 

RFTA   =  DIRECTION  COSINF  OF  MORMAl  VfCTOR  Tn  Y-A*TS 

GAMMA  =  DIRFCTION  COSINF  OF  NORMAL  VFCTow  To  7-AXIS 
STRIKK=  SlRTKE  nf    PI  ANfc'  (AzjmiiTh)  In  OE.GRFF.S 
DIP    =  OIP  OF  PLANF  IN  OEGWFFS 

OIRm   =  A/IMIJIH  (IN  DEGREES)  OF  Y-AXIS  OF  ELEVATION  GP1 


Pt'SCHIPTlON  OF  CONTROL  CAHUS- 
MOMF  RFOUIRFO 


RfOUIRFO  SUBROUTINES- 
MONF 


OPTIONAL  SUBROUTINES- 
NONE 


ANGG=ARCOS (GAMMA) 

AC  =  SORT(|AL_PHAJ»«2*  (HFTA)««^)«SlGN(l  .  0  .  ALPHA > *SIGN (  1  . O.HFTa) 

costh=heu/ac        XF        (  f)C.  £"Q  .  O- )      G>C    TO  /O 

■AHCOS (COS TH ) )»S7. 2957) +DIRN 

—    0>  o  ~nD  /G~  77 

wtr"PN  /O      STRiKK-    =0.O       774 


SfHTKF=( (] .S70B- 

(  c.  oir=anggJs7.?4sT"^/ 

~*'IRN  t- 

FNI) 


"<? 
J  4 

1* 

38 
40 
42 

4* 
46 
4M 

5n 
52 

54 
Sfi 

SA 
6,) 

*4 
66 

7T 
7? 
74 

Rfl 


73 


e)   FORSTAT    CROSS-REFERENCE    LISTINGS    FOR    SUBROUTINE    TILTS 


SUfWOUTINF 

TILTS 

label 

TYRE 

DFFINITION 

AC 

R 

4. 

5. 

ALPHA 

R 

1. 

4. 

A'^RG 

R 

3, 

7. 

ARfOS 

H 

1, 

FUNrTION 

6, 

BETA 

R 

1  . 

4. 

COSTH 

R 

5, 

f>, 

DTP 

H 

l  . 

7. 

DlRN 

R 

1  . 

6, 

G4MMA 

R 

1. 

3, 

SI  AN 

R 

4. 

FUNCTION 

SORT 

R 

4  , 

FUNf.T ION 

STRIKF 

R 

1  . 

6, 

SU«ROUT 

R 

?. 

REFERENCES 


13 

STMHOIS 

HEFTNED 

? 

MAXIMUM 

NUMRER  OF 

rfffHences 

1  1 

TOT/1 1. 

NUMRFR  of 

Rfffnfnces 

.  R 

AVFPAGE 

nIImhFR  OF 

RFFTKENCFS 
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f)    LISTING     OF    SUBROUTINE    WULFF 


7, 

B, 

9, 

in. 
11. 

1?. 
13. 

1*. 
1S. 

i  6, 

i  r. 


SUhROUTINF  *ULFF  (KOUNT2.  IX.TY.TZiPMOH.ITYPF ) 
SUBROUTINE  WULFF 
PURPOSE — 

TO  PLOT  A  SoaTTFR  DIAGRAM  OF  A  GROUP  OF  VECTORS  DEFINED  HY 
nIRECTlON  cosUEs.   THE  OIAGHAMS  MAY  HE  USFD  AS  OVFRl.  «YS  Oh    A  ?0- 
rfnti^ftfr  oiamftfr  kulff  net,  the  oip  anrlfs  have  ftefn  doubled  to 
increase  thf  oispfp^ion,  thf  number  in  fach  rectangle  imiicatfs  the 

MUMBKR  OF  POINTS  FALLING  IN  T^AT  RECTANGLE. 
REFERENCE--  SUhROUTInE   wulff  IN— 
T.  V.  LOUOOn,  lv64, 

computer  analvsis  of  orientation  data 

IN  sTRUC'UAl  GEOLOGY. 
UFFICF  UF  ^AVaL  RESEARCH,  i,FOr,,    HHANO 
TFCHNTCAL  REPORT  Hi  ONR  TASK  NO.  1P9-nb, 

129  PAGES. 


The  gcahs  system 
Purdue  university 


a,  K ,  IURNEBi 
CIVTl  ENGINEERING, 
JUNF  1969. 


USAGE  — 

CALI.  Wlil.FF  (KOUnT2,TX»TY.TZ.PHOR,ITYPF) 

WHERE-KPUNT7  =  NUMBFH  OF  TRIANGULAR  PLANES 

TX  =  DIRECTION  COSINES  OF  NORMAL  VFCfORS  TO 

TY  c  DIRECTION  COSINES  OF  NORMAL  VFClORC  TO 

T7  B  DIRECTION  COSINES  OF  NORMAL  VFCTOP?  TO 

pROH  b  0*Ta  SrT  f i T( e  (Up  To  flO  ALPHANUMERIC  CHaP*CTFRS 

1TYPF  =  TYPE  OF  TRIANGLES  BEING  ANALYZED 


■AXJ9 

Y-AXIS 

Z-AXT<; 


OESCHTPTlON  OF  CONTROL  CARDS. 
NONF  RFCUIHFD 


BFOUIREO  SUHHOUTlNfS- 

NONF 
OPTIONAL  SUPROUTlNfS- 

NONF 


DIMENSION  TR0W(4R»40)  ,  TX(qUO),  TY(HOO),  T7(S0O) 
OIMFNSION  PROB(H) 


ImTTIaLI^f.  PUTTING  -g  IN  THE  AREA  OUTSIOK  THF  STFPFOGRAM  , 
WRITF  (6ifl0) 


SFF  COMMENT  ON  CONSTANTS  hELOW 


P=l ,0/(3. 917*6,0) 
0=1  .0/(3. 917*5,0) 


FOR    POINTS    ON    THE    f.IHcUMFFRENct    OF    THF    PRTMITIVF    rlRCLf. 

(mp/?)««?*(no/?i««2  =1     .   ^  the  Trohfm  of  pyThagora? 


[10     IS    M=l  ,A« 

R=4HS(R.24,n) 

00    S    N=l ,40 

1R0W (M.N) eO 

M?  =  i;uRT((1.05..(R«p)«*2)/i1i<n>p)»0,=I 

N1  =  ?fl-M? 

Mpa?0*K|g 

00     1 0    N=l ,N1 

IR0W(M,N')a-q 
lid     IS    N=N2,A0 
IROw ( M,N) x-9 


10 

1? 

U 

16 
.   IP 

?n 
?? 
?A 
26 
?H 
30 
32 

14 
1ft 
1H 

40 

4? 
44 

4fi 

AH 
HO 
5? 
54 

Sft 

Hm 
M> 
62 
64 
66 
68 

7n 
?? 

74 
7ft 
7* 
Rfl 
R2 
P4 
rtft 
Rfl 
90 
9? 

9  A 
9ft 

j  on 
102 

104 
1  Oft 

I  OH 

I  in 
112 

114 
1  1ft 
IIR 

120 

122 

1?4 
1?6 
1?R 

no 
iip 

1  3A 

1  ift 

1  1" 

140 
14? 
144 

146 
l  4  q 
|«ifl 
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19. 
?0. 

2) . 

22. 
?3. 


24. 
S5, 


26. 

?7. 
?«. 
?9. 

30„ 

31. 

32. 
33. 
34. 

35. 
36. 
37. 

38. 


39. 

40. 
41. 


4?. 
43. 
44. 
45. 
46. 
47. 
«H, 
49. 
50, 
si. 


55. 
56, 

57, 

5B, 

*9. 

60. 
61. 

6?. 

63. 


rOMIJTf  THE  NllMpF*  OF  POINTb  Jm  EACH  HFCTAMil.E  OF  THF  «TERF0GR»M. 


DO  40  T=l ,KnuNT? 

I iP  =  T2  (  I) 

IF  (UP)  20,?5,2b 
20        TX(T)«-TX(I) 

TY(I ) =-TY(I) 
c 

C      IN  (HF  NEM  TWO  EXPRESSIONS,  3.C37  TS  ThF  RAUIUS  OF  THE  NET  IM  INCH 
C     6.0  IS  THE  MJmBfR  OF  LINE  OF  OUTPUT  PFR  TNCMi  5,0  IS  THF  NUMHFR  OF 
C     TwO-OTGTT  CHARACTERS  PEP  Itfn.  0.5  TS  A  ROUNUINC,  FACTrR. 
r 

25        M=((TY(T)/(1.0*ABS(T7(I))))*3.93  7«(-6.0))*1.0*3.93  7«6.0 
N=((TX(T)/(1.0*AHS(TZ(I))))«3.937«(-5.0))*1.0*3.9?7«5,0 

C 

C     OflUnLF  THE  OISPFRSlON 

C — 

c 

IF  (M.RF.25)  PO  TO  30 

IF  (N.LF.20)  N=2»N-?0 

TF  (N.fiF.21)  N=2»N-?1 

r,0  TO  35 
C 
30        m=2«M-2S 

IF  (N.LF.20)  N=2«N-?0 

IF  (N.GF.21)  N=2»N-?1 
35         IF  (H.LF.l)  Mol 

TF  (M.GE.4R)  M  =  4H 

IF  (N.LF.l)  N=l 

TF  (N.OF.40)  h  =  40 
40      IKOW(M.N) sTROW (M,N) *1 
c 

C   CONVERT  THE  NUMBER  Of  POlNTS/HtCTANGLE  TO  PERCENT 
C    

DO  45  Mb) ,48 
00  45  N=l ,40 
45     THOW(M,N)=( 1 00. *FL0AT ( TROW (M,N) )/Fl OATIKOUMT2) ) *0.5 

c      

C      PRINT  THE  RESULTS 

C    

00  50  M=) ,24 

50      WKITF  (6,60)   ( TROW (M,N)  ,N  =  1 ,20 )  •  ( IROW ( M, N)  ,N  =  2l  ,40 ) 

WKITF  (6,75) 

[iO    55  M  =  25,4P 

55  WRITF  (6,60)      ( T ROW ( M, N) , N= 1 » 20 ) f ( IROW ( M , N) , N=2 1 , 4 0 ) 

WKITF  (6,BS) 

WHITE  (6,100)     PROB.1TYPE 

WRITF  (6,6*)     KOIINT? 

WHITF  (6,70) 

WHITF  (6.9S) 

WHITF  (6,90) 


RETURN 

FORMAT  (1H  ,20X,20I2.1H.,2ol2) 

FORMAT  (1H0,T4,13H  MEASUREMENTS) 

FORMAT  (  lHO.UlHSTEREOfiRAM  SHOWING  DISTRIBUTION  OF  PCIk-TS,  IM  PFRC 
1FNT,  PROJECTFD  ON  LOWER  HEMISPHERE  OF  20-CENT  ImfT FR  WU|  FF  NFT) 

FORMAT  (19<,85HWES1 (3) 

1 EAST  (2)  ) 

FORMAT  (lHl ,57x,5HnORTh> 

FORMAT  (56X,BHsOUTH(l) ) 

FORMAT  (1H0.3BX  ,45H--  — -_.  —  .  — __ ---_...) 

FORMAT  (1H0.9BHOISTR1BIITION  HAS  BEFN  EXAGGFHATFn  By  rOUHLlNG  EaST- 
)WEST  4f»0  NORTH-SOUTH  DTSTA.gCES  ERO"  THF  CFMTER.) 

FORMAT  (IHOtllHSOBAREA  IS  ,bA10.17H  TRTANGtE  TYPE  *  ,Tl) 

FNO 


)  52 
l^o 
156 

15«l 

1*0 
16? 
16* 
)66 
1  6fl 
170 

i  n 

1  7« 
176 
1  7» 
1«0 

182 
]H4 

1H6 
188 

19n 

14? 
1<J4 

\QH 

200 

?n? 
204 

506 

208 

?1  0 
?1? 
21  4 
?16 
21H 
2?0 
??? 

??4 

226 
22fl 
?3n 
?32 
23  4 

?  *6 
?3H 
J4H 

?"»4 
?46 
?4H 
?80 

252 

?S4 
?>-,6 
?58 

2ft  n 
262 

?64 

?66 
?6H 

?7n 
?7? 
?74 
?76 
278 
?80 
?H? 
284 
?Hfc 
?HH 
200 
?R?. 
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h) 


?. 
3. 


4. 
5. 
6, 


7. 

B. 

9. 


1  0. 
1  1. 

]?. 
13, 
1*. 
)S. 

16. 


1  !. 
I". 
1  0. 


LISTING    OF    SUBROUTINE     HISTO 

SUBROUTINE  HISTO  (N,XMa* ,XmIN,X»PR0R»KK, ITYPE) 

C ... ..... .-- 

C  5UBPOUTINE  HISTO 

C 

r  PURPOSE— 

C         Tn  PLOT  HISTOGRAMS  OF  ST&IKFS  ANO/OR  DIPS  EACH  HISTOGRAM  HA* 

C  40  INTFRVALS  AND  THE  HIGHEST  CLASS  IS  SET  EQUAL  TO  100  ARBITRARY 

C  UNITS. 

r 

r 

C  ThE  GOARS  SYSTtM 

C  PURDUE  UNIVERSITY 

C 

C  A,  K.  TURKERt 

C  CIVII     ENGINEERING, 

C  JUNE     )*«>*. 

c 

r  . . .... ._- . — ... 

C  USAGE— 

C  CALL    HISTO (N»XM7X,ZmIN,X»PROH,KK,ITYPE> 

c  WheRe-n  »  nUmber  of   'NiangUla"   Pl.ANES 

C  XMflX   s  MAXIMUM  VALUE  OF  X-ARRAY 

C  X"IN   ■  MINIMUM  V«LUf  OF  x-ARRAY 

C  X      a  ARRAY  OF  OBSERVATIONS  BEING  GBAPHED  (nIPS  OR 

C  STRTRCS) 

C  Pl'OH   ■  UaTa  SET  TITLF 

C  KK      b  1  OR  ?  .   IF  KK  b  1 ,  HISTOGRAM  IS  OF  PIPS. 

C  IF  KK  a  ?,  HISTOGRAM  IS  OF  «TRlKFS, 

C  ITYPE  =  TYPE  OF  TRIANGLES*  BEING  ANALY7FD 

c        „, ... .... . ..... .. ..... 

C  OESCHIPTTON  OF  CONTROL  CAROS- 

C  NONF  RFOUIRFO 

C ... ..... ..... 

C  nfOUIRFO  SUBROUTINES. 

C  NOME 

C  OPTIONAL  SUBROUTINES" 

C  MONF 

r .... ... . . — 

DIMENSION  K(40),  X(N),  I  h  I  <;  T  (?)  «  ProB(B) 

DATA  STAR/JH*/, (IHIST(KK) tKK»l,2)/f,H  DIP   .6HSTRIKE/ 

C    - 

C     CALCULATE  SI7F  OF  INTERVALS,  SET  NUwRER  OF  OBSERVA T I  Of*  S  PER  INTfRVAL 

C    TO  0 

C    ...... . ... ... 

XTERVL=<AMAV*1 .-XMIN) /A0. 

00  5  Ib]  ,40 
S       K( l)sO 

c       - — 

C    PLACE  EACH  OBSERVATION  In  APP^OPRIaTf  INTFPVAL 
c         .- ... 

no  in  i»),n 

Mat  (Xdl-XMlNt/xTERvLI  *1  .0 
10  K(M)bk(M)«1 

C       - 

C    CINQ  HFIflHT  OF  MAXIMUM  INTERVAL  =  KMaX 

C    «FT  (<maX  e  mo  STnRS,ANO  CALCULATE  NUMBER  OF  OhSfRV  AT  I  ONe/sTaR , 

r    ...--..-_-._.-_._._-_......._.__...., 

KMAXbH ( 1 ) 

DO     IS     TbI, 40 
IS  IF     (K(  I)  ,GT.KM(,X)     KMAXoKlIl 

FACKalOO./FLOAKKMAX) 

00    ?0    m=1 ,40 
20  K  (M)  eF»CK«FLOAT  (K  (M)  )  *0.5 

FACKsl ,0/FACK 

c        

C    PRINT  TITLES 

C 

WRITF  (6.10)   IhI^T  (KK)  ,PRO(> 

WRJIF  (6, IS)   ITYPE 

WWI1F   (6,40)   F*CK 


? 

4 

6 

H 

in 

1? 

14 

16 
lfl 

?n 

?? 

?4 

26 

28 

10 

1? 
34 
36 
iH 
40 
42 
44 
46 
4H 
SO 

52 

S4 
S6 
SH 
60 
6? 
64 

66 
6fl 
70 

7? 
74 

76 
7fl 
HO 
R? 
P4 
H6 
8fl 

go 
92 
94 

«6 

1  no 
10? 
104 
106 
109 

no 
n? 

114 
116 

na 

120 

i?? 

124 
126 

1?P 
13(1 
1  1? 
134 
116 
l  3fl 

140 
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c 

c 

?0. 

21. 

22. 

?3. 

24. 

25. 

?5 

?6. 

c 

?7. 

30 

?s. 

35 

?9. 

4(1 

in. 

4S 

31. 

pRINl  OUT  hI?T0G°«f 
XBEGNsXMIN 
DO  25  1=1,40 

XEND=XRFGN*XTFHVL 

NN=K ( 1 ) 

WRITE  (6,45)  XBEGN.xENO, < STAR, Ja 1 ,NN) 
XBERMaXEND 
RETURN 

FORMAT  ( 1 H] , 13HHISTOGRAM  OP  ,Aft,20H  VALUES  FOR  SU8ARFA  ,PA10> 

FORMAT  <JH0,19hTYPE  OF  TRlANGLF  c  ,11) 

FORMAT  UHO.lflHFACH  •  PFPRFSENTS  ,F10.3,14H  OHSFRVaT IO» S . //// ) 

FOHMAT  (lHq,F12.fiilXilH-»lXiF12.6.1X,lnOAl) 

FN" 


1*? 

144' 

146 

l*S 

ISO 

152 

154 

156 

ISO 

160 

16? 

J64 
166 

]  66 
l7n- 


i)  FORSTAT    CROSS-REFERENCE    LISTINGS    FOR    SUBROUTINE    HISTO 


SUWHOUTINF  HISTO 
STATEMENT  NUMrtER    DEFINITION 


REFERENCES 


5 

in 

15 
20 
25 

3n 
35 
40 

4S 


6. 

5. 

9, 

7. 

1?. 

11. 

15. 

1*. 

25, 

21. 

27, 

FORMAT 

17. 

28, 

FORMAT 

18, 

29, 

FORMAT 

19, 

30, 

FORMAT 

2*, 

9 

STATEMENT  NUMBFRS 

DEFINED 

1 

MAXIMUM  NUMrtER  OF 

rfff«ences 

9 

TOTAL    NUMBER  OF 

RFFFKENCES 

1,0 

AVFOAGE  NUMHFR  OF 

RFFFRk'NCES 

SUwRDUTINF 

HISTO 

LAREL 

TYPE 

OEFIN 

ITION 

FACK 

R 

13. 

FLOAT 

R 

13. 

T 

5, 

IHIST 

2, 

0 

ITYPE 

1  , 

J 

?«. 

K 

s, 

0 

KK 

l  . 

KMAX 

if. 

M 

R, 

N 

1  . 

n 

NN 

?3. 

PRn« 

P 

1  . 

n 

STAR 

U 

3. 

X 

R 

1. 

n 

X5TGN 

R 

20, 

XENH 

R 

22, 

XMAX 

n 

1. 

XMTN 

a 

l  . 

X  If  rvi 

R 

». 

REFERENCES 


FUNCTION 
OATA 


DATA 


15. 

15. 
6, 
3." 

IB, 

6. 
3, 

1?. 
9, 
2, 

2*. 

? 
?*. 

??! 

4. 


16. 


7, 

17. 


9. 
17. 
1  3. 

7. 

IT. 


'4. 

25, 


19. 

1°. 
15. 

25, 


12.    ?1. 


1?.    15.    23, 


23. 
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4  4"  n  O  <"t  r-*;  rt*  O  -O  D  O  t»  — I  J"  <Y 
4-4"  J-4"-OOOOvOintA  nr-cOO 
4"     t     t^4nJ1-t44     4    -4-    4"    4-    -4- 


a: 
o 

h- 
> 


< 

cr 
o 
o 
q: 

Q- 

cr 
o 


3   ' 

r 


=) 

i              a 

■f 

CL 

2 

¥                 ■-" 

o 

I    •        1  i- 

LU 

_l 

P  in  m     •  a. 

Q. 

o            -r  ill 

< 

*o        f* 

(/> 

f\J     -H 
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K) 


Typical  output  of   program  VECTOR 


C0MPUI4T10NS  M*  P*0(»MAM  VECTU 


•  PHOO«AH  tO  »f<«Lr5f  TtHH«lN  MUDiMNtSS 


VALUC9  HfcAU  IN 

HO   Of   8U«S"   211   Si!   Slf    CUU5"   F  ■    ASTF.P 
OAT*    FOH«ATaU0F».I  I 


.00    T-AKIS    01KCCIION 


COWUTAtlU"    StSuUlS    FOU     SU8AHC«     ■ttllSfUMt     |.        UBOHLlN    FitLU    Nt"     TOHK     STAtt 


MI»1«»U  (ttV«'IVN  U»T» 
460.  »<40. 


»7t. 

All. 


•  It. 

III. 

!lll 

til. 

tTf. 

III. 

•  rj. 

ATI. 

AM. 

♦?»i_ 

»•«. 

♦♦♦. 

III. 

•  ♦I. 

tM7 

AJJ. 

ItT. 

AM. 

•  It. 
All. 

»»3. 

_»t#l_ 

*fl. 

•TATIST1CAL    HEASUMCS    OF    FLtVATlON    UATA 

"CAN"       At>.62<9  <.<llr.ni     lill.1000  STANIiANU    111  V  I  At  IONa  J6.170I  HtMN 


THl    MouHAM    lb    -.a*    FITTING    T»PB     ?    !HIANt»LE«. 


AM. 

til, 

*T|. 


til. 
til. 

«ll. 
tlf. 

tl». 

*•»• 

stl.oeii      Hi  m  ww.    iit.ntt 


ITAUSTICAL    HCASIJMfcS    OF    STkiai^    of    Tw|a*oi.f.s 

h(ANa  A.ftOO*  .*<lANre»    Mll.tVn?  ^TANUAMO    DtVlAIION* 

•  lATllllCAl.     «IIA>0*FS    OF     M||>>    .|F     f  M  I  ANULf.  !. 

•AtANA  A.A1a*£  WAMIAMrtA  V.Kul-l  SfANOAXU    UtVlATlriNA 


•O.01BH  HA.|NU«.  90.0101  >l>|«Ml        •tt.ttt* 

J. lilt  HAAIKUi.  IT. tilt  MINIMMa  ••. 
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„.w>„  B,    wtCTO«  STRENGTH  ADC  R|"     718.7,78.     «.  O.ttSAli.  «Ht«C  MRl/NUMIE"  Of  t»I»WCH.        

^Ho«^W."oN%£rs"?Da"hT^5s-01S>C»»10H  f ACTOR. II     »l 1.0.  V.LUCI  INOICATI  VARUKt  TfM.ll. 

THf  VARIABILITY  FACTOR.  V 


ELEVATION  RANGE/LOO  USTKI  .  II  A  Tli0**  

VECTOR  «R-  ALPHA.   -O.0O0S2   iCTA.    0.00011   OAKAIA.    1.00000 
0.JJ11 


THf  DIRECTION  COSINES  Of  THE  HE 

TMC  ORIENTATION  Of  THE  REAN  PLANE  Is  A-^OLtoi»'-»-~~ITHlKt.-ll.».at.  Dl* 


H!tfo<).AN  of DIP    VALUfS  'OR  SUBARCA  VEEOSPORT  1.   ORUHLIN  FIELD  NC»  Y0««  STATE 

TYPE  Of  TRIANSLE  •  I 

CACH     •    HERRE5ENT5  0.570    Q85(lMWA  '  ION*. 


1.4*1 ' . < 

4.S5214I 
S. 007355 
4.462569 
S. 917783 

».17<vvr 

T.28J426 

7.711640 

a.  l v3oin . 

1,649060 

9,l»<<«2 
9.56««98 
10.01*710 
10.66992* 

U.J.OJSI 
11.01*5*4 

11.7*5994 
1). 201200 

14.111636 
14.566850 
15,. 022066 


1S.A77270 
IS. 932491   i 

lo.-anoa  i 

16.1*2920 
I7.2t01)« 
II." 


S..6256* 

a-AUioi 

6.372997 
6. 828212 

l.illlil 
7.738640 
8.193454 
_Ji-l4»06.a 
9,104282 
9.SS9496 
10.014710 
10.469924 
10.929130 
11.3803S2 
11.835466 
12.290700 
12.7*S9Y» 
13.201200 
13.656*22 
11630 
l*.566at0 
10,022104 
-  IS.6TT2UL 
IS. 932492 
10.387706 
1A.R*Z9?0 
17.298134 
17.753348 


NOR7H 
-o-o-O-o-o-O-o-o-o-O-O- 0-0-0    00000    0.    00 

-O-O-O-O-O-O-O-O-O-O-O-O  I  1  I  I  I  I  I  I.  I  1 

-O-O-O-o-O-O-O-O-O-O  000000000  0.  00 
-0-0-0-0-0-0-0-0-0  0000000000  0.  00 

-a-a-B-O-0-0-0-0  lllllllllllli  00 

-0-0-0-0*0-0-0  000000000000  0.00 
-0-0-0-0-0-0  0000000000000  0,  00 

-n-Q-fl-ii-fl  Qouoauooooooooo.oo 

-0-0-0-0  000000000000000  0,  00 
-o-o-o-u  O0OOOOO00O0OO00O,  00 

-o-o-o  muh  n  i  n  t  H  h  o,  n 

.0-0    000000000000000000.00 

•  0-0    09000000000000000    0.    00 
afl-0   flaeOUOOO0OO0QOaQ0B.ee 

-OoOOOOOOOOOOOOOOOOOO.OO 
-0    0000000000000000000.00 

>a  uoooflD-OflflQ-fl-O  o  n  -o  1000.00 

-OuOOOOOOOOOOOOOOOOOO.  00 
-OOOJOOOUOOOOOOOOOOO^, 

ll  1  I  I  I  1  I  i  I  1  I  1  i  I  I  1  I 

-0  000000  OOOOOOOOOO 

•0  0000000000000000 

-OUQ.OQQuOUO.OQOOO 

-oeooooooooooooomo 

•  EST 

-ji  00000000.00000000 

-ooouooooooooooo 

•  0  0000000000000  0  Q)  0  (JJ  0  0.  0  0 

-ooooooauouooaoo 
•ooooooooouoooooooooo.  00 

-ODOOOOOUOOOOOOOOOOOO,  00 
•0  0  0  OQQQUQ0Q0QOQQ0000.99 

•  ooooooooooooooooooon.  go 

-OOOOOOOOOOOOOOOOOOOO.OO 

-.0-0  0  UQQQOflOOOfi  a_Q  JL_aajLa»_J  a. 

-0*0  OO0000OOO00O00QO0  0.  00 
-0-0  0O0O0O00O00000000  0.  00 

-a.u-a  QaQaaaooaaoouoaa  -tu  a  a 

•0-0-0-0  oooooooooooooooo.  00 

-0-0-0-0  000O0OOO00000OOO.  00 
-J. 0.0.0.0  0-  0  J-O.H-0-  I  I  I  I  I  I  I  I  G.  I  I 

-u-u-c-  j.0-0  uooooooooooooo.oo 
-0-0-0-0-0-0-0  u  00000000000  0.00 

-J1-J--J-.U-J3.0.0.J  o  o  a   ii- a  u  a  o  a.o  o  a.-u  a 

-0-0-0-0-O-O-o-u-O  0000000000  0.  00 

•o-o-o-o-o-o-o-o-o*o  aooooooooo.  oo 

-a-juo-o.o-o.o-o.u.u-o.o  o  a  aji.ao-a  a.  a  a 

•o>o-o-o-o-o-o-o-o-o-o-e-o-o  oooooo,  oo 

•o-o-o-o-o-o-o-o-o-o-o-o-o-o-o-o-o  ooo.oo. 

__  _  _         south m 

MWU    It    RCEOSPOHI    1,      ORUNLIN  P1ELO    -f    YORK    S7AIE 


0  0  o 

0  »  3 

0  0  0 

C  5  0 

V  0  9 

0  0  0 

0  0  0 

0  0  <i 

0  0  0 

0  0  ') 

0  0  0 

0  0  0 

0  0  0 

1  0  0 
0  0  c 
0  0  0 
C  0  0 
0  0  0 
0  0  0 
0  J  0 
0  0  0 
0  0  0 

0  0 

0  0 


o-o-o-o-o-o-o- 
a  a-o-a-o-o-o- 

BOO  0-0-0-0- 
0  0  0  0  0-0-0- 

a  o  o  o  o  o-a- 
ooooooo. 

0  0  0  0  0  0  0 
lllllll 

o  o  a  a  o  o  o 
o  o  o  o  o  a  o 
iingii 
o  a  a  o  o  o  a 

0  0  0  0  0  0  0 

a  a  (i  o  o  o  a 

0  0  0  0  0  0  0 
0  0  0  0  0  0  0 

a  a  o  o  a  a  a 

0  0  0  0  0  0  0 
0  0  0  0  0  0  0 

a  a  o  o  a  a  o 

o  o  o  o  a  o  o 

a  o  o  o  o  o  a 

a  a  o  o  a  a  a 

a  o  o  o  o  a  o 


0-0-0 
0-0-0 

0-0-0 
0-0-0 

«-»-« 

0-0-0 

0-0-0 
0  1-1 
0  0  0 
0  0  0 
«  0  t 
0  0  0 
0  0  0 
9  0  3 
0  0  0 
0  0  0 

til 
too 

e  «  5 
t  i  0 
0  0  0 

0  0  9 
C  0  0 

0  0  0 


llllllill 

000000000 

ooooooooto 
ooooooooeo 
oooooootoo 
ooooeooooo 
1111  111!  II 
oooooooaoo 
oooooooooo 

lltll  1L1II  « 

ooootttoti 

oooooooooo 
ooooo  a  o  a  a  a 
oooaotoooo 
oooottooto 

II  111  IIHI 
oooooteooo 

OOOOOOOOOO' 
g  a  a  p  a  a  a  a  fl-o. 
0  0  0  0  0  0  0  0-0-0- 

oooooo  o-o-o-o- 
4  a  a  a  a-a.a-a-o-a. 
0  0  o-o-o-o-o-o-e-o- 

.0-0-0-9-9-0-0-0-0-0- 


Ill 

0    0    0 

0   t  t 

0  0  0 
0  0  0 
0     0     9 

ant 

0  0  0 
0    0    0 

0   I   t 

9  t  0 
0  0  0 
0   10 

0  0  0 
a-»el 


0-0-0-0-0 

.t-o-B-v-« 
,0-0-0-0-9 
•0-0-0-0-0 
•l-l-l-l-t 
■o-o-o-e-o 

■0-9-0-0-0 

.g.g.g.|-g 

.0-9-0-9-0 

■o-o-o-o-o 

g.g.g.g.g 

0  o-o-i-o 

0    0-0-0-0 

>  l-l-l^l 

0    0    0-9-0 

0  0    0-0-0 

1  I   l-4?l   

0    0    0    0-0 

0    0    0    0-0 

J-IJLi-J 

HIM 

0  0    10-0 

III  l-a 

1  I  I  l-l 

lllllll 

a  a  jl  a-9 

t  0  0  0-0 
0   0   0   0-0 

10  0  0-0 
0    0    0    0-0 

a  o  t  t-t 

A  J  l-l-t 

0  0  0-0-0 
0   0    0-0-0 

XJiWsl 

0  l-t-l-l 
0  0-0-0-0 
l.o.n.fl.a 
0-0-0-0-0 
9-0-0-0-4 
t.o.a.o.o 


0-0-0-O-O-O-l-t 

0-0-0-0-0-0-0-0 

o.p-0-.t-o-.-p-t 


TR1AN-RC    TTPf    •   I 


711    ■tASORE»CNta 
HIM9»jM_JHO«INa_BIJTR|sllT|ON  Of   POINTS.    IN  PtRCENTi    PROJECTEO   ON  LORER   HtRltHatm   of   10-tEwT lutTtR   -lAIT  ItM 
01ITRIIUTION    HAS    BEEN   EAltOCRATCO    ST    00URLIN9   EAST-REtT    AW)    NORTH-KHITH    DIITAXCI1   MOD   TM   CtNTlR. 
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I)  LISTING    OF  PROGRAM    LOADS 


HKijfiWflM    LOAnsUMfU'l  .0U1HUT.IAPF2»TARE5sInPU1.TAPE6=OUTRUT) 
PROGRAM      |  OADS 


MURHOSF-- 

TO  LUAfl  X.Y.Z  CUUHUTNATf  DATA  FROM  CAR05  TO  TAPE. 
flATA  TAPt  IS  TAPF?   , 


THE  IjCAHS  bYSTtM 
PURI  UE  UNIVERSITY 


A,  KF1TH  T 
CIVTl  ENGI 
JUNF,  is^y 


UMNFRi 
NEEPlNG, 


general  Information 

THE  PROGRo"  CAN  hANhLE  T«0  TYPES  OF  Z-PATa.  FLOaTI 
INTEGER.  FLOATING-POINT  /-VALUES  (FOR  EXAMPLE  TOPOGRAPH 
ELEVATIOnS)  ARE  REaU  AND  PLACtO  ON  TAPE  ►ITHPUT  MOOlFlC 
TNTEGFR  /-VALUES  (HOH  LXaMPLE  SOTL  TYPF  rODFSI  APF  CONV 
CORRESPONDING  F Loa T I NG» POINT  COST  OR  HFNFFIT  VaLUF  SUPE- 
rONlHOl  RARAMEIE^.  THE  X  AND  Y  COORDINATES  AND  THFSF  V 
PLACED  O'J  IHF  lA>JF  ,  AND  THE  CODE  NUMBERS  ARE  DISCARDED 

TAPF2  CONTAINS  a  TITLE  Card  FOLLOWED  BY  a  SFPIES  0 
FaCH  CONTAINING  1000  -  OtSERVATIONS  (300(1  WORDS)  IN  PEL 
HlMANY  FORMAT.  THE  I  AST  hECOHQ  WILL  BE  OF  VaHIARIF  LEN 
DppENDltgG  ON  THE  nUKhEK  of  OHSERVAT  IONS  t 

The  program  CMN  COMBINE  SFVEHaL  sfts  of  siihaREa 
mEaSURFO  In  ITS  nwh  coordinate  system.  InTO  a  single,  la 
HAVING  a  Uniform  cOUHDlN/xIt  SYSTFM.  ALL  DaTa  MUST  HOWE 
kEFErENCFU  TO  AXFs  OF  CONSISTENT  DIRECTIONS.  «ITM  X  pile 
RIGHI  AND  Y  POSITIVE  Upw,\HU.  SUhAREaS  MAY  HAVF  OTFKERE 
DATA  SUhsETS  ARE  SFPaRaTFU  bY  ENO-OF-RFCORD  CARDS, 


NG-P01NT  AND 

IC 

AT  ION, 

FRIED  TO  THF 

LIED  AS  A 

ALUFS  ARE 


F  RFCORDS, 

OCATAHLF 

GTh 


DATA.  EACH 
Pr-F  DATA  SET 
Vh,R  BE 

ITIVE  TO  THE 
NI  SCALES. 


TAPt  REQUIREMENTS  - 

TaPF2  -  ThE  OUTPUl  DATA  lAPF  ON  WHICH  THE  CaPD  DATA  IS  IOADFD. 

INSCRIPTION  OF  CO'  THliL  C/iRUS  - 
(1)  »MYOAln  r«i<u- 

•values  (cols  ,-7)  -  imuicotfs  thf  z-values  a 

POINT. 

•codes     (cols  i  — o >  -  indicates  the  /-values  a 


RF  FLOAT  INT, 
RF  INTFGFP. 


(?)  title  card  - 

eactor-lculs  1-30)  alphanumfr1c  factor  identification. 
apfa  -(('uls  <n-?o>  alphanumeric  area  toe  n  f  1  e  icat  ion, 
itypf  -(col  78)     if  data  measures  costs   .  ttypf  =  ? 
if  data  measurfs  hfmefiis.  itypf  =  3 

the  following  cards  should  pe  hfpeated  fur  fach  dak,  suhsft  - 

(3)  surtttlf  card  - 

suhtiticols  1-70)  subarea  identification  (alphanumeric). 

U)  origin/scale  Card  - 

XORG  -(COLS   1-10)   X-COORDINATE  PF  SURAPEa  ORIGIN  IN 

ovfr-ai  l.  coordinatf  system'. 
yoro  -icois  11-20)  y-coor01nate  of  surarea  origin  in 

ovfr-all  coorotnatf  system, 
scalemcols  21-30)  relativf  sralf  of  siirarea  (if  blank, 

scale  is  sft  eoual  to  one), 

the  following  cards  are  used  with  fi oaTing-point  z-valufs  - 

(bA)  FORMAT  CAHD- 

NPT   -((OLS   1-2  )   NUMHER  OF  OHSE.RVATTON  SETS/DATA  CARD, 
FORM  "(COLS   /-lb)        FORMAT  OF  DATA  CARDS. 

i*A)  data  cards. 

CONTAIN  DPI  SETS  OF  X.Yt  AND  Z-COORD 1 Na TF S . 


10 

12 

1* 
lis 
IS 
?o 
22 
?ft 
?6 
28 
in 
32 
3ft 
3<S 
3* 
40 
42 
4ft 
4fe 
4H 
SO 
S? 
S4 

s* 

SH 

60 
62 
6* 

6b 
68 
70 
7? 
74 
76 
78 
80 
82 
RA 
HIT, 
HH 
111 
■>» 
■>4 
-)(S 
9A 
I  Ol) 
1  1? 
1  04 

10* 
10" 

110 

)  1? 

1  1ft 

1  16 

i  ia 

1  20 
122 
124 
l?6 

126 

llo 
IV 
1  1ft 
1  ^6 
138 

140 

1*2 
1  Aft. 

I  on 
IftH 
ISO 
IS? 
ISA 
ISA 
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10, 

1 1 , 


12. 
1  3. 

16. 


(74)  ENP-0F-kf  C(,HlJ  CAHU  - 

THIS  SIGNIFIES  THE  END  OF  ONE  SUbSFT  OF  OATA. 

THE  FOLLOWING  CaROS  ahe  used  WITH  INTEGER  Z-VALUES  - 

(58)  FORMAT  CARU- 

NPT   -(COLS   1-2  )   NUMBER  OF  OBSERVATION  SETS/DATa  CAPO. 
NFAC  -(COLS   3-*  «   NUMBER  OF  ASSIGNED  z-codes. 
FORM  -(COLS   7-76)   FORMAT  OF  OATA  CARPS. 

(6B)  VALUE  ASSIGNMENT  CaRO(S)  - 

THFSf  CAHDi  CONTAIN  A  cEQUFNCE  OF  VALUFS  TO  HE  ASSIGNED 
TO  THE  INTEGER  Z-VALUES  (COOES).  EACH  Z  WILL  BE  ASSIGNED 
THF  *FaC*  VALUF  AT  I  HE  POSITION  IN  IhF  SEQUENCE  CORRFSPONO 
-ING  TO  ITS  COPfcU  Z-VAIUE.  (I.E.  IF  Z=S.  Z  BECOMES  FaCIs)) 

THF  VALUES  '.HE  READ  FROM  5. COLUMN  FIELDS  ( 1-5.6-10. 11-1 P. 

lf.-2d.ETC) .THERE  bEING  A  MAXIMUM  OF  16  VALUES  PER  CARD. 
THE  TOTAL  NUMBER  OF  VALUES  SHOULD  NOT  FXCEEO  FIFTY, 

(78)  DATA  CARDS- 
CONTAIN  hPT  SETS  OF  X,Y»  AMP  Z-COORDINATES. 

(88)  ENP-OF-RECORO  CARD  - 

THIS  SIGNIFIES  THE  tNO  OF  ONE  SURSEl  OF  OATA. 

THE  FOLLOWING  CaRD  MUST  r>E  ADDED  AFTER  Thf  LAST  OaTa  sursfT  - 

(9)   END-OF-RFCUHD  CAHQ  - 

THIS  SIGNIFIES  IHE  END  OF  ALL  OATA. 


SPECIAL  nUIE  - 

whFm  MULTIPLE  OBSERVATION  SFTS  PER  CARD  ARE  USED  AND  THE  LAST 
C'RU  IS  NUT  COMPLETELY  FILLED  ONF  ADDITIONAL  DUMMY  OBSERVATION  MUST 
RE  ADDED  TO  INDICATE  This  SITUATION. 
THIS  OHSFHVATION  SHOULD  Mt  - 

(A)  X  =  o.U»  Y  =  0.0.  7    =  -999,0  FOK  FLOA T ING-PO I Nl  7-PATA 

(R)  A  BLANK  RECUHD  FOR  INTEGER  7-0aTa 


PEOUIRFO  SUBROUTINES 
NONE. 


OPTIONAL  SUBROUTINES 
NONF  , 


DlMENSIUN  NVAL(ln)'  EAClOR(J).  AREA(3).  SUBHT(7)»  FCRm(7>,  X(J000 
1).  Y(inUO).  Z(innO).  XT(IO).  YTIIO).  ZTIIDIi  FACI50).  IPATa(?) 
JNTFOFH  FORM 

INTEGER  SYSOT.SYSIN.SDE2 
DATA  SYS0T,SYSIN.SDF2/f.»b»?/ 
DATA  (I  DATA  (  I)  .1  =  1  »2)/6H«VALUE ,bH»CODEs/ 

REMIND  SDF2 
KUUE=3 

KOUNT=l 
KT  =  (i 
KNUM=n 


-FAD  OMYDATA  CARP 


REAP  (STSIN.110)  MYDATA 
IF  (mYoaTA.EQ.IPATA ( 1) )  KOOt=l 
If  (MYDATA. EQ.IDATA  (2))  K0ut=2 
IE  (K0PE.F0.3)  GO  TO  jn0 

<FAD  TITLE  CARP 

REAP  (SYSIN.HS)  FACTOR. AREA, ITYPE 
lUTPUT  TULF  CARP  A<;  FIRST  RECORD  ON  TAPE? 

WHITE  (SOF?)  FaCTOx.AHEA. ITYRF 


A 

158 

A 

lf>U 

A 

16'2 

A 

164 

A 

1  66 

A 

168 

A 

1711 

A 

172 

a 

1  7* 

A 

176 

A 

178 

A 

1  wi) 

A 

1"? 

A 

1  84 

A 

1  86 

A 

1  B8 

A 

190 

A 

192 

A 

194 

A 

196 

A 

198 

A 

2  0(1 

A 

202 

A 

204 

A 

206 

A 

?0B 

A 

210 

A 

212 

A 

214 

A 

216 

A 

218 

A 

220 

A 

222 

A 

22* 

A 

226 

A 

22H 

A 

2  3  0 

A 

212 

A 

234 

A 

236 

A 

238 

A 

240 

A 

242 

A 

244 

A 

246 

A 

248 

A 

?50 

A 

2S2 

A 

2^4 

A 

2Sh 

A 

258 

A 

260 

A 

262 

A 

?64 

A 

266 

A 

2  68 

A 

2  70 

A 

272 

A 

274 

A 

276 

A 

276 

A 

2*n 

A 

242 

A 

2«4 

A 

286 

A 

288 

A 

?90 

A 

?92 

A 

2«4 

A 

296 

A 

,04 

A 

1(111 

A 

102 

A 

104 

A 

3  0  6 

A 

108 

A 

110 

A 

112 

A 

314 

A-44 


18. 
19, 


20. 
21, 


22. 
23, 
?4. 
25. 


Ph. 


?7, 
?8, 
29, 
30, 


31  • 

3  *. 

'*  • 

35. 

lb. 
37. 

IB. 


4  I  . 

42, 

43. 


Ml. 
'  I. 


e 

C    SET  UP  PnbF  HFADIMG  FOR  SlrtRT  OF  OUTPUT 
C    

WHITE  (SYSOT.l?0>  FACToR.ARfcA, TTYPF 

WHITE  <SYSOTil25> 

C    

C    WORK  WITH  EACH  SURAHFA  In  TURN 

C    

C    RFAO  SUBTITLE  CAPO 

C    

5  RtAO    (SYsir*tl3n)    SUPTIT 

IF     (EOF.SYSIM     7(/il0 

C         

C         READ    ORlr.lN/SCALF    CAHD 

c         

1(1  READ     (SYSINH35)     yORG»YOHG,SCA|  E 

IF     (SCALE, EG. 0,0)     SCALK=1.(I 

WRlTF  (SYSOTtUO)  bUHT  IT  i  XOhG.  YOrG,  SC  ALE 

GO  Tn  (  15t40)  »  KOOE 

c 

C    uFaD  DATA  wlTH  Z  At,  FLOATING  POlNl  VALUES 

C    

C    RFAD  FORMAT  fARD 

C    

1?     RtAO  (SYSIN»145)  NHT.FORM 

C    

C    RFAO  OAT*  CARDS 

c       

20     Rt«l)  (SYSIN.FORM)   < X t  ( I )  , YT <  I  )  ,  ZT ( I )  »  I  e 1 , NP  I  ) 

IF  (FOF»SYSIN>  35t25 
25     DO  3d  N=l ,N'PT 

IF  (ZT  (N)  .EQ.-999.)  GO  10  3fi 

C    

C  FIX    UP    THE    COORDIMftTtS 

c . 

X(KOUNT)=XT<N)»SCALF*XOWG 

Y(KOUMT)=YT(M)*bCALF*Y0RG 

7<K0UMT>  =7T <M) 

HOOMl  »KOUf»T*  1 

IF     <KOIINT.GT.1000>     GO    TO    65 
30  CONTINUE 

GO    To    ?0 
C 

35             GO    TO    S 
C 
c  ... — _ 

C    READ  I    AS  INTEGER  CM'ES  AND  CUNVFHT  TO  VALUF5 

C    

C    RFaO  FOR^Al  CARD 

C    

40     RtAO  (SYSIN»15()I  NHT»NF«CiFORM 

C    

C    RFaU  COOFh-VALUES  THaNSFOHMAT IONS 

c    

READ  (SYSIN»155)  (f AC  (  I  )  . I  =  J  ,NFAC) 

C    

C    uFAU  DATA  CArOS 

C    

45     RtAO  (SYSIN.FORC)   ( XT ( I ) • YT ( I ) , NVAL ( I ) » I»l »NPT) 

IF  (FOFtSYSlN)  60, 50 
50     00  55  w=l «MPT 

IF  (MVAL  (M  .FCj.O  I  GO  [0  S5 

C    

C    FTX  UP  THE  COORDTNATcS 

c         

N/sNVAL (N) 

X(K()UNT)=XT(M)»SCALF*XORG 
Y(mtUM)=YT(M)«SCALF  +  YORG 
I (KnUNT) =FAC (r?) 
M)JnI=KOUNT»  i 

IF  (KOUNT.GT.1000)  GO  To  65 
55     COnTTNUE 
r-,0    fn  45 


60 


GO  10  5 

l>  NUI<l3<MUH  ♦  1 


A 

316 

A 

31« 

A 

3.20 

A 

•322 

A 

324 

A 

326 

A 

32H 

A 

3  3(1 

A 

332 

A 

334 

A 

336 

A 

3?H 

A 

341) 

A 

342 

A 

344 

A 

946 

A 

34H 

A 

350 

A 

3^2 

A 

154 

A 

356 

A 

3^» 

A 

36-) 

A 

162 

A 

364 

A 

366 

A 

36-1 

A 

3  71) 

A 

372 

A 

374 

A 

376 

A 

37H 

A 

980 

A 

3*2 

A 

1M4 

A 

386 

A 

^H« 

A 

390 

o 

3<»2 

A 

394 

A 

396 

A 

398 

A 

400 

A 

402 

A 

-04 

A 

406 

A 

4fH 

A 

410 

A 

412 

A 

414 

A 

416 

A 

41  1 

A 

420 

A 

422 

A 

424 

A 

426 

A 

42M 

A 

430 

A 

432 

A 

4  14 

A 

4  16 

A 

43* 

A 

4  4l) 

A 

442 

A 

444 

A 

440 

A 

44rt 

A 

4^" 

A 

4^2 

A 

454 

A 

.Sft 

A 

458 

A 

460 

A 

462 

A 

4*4 

A 

466 

A 

468 

A-45 


55, 
56. 

=  7. 
58. 


59. 
60, 
61. 


62, 

83. 
8*. 
65. 
66. 
87. 
68, 
6^. 

TO. 
71. 

7?. 
T3. 

74. 
76. 
78. 
77. 

re. 

79. 
BO, 

fil. 
R?. 
R3, 
84. 

85, 
R6, 

B7, 
B8, 
B9, 


90, 


"J2, 
93, 
94. 
95, 
96, 
97« 


98, 

4<v. 


lno. 


c   

C    WHITE  BLOCK  OF  1000  UbSERVAlIUNS  ON  TAPE  2 
C    

KOUNT=l 

WHITE  (SDF2)  (X  (I)  .Yd)  ,Z(  I)  .1  =  1.1000) 

GO  TO  (20.45) »  KODt 
70     KOUNT=KOUNT-l 

c    

C    WHITE  FINAL  RECORD  On  TAPE2  aid  THEN  AN  FND-OF-FILE 
C    

WHITE  (SDF2)   (X(I) ,Y (I) tH I) ,l=l,KOUNT) 

END  FILE  SOF2 

rEwInd  SDF2 

c    

C    CHECK  DATA  ON  TAPE  n»  LISTING  IT 

C    

READ     (SOFp)     FACTOR. APEAtI TYPE 

WHITE     (SYSOT.120)     FACTOR. AHtA.HYPE 

NTOT=KNUH»1000+KOUNl 

WHITE    (SYSOT1I6O)    MOT 

IF     <  KNIJvt.LT.  1)     GO     TO    90 

00    85    KKsi.KNUM 

REAO     ISHF2)     (X  (I)  »Y(I)  »Z(I)  ,1  =  1.1000) 
IKT=250 
75  I8Fr;=] 

DO    hO    IlaltlKT 
IEND=IBEG*3 

»i«TTE     (SY50TU65)      (X(I)  lYIItTfl  I)  »  I  =  THEffrTE  NO) 
80  IRFO=18EG+4 

R5  CONTINUE 

IF     (KT.tO. 1 )     GO    TO    ^5 
90  KT=l 

READ     (SUF2)      (X  I  I)  .  Y d)  ,Z  (I)  1 1  =  ] .KOUnT) 

IKTeKOUNT/4 

GO    TD     7b 


c 

ion 

105 

C 

110 

lis 

120 


125 


130 
135 

141 
US 
150 
155 

160 


165 

170 


ILFT=K0UNT-IKT«4 

IF     (ILi-T.NF.O)     WRITE     (SYS0T.165)      (  X  (I  )  ,  Y  (  I )  •  Z  (  T  )  •  I»  IREG.K  OI1NT  ) 

RE*Ini)    SOF? 

GO    TO    105 

WHITE  (SYSOT.170)  MYDATA 
SHlw 


FORMA 

FORMA 

FORMA 

1ME  GC 

2»29H- 

3ITY  - 

4ATA  F 

58HIIY 

6YHE  = 

FOHma 

IE  FOL 

2TESSC 

FO*ma 

FORMA 

FORMA 

FUHMfl 

FOkMA 

FORMA 

FORM* 

1  IS  . 

2X.25M 

3Z-VAL 

F0KM6 

F0RM4 

IT  CHF 

FN 


T  (Afe) 
T  (3A10 
T  I1H1. 
a^S  SYS 


.10X.3A 
55<»29h 
TE« 


LOAD  0 
OR  --,/ 

PF   -.1 

3  WHFN 

T  UHOi 

LOSING 

A|  E,/, 

T  (7A10 

OF10 

(1H0. 

d?.4 

(21?, 

H6F5 

(  JH  » 


I 

(  8£AD 
IF   )  ,/ 
T  (1H  , 

r  oho, 

CK.  »MY 


ATA  KRO 
/•40X.9 
2.36H  I 
DATA  m 
1  OX  .f,9H 
jyRARtA 
6X,1HX, 
) 
.0) 

10X.751 

X,7A)0> 

2X.7A10 

.0) 

55mT(iTa 

,39H    LI 

ACROSS 

/) 

12F10.3 

DATAx     . 


10.6A, 12) 

,/,56X,29HI  T 

I,/,56X,29Hl  PUROUF  UNIVFP9TTY  VERSloN  I»/,56X 

—  _  —  —  _.«. »///3?X,56r(DATA  PREPARATION  ACTIV 

m  CARDS  TO  TAPt,//V,40X,30HTHlS  TAPE  CONTAINS  D 
HFACTOR  -  ,3A10,//,40X»9MARFA  -  ,  1A \ 0 , // , 4  OX , 
TYPE  =  2  WHFN  DATA  MEASURES  COSTS .,/, SOX , 39M  IT 
tASUKES  8ENFF ITs. ,/////) 

'he  tape   is  being  loaded  with  card  data  from  th 

S--,//»36xl?hSUR*"EA  NAME,3^x,2,H0RIGINCOOR[)TNA 
IoX.IHY,/) 


0i5Xt3Fl0,3) 

) 

L  NUMBER  OF  POImTS  FOR  THIS  FACTOR  IN  THIS  AREA 
STING  OF  DATA  PL«CEn  ON  TaPf  FOLLOWS »//// 1 h  ,«6 
|HE  LINES  ).////. 1H  ,4(30H  x-COORC,   Y-COOR(i. 

) 

kHOR««»  CONTROL  CARD  ERROR , umyDA T A  CARD  DOFs  NO 
A6) 


A 

470 

A 

472 

A 

47<» 

A 

4?i 

A 

47* 

A 

480 

4 

46? 

A 

4«4 

A 

4H6 

A 

4*H 

A 

49  0 

A 

49? 

A 

494 

A 

496 

A 

49H 

A 

son 

A 

502 

A 

504 

A 

506 

A 

508 

A 

510 

A 

512 

A 

51  4 

A 

516 

A 

sis 

A 

520 

A 

522 

A 

524 

A 

5?" 

A 

5?6 

A 

S30 

A 

532 

A 

534 

A 

536 

A 

53H 

A 

54  0 

A 

54? 

A 

544 

A 

546 

A 

54* 

A 

550 

A 

55? 

A 

55h 

A 

556 

A 

55H 

A 

560 

A 

56? 

A 

564 

A 

566 

A 

66* 

A 

570 

A 

57? 

A 

574 

A 

576 

A 

57H 

A 

580 

A 

si? 

A 

SH4 

A 

6*6 

A 

5** 

A 

S90 

A 

64? 

A 

6«4 

A 

64* 

A 

6<JH 

A 

60  0 

A 

60? 

A 

6  04 

A 

606 

fi 

60* 

c/> 
a 
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IMPORTANT  NOTICE 

Program  TREND  involves  matrix  inversion  procedures.  The  CDC 
65OO  has  a  60-bit  word.  Due  to  this  long  word-length  single  precision 
gives  sufficiently  accurate  results. 

However  on  other  computers  with  shorter  word  lengths  (including 
all  IBM  360' s  and  709^'s  and  similar  machines)  double  precision  must 
be  used. 

The  following  variables  need  to  be  declared  double  precision :- 

PAR,  XX,  A,  B,  BB,  XINV,  SSZ,  SSQZ,  SSR. 

BBB  in  Subroutine  RITER 
AM  in  Subroutine  INVER 
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o)    LISTING    OF    PROGRAM   TREND 


PROGRAM  TREND ( I NPUT, OUTPUT , PUNCH i TAPE  1 »TAPE2, TAPE4, TAPE5* INPUT, TAP 
1E6-0UTPUT)  : 

PROGRAM  TREND 


PURPOSE  -- 


TO  PERFORM  TREND  SURFACE  ANALYSIS  USING  POwER  SERIES  EXPANSION. 
POLYNOMIALS  UP  TO  THE  FIFTH  DEGREE  CAN  BE  REQUESTED.  THE  PROGRAM  HAS 

EVOLVED  FROM  AN  EARLIER  SERIES  OF  PROGRAMS  CALLED  *   PULYFIT  *  . 

THE  PROGRAM  IS  USED  WITHIN  THE  GCARS  SYSTEM  TO  CHECK  DATA  FOR 
GENERAL  VALIDITY  .  AND  TO  CONVERT  TOPOGRAPHIC  MEASURES  TO  EARTHWORK 
COST  VALUES. 


THE  GCARS  SYSTEM 
PURDUE  UNIVERSITY 


A.  KEITH  TURNER. 
CIVIL  ENGINEERING, 
JUNE  196V. 


GENERAL  INFORMATION  _  _    

THE  PROGRAM  CAN  FIT  A  SURFACE  TO  AN  ESSENTIALLY  UNLIMITED  NUM- 
BER OF  OBSERVATIONS.  APPROXIMATELY  lliOOO  OBSERVATIONS  HAVE  BEEN 

SUCCESSFULLY  ANALYZED  AT  ONCE.  HOWEVER,  CaRE  SHOULD  BE  TAKEN  TO  SEE 

that  values  representing  sums  of  powers  do  not  exceed  thf  maximum 
size  that  can  be  handled  by  the  computer.       

observations  must  be  defined  in  terms  of  two  independent  vari- 
ables (  x  and  y  coordinates  )  and  one  dependent  variable  t.z  )..  the 
program  computes  a  series  of  least  squares  regression  surfaces 
(  trend  surfaces  )  to  these  data.  program  options  include  - 


(1)  CALCULATION 

(2)  LISTING  OF 

(3)  LIST  OF  COE 

(4)  ESTIMATE  OF 

(5)  ANALYSIS  OF 

(6)  LISTING  OF 

ALL  OBSE 

(7)  GRIDDED  VAL 
ANDI8)  CONTOUR  MAP 

PRODUCED 
THE  VARIOUS  OPTION 
CONTROL  CARDS  ARE 
THUS  DESCRIBED  IN 


OF  REQUESTED  TREND  SURFACES 
MATRICES  USED  IN  DETERMINING  THE  COEFFICIENTS 
FFICIENTS  

ERROR  IN  EACH  COEFFICIENT 

VARIANCE  STATISTICS 
ORIGINAL  DATA,  COMPUTED  VALUES,  AND  RESIOUALS -FOR- - 
RVATION  POINTS 

UES  OF  SELECTED  SURFACES  On  A  DESIGNATED  SPACING 
S  OF  THE  TREND  SURFACES  ON  THE  PRINTER  {  THESE  ARE 

BY  SUBROUTINE  PRTC0N4  ), 
S  ARE  REQUESTED  BY  APPROPRIATE  CONTHOL  CARDS.  THE 
READ  BY  A  SPECIAL  SUBROUTINE  (  IREAU  )  ANO  ARE 
DETAIL  IN  THAT  SECTION. 


TAPE  REQUIREMENTS  — 

THE  PROGRAM  UTILIZES  UP  TO  THREE  TAPE  FILFS  IN  ADDITION  TO  THE 
NORMAL  INPUT,  OUTPUT,  AND  PUNCH  FILES.  THF  NUMBER  REQUIRED  DEPENDS 
ON  THE  OPTIONS  REQUESTED.  THE  TAPES  ARE  - 

A)  TAPE1  -  A  SCRATCH  TAPE  USED  TO  STORE  OBSERVATION  DATA  OUHING 

THE  ANALYSIS  (  USE  REQUIRED  ). 

B)  TAPE2  -  INPUT  TAPE  SUPPLYING  OBSERVATION  DATA  IUSF.  OPTIONAL). 

o  tape*  -  output  tape  containing  original  and  computed  values 

FOR  ALL  OBSERVATIONS  (  USE  OPTIONAL  ). 

THE  PROGRAM  ALSO  CAN  PRODUCE  PUNCHED  CARDS  FOR  THE  GRIDDED 
VALUES  OF  THE  TREND  SURFACES.  THESE  CAN  BE  SUBMITTED  AS  DATA  WITH 
PROGRAM  CONTUR  TO  PRODUCE  CONTOUR  MAPS  ON  THE  CALCOMP  PLOTTER. 


DESCRIPTION  OF  CONTROL  CARDS  -- 

THE  CONTROL  CARDS  ARE  DESCRIhED  IN  SUBROUTINE  IRfAO  ,  WHICH  READS 
ALL  OaTa  REQUIRED  BY  THIS  program  FROM  CAwDS  and  tape. 


2 

4 
6 

s 

-40 

12 

14 

-16 

ie 

20 
22 

24 
26 
28 
30 
32 
34 
36 
38 
40 
4? 
44 

-  46 
48 
50 

52 

54 
56 

-  58 
60 
6? 

_64 
66 
68 
70 
7? 
74 

-76 
78 
8  0 

-  82 
84 
86 
88 
90 
9? 
94 
96 
98 

100 
102 
104 
106 
108 

no 

112 

114 
116 
118 

120 
122 
124 
126 
128 
130 
132 
134 
136 
138 
140 
142 

144 
146 
148 

i  mi 
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2. 
3. 

4. 
5. 
6. 

7. 
8. 

9. 
10. 


11. 
12. 

13. 
14. 

15. 
16. 
17. 

lb. 

19. 
20. 

21. 

22. 

23. 

2*. 
2b. 


26. 
27. 
28. 
29. 
30. 
31. 


32. 


33. 
34. 
J5, 
36  ■ 
37. 
38. 

39. 

4  1. 
4  1  . 


REQUIRED  SUBROUTINES  - 

-  IREAD  -  READS  IN  DATA, 

-  R1TER  -  WRITES  MOST  OF  THE  OUTPUT-. 


.ANOVA  -  COMPUTES  AND  WRITES  STATISTICS. 

-  INVER  -  INVERTS  MATRICES  TO  COMPUTE  COEFFICIENTS. 

-  POLY  -  EVALUATES  EQUATIONS  AT  GIVEN -X-Y-COOR&INAT€S. 

OPTIONAL  SUBROUTINES  - 

-  PRTC0N4  -  PREPARES  CONTOUR  MAPS  OF  THE  -SURFACES.— — 

-  STATS  -  COMPUTES  STATISTICS  REQUIRED  BY  PRTC0N4  . 


DIMENSION  ISURF(5)i  IPUNI5) 

DIMENSION  TITL (8)  

DIMENSION  C(50.50> 

DIMENSION  X(IOOO).  Y(IOOO).  Z(IOOO) 

OIMENSION  XXI21.21).  A(2l,5).  XINV (21 , 21 ,5) .  PAR(21»,  B8(21).  LlMl 
1T(5) 

DIMENSION  SYMTABI21) 

DIMENSION  AREA (3) i  FACT0R13) 

DIMENSION  C0N(5)»  T0UR(5).  LINESI5).  INCHESI5).  ZMIN(5)|  ZMAX(5) 

COMMON  /INFO/  KNUMt JNUMiLlMI TiNCAROS. ARE A .FACTOR t XMAX. XMIN. YMAX. YM 
1 1 N. ROWS » COLS » I NCK.NRES.NnCARD.KTR. I  SURF, I  PUN .CON. TOUR, LINES. INCHES 
2.ZMIN.ZMAX.ITYPE.NTAPE.K 

COMMON  /MATRIX/  X.Y.Z.A.XX.XINV.BB 

COMMON  /MAPS/  K0N.NT0UR.LINE.*IDTK,lHYPS0.INCON(8) .URANG (25) »Z»ilNI 
l.ZMAXI 

INTEGER  SYSIN,SYS0T,SDF1,SDF2.SDF* 

INTEGER  TITL  

INTEGER  SYmTAB 

INTEGER  ROwS.COLS 

INTEGER  CON.  TOUR  —  — 

HEAL  INCHES 

DATA  SYSIN.SYS0T.SDF1.SDF2.SDF4/5.6.1.2.4/ 

DATA  TITLUI/lOHCONTOUR  MA/, TITL (2) / 10HP  OF  TREND/ • T I TL ( 3) / 10H  SUR 
1FACE  F/,TITL(«)/10HOH  DEGREE  / .  (  T  ITL (J) . J»S» 8) /l OH  / 

DATA  (SYMTAB(I) .  I»l  .5 > /l OHONE  .lOHTWO  .10HTHREE  .1 
10HFOUR       .10HFIVE       / 

EQUIVALENCE  (X.C) 

REWIND  SDF2  -  - 

REWIND  SOF4 
REWIND  SDF1 

INITIALIZE  ALL  COUNTERS 


IHYPSO»0 
XD"0.0 
YD'0.0 
2D»ft.0 
ZSUM"0.0 
SSZaO.O 
c        ..... ..............       

C    CALL  SUBROUTINE  IREAD  TO  OBTAIN  CONTROL  PARAMETERS 

c        .. 

CALL  IREAD  (1) 
C 
c        .................................... 

C    FIND  LARGEST  DEGREE  FIT  REQUESTED 

C    . . 

DO  15  laltS 

J«6-I 

IF  (LlMIT(J)l  10.15,10 
10        N»( (J*1)«(J*2) )/2 

JUMPs I 

JACK«b-JUMP 

GO  TO  20 


CONTINUE 
GO  TO  2U0 


15? 
154 
15* 
158 

-169 
162 
16* 

-166 
168 
170 
172 
174 
176 
178 
180 
182 
184 
186 
188 
190 
192 
194 
196 
198 
?00 
202 

?04 
206 
208 
210 
?12 
?14 
?16 
?18 
220 
?22 
224 
?2(S 
228 
230 
?32 
?34 
?36 
?38 
2*0 

?4? 
244 

?46 
248 
?50 
?5? 

?54 

256 

?58 

?bO 

*62 

?64 
?66 
268 
?70 
?7? 
?T4 
?76 
278 
280 
?8? 
?84 
?86 
?88 
?90 
292 

294 
296 
MA 
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♦  2. 
13. 

44. 
45. 
46. 


♦  7. 


48. 
49. 
50. 
51. 
52. 

53. 
54. 
55. 
56. 
57. 
58. 
59. 
60. 


61. 
62. 
63. 
64. 
65. 
66. 
67. 
68. 
69. 
70. 
71. 
72. 
73. 
74. 
75. 
76. 
77. 
78. 
79. 
80. 
81. 
82. 
83. 
84. 
85. 
86. 
87. 
88. 
89. 
90. 
91. 


92. 

93. 

94. 


C    ......... — — ........... 

c   initialize  parui  and  zero  xx  and  bb  arrays 
c   ................ — — — 

20     PAR(1)«1,0 

DO  30.  I«liN 
00  25  J«I«N 
25       xx(l,j)=o.O 
30     BB(I>»0.0 

C    .....-.....—.—.— ........... 

C    PREPARE  JOB  TITLE  PAGE 

C 

CALL  RITER  (JACK.V) 

C — — 

C    READ  DATA  FROM  TAPE3 

C 

IF  (NCAADS.GT.1000)  GO  TO  35 

KNUM»1 

KTR«NCAKDS 

KTSW"1 

GO  TO  40 


C 
35 


45 
C 

c 
c 

c 


50 


55 


60 


65 


C 

C 

c 
c 

70 
75 


KNUM»NCARDS/1000 
JNUMaNCARDS-KNUMMOOO 
KTR»1000 
KTSw«0 

NNCARD=NCARDS 
NCARDS»0 
DO  85  KaliKNUM 
CALL  IREAD  (2) 


CHECK  EACH  DATA  OBSERVATION  TO  SEE  IF  IT 
FALLS  WITHIN  LIMITS  OF  AREA  OF  INTEREST 

DO  85  IKT»1»KTR 

IF  (X(IKT) .GT.XMAX.OR.X(IKT) .LT.XMIN)  GO  TO  85 

IF  (Y(IKT) .GT.YMAX.OR.Y(IKT) .LT.YMIN)  GO  TO  85 

NCARDS«NCARDS*1 

XDaX(IKT) 

YD»V(IKT) 

ZOaZ(IKT) 

PAR<2)«XD 

PAR(3)«YD 

GO  TO  (50.50.50.50i70) i  JUMP 

PAR(4)*XD*XD 

PAR(5)»XD*YD 

PAR(6)»YD*YD 

GO  TO  (55.55. 55.70.70) .  JUMP 

PAR(7)«PAR(4)«X0 

PAR(8)*PAR(4)»YD 

PAR(9)»PAR(6)*XD 

PAR(10)"PAR(6)»YD 

00  TO  (60. 60.70. 70.  70)  .  JUMP 

PAR(11)«PAR(7)*XD 

PAR(12)«PAR(7)«YD 

PAR(13)«PAR(4)*PAR(6) 

PAR(14)«PAR(10)*XO 

PAR(15)«PAR(10)«YD 

00  TO  (65.70.70.70.70).  JUMP 

PAR(16)»PAR(11)«XD 

PAR(17)«PAR(H)*Y0 

PAR(18)«PAR(7)«PAR(6) 

PAR(19)aPAR(10)«PAR(4) 

PAR(20)»PAR(15)«XO 

PAR(21)«PAR(15)*YD 


DEVELOP  XX  MATRIX.  -• 

IT  CONTAINS  THE  XY-TRANSPOSE.XY  VALUES 

00  75  I"liN 

DO  75  JaltN 
XX(I,J)«XX(I,J)*PAR(I)»PAR(J) 


348 
750 
352 
354 
356 
358 
360 
3ft? 
364 
366 
368 
370 
372 
374 
376 
378 
380 
382 
384 
386 
388 
390 
392 
394 
396 
398 
400 
402 
404 
406 
408 
410 
412 
♦  14 
416 
418 
420 

422 
424 

426 
428 
430 
43? 

434 

436 

438 

440 
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9b. 

9b. 

97. 

98. 

99. 
100. 
101. 
102. 
103. 

104. 
105. 
106. 
107. 

106. 


109, 


110. 
111. 
112. 
113. 
11*. 


115. 
116. 
117. 

lie. 

119, 
120. 
121. 
122. 
123. 
12*. 
125. 


126. 
127. 
128. 
129. 
130. 


131. 
132. 
133. 
13*. 
135. 
136. 
137. 
138. 


13V. 

l*u. 


C 
C 
C 

c 

80 


85 


DEVELOP  BB  MATRIX.  — ■ 

IT  CONTAINS  THE  XY-TRANSPOSE.Z  VALUES 

00  80  I»ltN 
B8(I)-BB<I>«ZD«PAR(I) 
SSZ«SSZ*ZD«ZD 
ZSUMoZSUM*ZD 
CONTINUE 

IF  (KTSw.EQ.lt  60  TO  90 
KTHaJNUM 
KTSW»1 
GO  TO  45 
C 

90     00  100  11=1.5 
K»6-II 

IF  (LIMIT(K))  95i100i95 
95        GO  TO  105 
C 
100    CONTINUE 

C   ..................................... 

C    CALL  RITER  TO  PRINT  UPPER  RIGHT  DIAGONAL"  or'THE  XX  ANHA> 
C    .................................... 

105    IF  (INCK.E0.2)  CALL  RITER  (Ki2) 

C    COMPLETE  DEVELOPMENT  OF  SYMMETRIC  XX  ARRAY 
C    ................ ..... 

NM1«N-1 

00  110  I»1,NM1 
IP1»I*1 

00  110  J»IP1»N 
110  XX(J,I)aXX(I,J) 
c    ................... ................. 

C    CALL  RITER  TO  PRINT  ENTIRE  SYMMETRIC  XX  ARRAY 
C    .................................... 

IF  (INCK.EQ.2)  CALL  RITER  <K,3) 
DO  125  K"l,5 

IF  (LIMIT(K))  115,125,115 
115      NTERMS»((K*l)»{K*2))/2 

KSOBaO 

CALL  INVER  (NTERMS.K.KSOB) 

IF  (KSOB)  120,125,120 

CALL  RITER  (K.6) 

LIMIT(K)»0 
CONTINUE 
IF  (INCK)  130.1*0.130 


120 

125 

C 
C 

c 

130 


135 
1*0 

C 

c 
c 


145 


ISO 

155 

c 

c 

c 


CALL  RITER  TO  PRINT  IDENTITY  MATRICES 

CALL  RI TER  (JACK. 8) 
DO  140  K»l,5 

IF  (LIMIT(K))  135,1*0,135 

CALL  RITER  (K.l) 
CONTINUE 


COMPUTE  COEFFICIENTS 


DO  155  K»l,5 

IF  (LIMIT (K) )  145,155,1*5 
NTERMS«M  (K*1)«(K*2) >/2 

DO  150  I"1,NTERMS 

A(I,K)*0.0 
DO  150  J«1»NTERMS 
A<I.K)»A(I.K)*XINV(1«J.K)»BB(J> 

CONTINUE 


CALLSU8R0UTINE  ANOVA  TO  COMPUTE  AND  L 1 5T"  STATISTICS 

CALL  ANOVA  (SSZ.ZSUM, JUMP) 
IF  (NRES)  160.165,160 


♦  42 

44* 
446 
448 
450 
45? 
454 
456 
458 
460 
462 
464 
466 
468 
470 
472 
*7* 
476 
478 
480 
482 
484 
486 
488 
490 
492 
494 
496 
498 
500 
502 
504 
=.06 
508 
510 
*12 
514 
=■16 
518 
=■20 
522 
52* 
526 
528 
530 
532 
S34 
536 
538 

=■40 

5*2 
5** 
5*6 
5*8 
550 
=,5? 
554 
556 
55H 
560 
56? 
564 
566 
568 
S70 
■•72 
574 
=•76 
578 
«,S0 
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1*1. 

142. 
143. 
144. 
145. 
146. 
147. 


14B. 
149. 
150. 
151. 
152. 
153. 
154. 
155. 
156. 
157. 
158, 
159. 
160. 
161. 
162. 
163. 
164. 
165. 
166, 
167. 
168, 
169. 
170. 
171. 
172, 


173, 
174. 


175, 


c    ........ . .- 

C   OUTPUT  RESIDUALS  ON  TAPE  OR  JUST  PRINT  THEM 
C — 

160 


165 


170 
175 


180 
185 
190 

195 

C 

c 
c 
c 
200 


CALL  RITER  (JACK, 8) 

CALL  RITER  (0.4) 

CONTINUE 

DO  195  K»1.5 

IF  (ISURF(K) .NE.l.ANO.IPUN(K)  .NE.l)  GO  TO  190 

XINC» ( XMAX-XM IN) /FLOAT (COLS-1) 

Y I NC»(TMAX-YMIN) /FLOAT!  RO"S- 1) 

DEFINE  STARTING  X  AND  Y 
XDaXMlN 

YO»YMAX 

DO  175  JJ»l.RO*S 

DO  170  KKal.COLS 

CALL  POLY  (PT,XD»YD»K) 
C(JJ»KK)aPT 

XO«xD«XlNC 

X0"XMIN 
YDbYD-YINC 
CALL  RITER  (JACK, 8) 
CALL  RITER  (K,5) 

IF  (ISURF(K).NE.1.0R.C0LS.GT.64)  GO  TO  185 
TITL(5)=SYMTAB(K> 
DO  180  I»1,NMAP 

KONsCON(I) 

NTOURpTOUR(I) 

LINEoLINES(I) 

WIDTH«INCMES(I) 

ZMINI»ZMIN(I) 

ZMAXIaZMAX(I) 
CALL  PRTC0N4  (C i T I TLiROWS,COLS,50 ,50 ) 
CONTINUE 
CONTINUF 
CONTINUE 
GO  TO  5 


NO  FITS  ASKED  FOR.  CALL  RITER  TO  PRINT  STATEMENT. 


CALL  RITER  (0.7) 
GO  TO  5 


END 


58? 
584 
586 
588 
590 
592 
594 
596 
598 
600 
602 
604 
606 
608 
610 
612 
614 
616 
618 
620 
622 
624 
626 
628 
6.10 
632 
634 
636 
638 
640 
642 
644 
646 
648 
A5D 

654 
656 
658 
660 
662 
664 
666 
668 
670 
672- 
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p)  FORSTAT    CROSS-REFERENCE    LISTING    FOR    PROGRAM    TREND 


PROGRAM     TREND 

STATEMENT  NUMBER 

DEFINITION 

5 

25. 

10 

36. 

15 

40. 

2n 

42. 

25 

45. 

30 

46. 

35 

53. 

40 

57. 

M 

60. 

Sf 

71. 

*S 

75. 

*c 

80. 

*» 

86. 

10 

92. 

75 

94. 

en 

96. 

•S 

99. 

90 

104. 

95 

107. 

100 

108. 

105 

109. 

no 

114. 

115 

118. 

120 

122. 

125 

124. 

130 

126. 

135 

129. 

140 

130. 

145 

133. 

150 

137. 

155 

138. 

160 

141. 

165 

143. 

170 

154. 

175 

156. 

180 

168. 

1A5 

169. 

190 

170. 

195 

171. 

200 

173. 

REFERENCES 


172. 
35. 
33. 

39. 


174. 


35. 


43. 
4t. 

52. 
103. 

70. 
74. 

79. 

85. 

70. 

92. 

95. 

50. 
100. 
106. 
104. 
107. 

111.  113. 
117. 
121. 
116. 
125. 
128. 
125. 
132. 
134. 
131. 
140. 
140. 
151. 
150. 
161. 
159. 
145. 
144. 

41. 


74. 

93. 


61, 


106. 


-*3^~ 


117.      121. 


127. 


136. 
13?. 


128. 


40 

STATEMENT  NUMBERS 

DEFINED 

4 

MAXIMUM  NUMBER  OF 

REFERENCES 

57 

TOTAL    NUMBER  OF 

REFERENCES 

1  .4 

AVERAGE  NUMBER  OF 

REFEHENCES 
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q)    LISTING    OF    SUBROUTINE    IREAD 

1,  M'wHJII'it     [HEAD     (KUOF.)  H  2 

C        - - —  R  * 

C  SUBROUl  InE  IREaD  8  ft 

c H  fl 

C  PURPOSE--  H  in 

C  ro  PfmU  KM  DaTa  and  CONTROL  OPTIOnS  FROM  CaRDS  And  TAPE  AS  HE-  R  12 

6 SlHEtl.  H  14 

C  H  lft 

C  THE  GCmPS  SYSTEM  R  18 

e  __  .....__   pUKUU£  UNIVFrSITY  fl  20 

C  8  2? 

C  A.  K.  TURNER*  fl  2* 

C : CIvlL  ENGINEERING.  Ft  2ft 

C  JUNE  1969.  R  28 

C  R  30 

e      as-A-ee— < r  32 

C  CALL  IREAU(KOI'E)  8  34 

C  iHERF-KODE  =  A  PUlNTER  H)  CORRECT  ENTRY  PolNT.  H  3ft 

e— '  IF  KOOE  =  ll  CONTROL  CARDS  AND  DATA  TU  RE  REAP  IN.  R  38 

C  DATA  Tu  BE  PLACED  ON  TEMPORARY  STORAGE  8  40 

C  I  APE  (TAPE  1)  R  42 

C It-  *OUE  =  2»  HEAD  DATA  OFF  TEMPUHARt  STOWAGE  TAPF  8  44 

C  (TAPE  1)  8  4ft 

C  NOTE —  THIS  SUbPOUTINf  UTILUeS  THE  FoLLUWlNO  LAHFLLEU  COMMON  8  4H 

(V . /MAlWlA/  fl  50 

C  /InFO/  R  52 

C fl  54 

C  Ot  SCR  [PT  1 1 1Mb  OF  CONTHOL  fAi.US —  R  5ft 

C  8  58 

C  (  1  )  I.AOLE  CARD  -  8  ftO 

C  -■             IHF  WORD  ^TRENDS*  MUST  APPEAR  I  r-j  COLUMNS  1-6.  THIS  CAHD  R  62 

C  INDICATES  IriE  5TAR1  OF  A  DATA  SFT.  R  ft4 

C  R  ftft 

C  -  <2)  FASTER  COnIRuL  CAwD  -  R  68 

C  ("IS  CAPO  MUSI  HAVE  TmE  "OrD  *TPENOS*  IN  COLUMNS  5-10  .  8  70 

C  ThF  CARD  CAN  HAVE  EITHER  OF  "1  wO  FORMATS  -  8  72 

C" (A)  ALL  TRENDS  R  74 

C  OR  (b)      TRENDS  1»2»3«4»5  whERE  ANY  OF  THESE  NUribERS  MAY  BE  R  7h 

C  OMITTED.  B  78 

C  In  CASE  (A)  *ALL  TRENDS*  PEUUESTS  i  SET  OF  STANDARD  OUTPUT  8  80 

C  OPTIONS  AnD  The  FITTIm,  UF  ALL  FIVE  1  kEnU  SUPFACESlSEE  B  82 

C  Nt*l  SECTION  Flip  DETAILS).  In  C*SE  I")  FIVE  TPfNDS  WOULD  R  84 

C  -  «fc  FITTED,  rtUI  NO  OPTIONS  have  BEEN  SPECIF  UD.  CASE  (H)  R  8h 

C  CAN  ALSO  HE  USED  ID  RfcUUEST  ONLY  CfRlAlN  SURFACES.  8  88 

C  fl  9(1 

Ci T3)  "PIIUN  CAPU(S)  -  R  92 

C  THESE  ALL  HAVE  A  STAR  In  COLUMN  1  ,  FOLLOwEU  HY  A  NA"E.  R  94 

C  THf  PROGRAM  Only  CHr.Cr\S  ThE  FIRST  SI*  COLUMNS  OF  THE  CARD.  R  9ft 

C  (HF  REMAINDER  CAN  HAVE  COMMENTS  DESCRIBING  I  HE  REQUEST.  8  9H 

C  ThF  FOLLOWING  OPTION  CAPUS  CAN  PE  USFU  IN  ANY  COMB IN4 1  ION-  B  100 

C  ohAPS      -  PP1NI  DAPS  OF  FITTED  TREND  SURFACES,  fl  102 

C  »PUnCH     -  PUNCH  CAPDS  CUivlAlNlNu  (imID  VALUES  FOR  EACH  OF  fl  1  04 

C  The  FITTED  TREND  SURFACES.  fl  10ft 

C  »TAPE2     -  INPUT  OBSERVATION)  UAru  SUPPLltD  ON  A  UPF  B  108 

C  (lAPE2)iTHfc  PROGRAM  EXPECTS  3000-wOHO  LOGICAL  B  110 

C  RECORDS  WITH  IhE  DATA  STORED  Al  •  Y  1  ,Z)  ,X2« Y?»7?  8  112 

C  UP  lo   AlOOOtY  1000,?10n0.  ONE  LEAD  RFCOPO  8  114 

C  CONTAINS  IDENTIFICATION  InFORMAIIOm.  B  lift 

C  »lAPE4     -  STORE  ORIGINAL  V  ALOES,  CO'-  PUTED  VALdESi  And  8  118 

C  RESIDUALS  FOR  ALL  OBSERVED  DATA  POINTS  On  TAPfc  B  120 

C  I1APF4)  .  IHE  DAIA  IS  IF  THE  SAME  SEQUENCE  AS  B  122 

C  DISPLAYED  HY  IhE  'RESIDUALS  OPTION.  B  124 

C  "IDENTITY  -  PPINI  IDENTITY  MATRICES.  THIS  IS  A  USEFUL  fl  12ft 

C  -  CHP-Cts  UN  -IhE  ACCURACY   )F   IHE  NAIPIn  INVERSION,  B  12H 

C  »UEbUG     -PRInT  ALL  MATRICES  USED  IN  THE  COMPUTATION.  8  130 

C  «HES1DUALS-  PHINT  ORIGINAL  VALOFS.CO-'PUTEU  VALUES.  AND  B  132 

C  RESIDUALS  FOP  ALL  OhSEkvFU  DAIA  POINTS,  B  134 

C  THt SE  OPTION  CARDS  CAN  dE  USED  TO  INCREASE  IHE  OPTIONS  H  136 

C  SPECIFIED  HY  IhE  USE  OF  AN  *ALL  TRENDS*  maSIEP  CARD.  8  138 

C  wllHUUT  MODlFlCAl Ion,  THE  USE  OF  *«IL  TRENDS*  IMPLIES  THE  B  140 

C  USE  OF  «MrtPS  AimU  "IDENTITY  OPTION  capuS.  fl  ]42 

?  8  144 

C  (4)  »ElvD  CAPO  -  8  14ft 

C  IHF  OPTION  CARDS  MUST  HE  1ER-UNATEO  ivlTH  A  "END  CARD.  «En|i  8  148 

C  IS  PLACED  In  COLUMNS  1-4  .  IF  NO  DPI  I  UN  CARDS  FOLLOW  ThE  8  150 

C  MASTER  CONTROL  CaRDi  ThE  "END  CARD  FOLLOWS  IMMFDIATELY.  8  152 

C  In  OTHER  WORDS,  A  I  LEAST  l«u  CAPOS  ArE  REUUlPEn,  fl  154 

C  FOH  EXAMPLE  -  fl  15ft 

C  A|L  TRt.vDS  fl  )i;h 

c  »f:nu  h  iftn 
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C  H  Id? 

C  <S>  okEh  DEFINITION  CahO  -  H  164 

C  .  N     -ICOLS   1-  S>-  NUMBER  OF  DATA  PUlNTS  HtltvC,    SUPPLIED,  R  1  bb 

C  AMAA  -(COLS  11-20)-  MAXIMUM  X-COOROINaTE  OF  ANALYSIS  AREA.  R  168 

q          —  AI-)  1  **  -(COLS  21-30)-  MINIMUM  X-COOHOINATE  OF  ANALYSIS  AREA.  B  170 

C  YMAX  -(COLS  31-40)-  MAXIMUM  Y-COORDlNATE  OF  ANALYSIS  AREA.  H  172 

C  YMIN  -(COLS  41-50)-  MINIMUM  Y-COORdInaTE  OF  ANALYSIS  AREA,  B  174 

G __.._,  R  176 

C    — --- — — -  —  --- —  R  17H 

C  R  ISO 

C fne-FOLLOwlNft  CARWi  ARE  REQUIRED  ONLY  IF  THE  »"APS  OR  "PUNCH  OPTIONS  R  182 

C  HAVE  BEEN  SPECIFIED  (  ThEY  ARE  REQUIRED  WHENEVER  THE  fALL  TRENDS*  R  184 

C  MASTER  CONTROL  CARD  IS  USED).  H  186 

C- -  —  R  18R 

C  (6)  (,R1D  SPECIF  I C  A  r  1 0 1  ■  CaRU  -  R  190 

C  KOviS  -(COLS  1-3)  -   DESIRED  NUMBER  OF  HOWS  IN  GRID  OF  COM-  S  1 92 

C • -  PUTEO  1RENU  SURFACE  VALUES.  R  194 

C  COLS  -(COLS  4-6)  -   UESIREO  NUMBER  OF  CULUMNS  IN  GRID  OF  H  196 

C  COMPUTED  TREND  SURFaCF  VALUES.  R  198 

C  —  -  -     NMAP  -(COLS  7-9)  -   NUMBER  OF  COPIES  uF  MaPS  OF  EACH  SUR-  R  200 

C  FACE   (MAXlMUMaS) .THE  MAXIMUM  NUMUFR  OF  MAPS  R  202 

C  That  CAN  HE  PRODUCED  RY  ONE  ANALYSIS  IS  R  204 

C  •—                         2S  (  b    TRENDS  X  S  MAPS  FOR  EACH  TREND).  R  ?06 

C  NOTE  -  ROWS. COLS  •  AND  NMAP  ARE  INTEGERS.  R  ?08 

C  R  210 

C  (7)  PWO-T  CONTROL  CARD(S)  -  R  212 

C  NMAP  SETS  OF  IhESt  CARDS  ARE  REQUIRED.  EACH  SET  CONTAINS  -  R  214 

C  R  216 

C  (A)UNE  CARD  WIIH  -  B  21 R 

C  CON  (COLS  1-2)  -  NUMBER  OF  CONTOUR  INTERVALS  (MAXIMUM  =19)  R  ?20 

C  (CON  =  DATA  HANGE/DES1RED  CONTOUR  INTERVAL)  R  222 

C  FOUR  (CuL3l-    IF  TOURbOi  PROGRAM  USES  MAXIMUM  AnD  MJnIMUM  R  224 

C  OF  SURFACE  FOR  UPPER  AND  LOWER  CONTOUR  B  ?26 

C  LIMI IS.  B  ?2H 

C  IF  ToUH=l,  DESIRED  TOP  AND  HUTTOM  CONTOUHS  H  230 

C  READ  IN  FROM  NEXT  CARD  (SEE  bELOw),  R  232 

C  IF  TOUR=3. CONTOUR  I mTFHVAL  BECOMES  \/g  fl  ?34 

C  STANDARD  DEVIATION  (max  b  +3  SIGMA.  R  236 

C  MIN  a  -3  SIGMA  )  .  8  ?38 

C  NOTE  -  TUUR  =  2  (lHHEGDLjR  CO-lTnUR  INTERVALS  OPTION)  R  ?40 

C  CANNOT  HE  REUUESIED  RY  TREND.  R  242 

C  Llt«ES-(COL  4)  -IF  LINES  =  0,  CONTOUR  LINES  ARE  PRIMED.  R  244 

C  IF  LINES  =  If  ALTERNATE  (EVEN  VALUED)  R  246 

C  CONTOUR  HANDS  ARE  PRIMED.  B  248 

C  IF  LINES  =  2«  ALL  RAMOS  PRINTED  (ALL  PRINT  R  250 

C  POSITIONS  ) .  R  252 

c  inches-(cols-b) i.mches  is  width  of  map  in  inches. if  gheaiep  r  254 

C  THAN  12.7  UR  BLANK.  SET  TO  1U.0  )  R  256 

C  R  ?58 

C.  (Ri  AUDI  F  IonAL  CAuo-   PEOI-HEO  ONLY  IF  TO"h  b  i  b  3«,q 

c  This  card  specifies  the  minimum  and  maximum  contour  b  262 

C  ELEVATIONS.   <!m1N  (COLS   1-10)  MINIMUM  ELEVATION.  B  264 

C  -/MAX  (COLS  11-20)  MAXIMUM  ELEVATION.  R  ?66 

C  NOTE  -  DECIMAL  POINTS  SHOULD  RE  PUNCHED.  B  268 

C  B  270 

C  fl  ?7? 

C  R  ?74 

C  li'E  FOLLOWING  CARDS  AWE  REQUIRED  ONLY  WHEN  TmE  DATA  lb  HFlNG  R  776 

C  SUPPLIED  FROM  CARDS.  THEY  MUSI  RE  OMI1TEO  IF  A  »lAPE2  OPTION  CARD  R  278 

C  Is  USFO,  r  ?80 

C  P  782 

C  (B)FOKMAT  CARD  -  B  ?B4 

C  NH  -(COLS  \'b    >-  NUi-mtR  OF  OBSERVATION  SETS  /  CARD.  R  286 

C  EORM-(COLS  l-b^)-    n.mrwvT  uf  data  CARDS,  r  ?bh 

C  R  ?9  0 

i-  ('■•)  TITLE  CARD  -  R  ?97 

C  EaCIOR  -(COLS   1-30) -ALPHANUMERIC  IDENTIFICATION  OF  FACTOR  R  ?94 

C  AREA    -(COLS  <f  1-70)-«lPHaNUmEHIC  I DEmT IF  1  CAT  I  ON  OF  AREA.  B  ?96 

0  HYPE  -  (COL  78)     -IF  DATA  MEASURES  COS  I S    .IFYPE  =  2  H  ?98 

C  -IF  DA1A  MEASURES  HENEF 11 S, 11 YPk  i  3  R  "500 

<-  R  ■*(}? 

C  (  I  0)  DATA  CARDS  -  R  ?Q4 

C  a  SERIFS  Of  CARDS  R  306 

C  A  SERIES  OF  CARDS  CONlAlMNG  X.Y.  AND  ^-COORDINATES  FOR  R  108 

C  ALL  OBSERVATIONS.  R  qo 
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PURPOSE — 

TO  PRINT  TRFNC)  SUhFaCE  EQUATION  COEFFICIENTS  ANO  10  COMPUTE  AND 
PWINT  ANALYSIS  Of-  VAHiANCt  STATISTICS.   STATISTICS  CO'"PUlE(>  INCLUDE 

(1)  ERROR  ESTIMATES  FOR  COfFFICIFNTS 

(?)  IOTaL  K-TESTS  FOR  ALL  EQUATIONS 

(3)  PARTIAL  F-TESTS  FOR  Hl(,nFST  DEGREE  FCUaTIOn 

(4)  PE^Cf  NT  VARIATION  F*PLaINEO  (H-SOUftPEU) 

iiSflGt  — 

CAl L  ANOV/A  (SSZ.ZSUM.  JACM 
WHERE  -^SZ   =  UNCORRECTED  SUM  OF  SQUARES  OF  7 
ZSUM  =  SUM  OF  L 

JACK  =  HIGHEST  DEGREE  TRFNO  SURFACF  RFOUFSTEO 
NOTE-  ThIS  SUHKOUTInE  UT1LI7ES  THF  FOLLOWING  LABELLED  COMMON 
/IM   0/ 
/MATRIX/ 

OFSCHIPTlOO.  OF  CONTKOL  CaKOS— 

nOmF  PFoiJIHKi 

REQUIRED  SUPROUTINT S- 
NONE 

OPTIONAL  SUHWOUTTM- S- 
NONF 

DIMENSION  AHFA(T),  FACTUrt(i) 
DIMENSION   rSUWF  (S)  ,   IPIIN(S) 

DIMENSION    nUMJ(3000)»    0UM2(4411 

DIMENSION    |.lMlT(b)f    A(?ltb)i     bH(21>»    SSR<5)»    SSF(b)«    DFRFE(S>«    ERR 
1(21),    x  1 1\  v  (  ?  1  .  ?  1  .  b ) 
IN1EGER    SYSINiSYSOTtSOF ltSOF2»sOF4 

CU'^MfiN  /INFn/  KNlJM,JNUM»LHl  r,MCARllSiAREAtFACTnHtXMAX,XMIN,Y,"IAX,YM 
IINiROwS.COI.S.UCK  .NRES.NNCAHD,KlR,TSURF.IPUN,rOM,TOUR,L  INFS.lNCHf S 
2»2MTN,ZMAX.ITYPEtNTAPE,K 

COi-imon    /N-ATRIX/    DUiiJ  »A.UUM2»XINV,bR 

uata  sysin.sys0t.s0f1 ,suf2.sdf4/5t6,  ] .2.4/ 
calculate  Corrected  sun,  ok   squares 

ZbAH=Z5UH/FlOAT(rCAKUS) 
SSZHAH=2SUm»Z8aR 
S5Z=SSZ-SS7HAR 
JACK=6-JUMP 
...  r   g4 

PRINT  COEFFICIENTS  C   Ri 
. (-  q^ 

CALL  HIlER  (JACK.p)  C  100 

WHITF  (SYSOT.10S)  NCAH.iS  C  1  0^ 

WNITF  (SYSOTfllO)  C  104 

00  in  K  =  l  .S  r  1M 

TF  (LIMIT (K)  )  S,10,S  C  1  OH 

NTE»MS=( <K*1 )«(K+2> )/£  C  110 

WRITE  (SYS0T.136)  C  11? 

WRITE  fSYSOT.llb)  K,(A(I»K)  iIsl.NTERMS)  f  11* 

CONTINUE  C  lift 

WHITF  (SYSOTil35)  C  1 1 « 


?5 


FHHnR  FSIIMATES  FpH  CoKKFIClENTS 

CALL  RllER  (JACK.R) 
WRITE  (SYSOT,1?0) 
Do  30  K=1,S 

IF  (LIMTT(K))  IS, 30, lb 

NTE«MS= ((«♦!)» (k  +  2) ) /£ 

SSR(K) aO.O 

00  ?0  J=l ,NTFRMS 

SSR<K)=SSR<K)*AlJtK)»riH(J) 

SSR(K) =SSR (Kl-SSZHAP     , 

SSF(M  =SSZ«SSR (K) 

DFREE(K)*NCARDS-NTEPKS-1 

SIGMAF  =  SORT(SSE(M/nFHEF(K)  ) 

00  ?5  J=l  ,NTEf.MS 

FRR(J)=SlGMAF»Srjnt(xIN>/(J,J,K)) 

'-RITE  (SYS0T.13b) 

WRITE  (SYS0T,)2S)  K, (ERR ( J) , J»l .nTERMS) 
ru.~  TlMi'E 
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44 
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4M 

c 

SO 
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Si 
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SR 
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so 
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62 
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64 
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11 
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70 
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72 
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74 
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71 
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7H 
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HI 
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124 

12ft 
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WKITF  (SYS"T,I^S) 


ANALYSIS  OF  VARIANCE 


TOTAL  F-TtSTS 


41  . 
42. 
43. 
44. 
45. 
46. 

47. 
AB, 

49, 

50, 
51. 
52, 
53, 

^4  . 

55. 
156, 
57, 
58. 
59, 
60, 
61, 


6?. 

63. 
64, 
65, 
66. 

67. 
68. 
69, 

7  1). 
71  . 
7?. 
73. 
74. 
75. 
76, 
77. 
78, 
79. 

HO. 
SI  . 

R?. 
R  J. 

R4, 
R5, 
96, 
87, 
RH. 
R9, 
90, 
91, 


42. 
93, 
9*, 
95, 

96, 
9  7  . 
9H  . 
99, 
100. 

Inl. 
In2, 
10  3, 

If  4. 

IPS. 


4S 

50 


60 
C 
C 
C 
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ro 


CALL  RTTER  <JACK,t<) 
WHITE  <SYSOTil30> 
WKITF  (SYS0T»j35l 

WR1TF  <5YSOT,140) 
WKITF  (SrSOT,135) 
MjF=nca*0S-1 

WHITP   (STSOTiTS)   KlidSS/USR 
WRITF  (SYSOT.135) 
DO  4(1  K=l»5 

IF  (LIMtT(K) )  3b, 40, Jb 

NTFRMSx(  (K*l )*(K*2> )/£ 

XSEsSSEIKl/OFREEIK) 

ASW=SSR(K)/KLOAT(NTFHMS) 

XSE=SSE(K) /OFREF (K) 

XSR=SSR(K) /FLOAl UvThRMS) 

FxXVk/XSF 

uOFsOFPEE(K) 

WRITE  (SYS0I«)4b>  K,NTERMs,nDF»SSR(K) .XSRtSSFfK) »XSF«F 

WRITE  (SYS0T,13S) 
CONTINUE 
WrI TF  (SYSOT .80)  SS^ 


PARI  ML  F-TF^TS 


CALL  RlTFJR  (JACKia) 
WHITE  (SYSOT.KO)  JACK 
WHlTF  (SYSOT,155) 
wRITF  (SYSOT,l-sn) 
wkiTF  (SYSOtjlSS) 
Nl'F  JxOFKEF  (  JACK) 
Kh-T=l 
DO  6n  K«]i5 

IF  (LlMpiM)  4",, 60, 45 

NTFWMS=(  <K*  1  )»(K»2)  )/2 

(JO  TO  <<Ml,bK)  ,  KSwT 

KSWT=2 

NTf*»mS  =  nTF  RwS-1 

XMSO=SSP(K)/FI OAT(NTtHMS) 

F  =  XMSO/(SSK  ( JACK) /DFkEEl JACK)  ) 

WRITE  (SYSOI,16b)  K.NTEHMStNDFJ.SSK (K) tXMSQ,F 

WRITt  (=;YS0T,l5b) 

M TEaNTERMS+l 

J»K 

GO    TO    6n 

DFsMTEPmS-NTE 
SS1=SSR(k) -SSH(J) 
XMSQ=SSO/DF 

F=XmS0/(5SF ( JACK) /OFKtt( Jack) ) 
NOFsilf 

WRITE     (CYSOT.  )&b)     K,Nl)F,^UFj,SSo,XMSO,F 
WRITE     (SYSOT.lbb) 
J=< 

MTF=nTFRHS 
CONTINUE 


RFHCtNT  VARIATION  FXFLAlNtd 

CALL  RlTER  (JACK, a) 
wm I TF  (SYSOTtBS) 
WHITE  (SY50Tt9f)l 
WHITE  (SYSOT.9S) 
WHITE  (Srsoi,on( 

UO  7  0  KX  1  ,r, 

IF  (LIMIT!*))  b^,70,ob 

HR=SSR(K)/Sq? 

PCENTaRfl'lOO.O 

R=SOH1 (RR) 

WRITE  IPYSOTilOO)  K.t'CEMT  ,RR.  R 

»nire    (<;ysot»vo> 
CONTl  ^l)E 

REH.WM 
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InH, 

109, 
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117. 

118. 

II1*. 
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1        I 

2(1"!  , 

FOKM4 

lUiu.h 

FOKMA 

INT    si) 

2H-S,lH 

3 1 E  N  T . 

F  UK'-"> 

1  ) 
FOKMA 

lPt"CF 

2  1  i? 
FOkMA 

i    .Fin 

FOKMA 
IS     IS 

FOkMA 

)NL)     Sll 

2$X«'Dn 

3A7(X« 

4*A12< 

s««4>* 

ft»»3) ♦ 

FOKMA 

FOKMA 

IaTf     n 

FOwMA 

FOKMA 

1  F-TF 
FORMA 

1 

FOKMA 
IE       I 

2Rok 
llrtl.l 

FORMA 
1I1.3H 
2A,2(l 

FOKMA 
1AL  F- 
20  EQU 

FOKMA 
1 

FOKMA 
1GHOUP 

2  F-K 
3l»/t2 

FOKMA 
)  HE 
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END 


r  ( \  m  ,«*,, 

1 , » I4.3H  1 
1  -I  X  )  ,  \  H  I  ,  1 

r  (///3ox. 
i 

T  (1M0.54X 
Ms  OF  SQUA 
A^tO,  ThE 

T  I 1H  22X. 


2(lhl,lnX)  ,4  1  lMlilBX)  ,  1      r.-l?A,lMT,/,SX.  1KHI  MFAIM 
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r  (  lHO>A0X,t,2hlhe  imUMhlK  OF  ORSFH^at  [ii.-iS  USFO  U  IhIS  AMALYSI 
.IS) 

T  (  lH0.37x,(-,7HlHF  COMPUTEn  COEFFICIENTS  FOO  ThF  RF.OUESTED  TKE 
^FmCF  FQUaI  IUi.S»//.4 JA«30hThE  EQUATIONS  HaVF  THE  form  — .//.4 
HFST.7  =  A1*«21A) ««J(Y) ♦A41X«»2) *AS(X) (Y) .Aft ( Y««2> *»/.53X.47H 
«-))*Afl(X»«?)  (  Y  )  »AW  (  X)  (Yo»?)»AlO(Y»»3)»  i/i'ilXi  +  THJU  |K»«»I 
X»»3) (Y)*A13(A*»2) <Y*»*) ♦  ./.S3A.47HA14(X)(Y*»3)*A15(Y 

AlS(X««5)+A]7(X*»i,>  IY)  ♦   ./•S3X.47HAim^o«1)  (Y»«?)*A19(X«»?)  (Y 
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T     (lHntb4X.3lMAiMAI.YSIS    OF     VAKIAnCF     ST  a  1  [  ST  T  C  S  ,  //  »  4ft  X  .  4AHE  S  T  I  M 

F    tHROR    IN    1Mb     POi  Yi\0*1aL    COEFF ICIFnTS.//,  1  OH  DEGREE) 

T     (1H0.I6.1X. ?El6,2t//E24,2,6Fl6,?,//F24,?,f>Elb.2) 

T     ( lH0.b4X.31hANAI.YSIS    OF    VARIANCE    ST  A  I  I S TT CS , // , S2X » 3ftHT OT AL 

STS    FOR    ALL    TRt.NO    SURFACES.//) 

T     <1H    .HX.H2M 


T  (1H  .BX, 
0  .  F  .  1 
I  MS  nuE  T 
2X.1HI) 
T  (1H  .BX. 
T  .I?.]H, 
HI  . 1  OX)  .4 ( 
T  <lHf),54X 
TF5IS  FOR 
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T  >1H  »20X 


2 (lhl .loX) 

SS  DUF  I 

O  EKKOH   1 


) 
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33. 


36, 
37, 

3H, 


LISTING    OF    SUBROUTINE    RITER 

su.iwoutine   piter    iKisnnt)  n        2 

c       ' n       4 

C                         SUBROUTINE  wITEw              -  D    6 

c n   r 

C    PURPOSE —                                    «  D   1» 

C          Tn  PRINT  ANSWERS  ACCORUInG  TO  REQUESTED  OPTIONS  D   12 

C    D   1* 

C    HSAbE —  0   16 

C                CALL  RlTERUtKODt)  P   \H 
C          fJHEWE-K      =  DEGREE  OF  E'JUATION  FOR  WHICH  ANSWER  IS  BFlNG  PRO-    p   20 

C                       t'UCtO  p   ?2 

C                KOPE  s  POINTER  US>-u  To  SELECT  PROPER  RFOUEST  ACCORDING  TO   D   2* 

C                       COMROI  OPTIONS  IN  CALLING  ROUTINES  D   26 
C         NOTE-  THIS  SUBROUTINE  uriLl7ES  THF  FOLLOWING  LABELLED  COMMON-     P   28 

C                               /INFO/  D   30 

C                             /maIkix/  p   32 

C    p   34 

C    PFSCkIPTTONS  OF  C'.lNlROI.  CARL'S  -  D   36 

C         NONE'  REQU1RFP  p   3B 

C    p   40 

C    RFQUIREl)  SUBROUTINES-  P   42 

C          POLY  P   44 

C    OPTIONAL  SUBROUTINES-  0   4ft 

C          NONF  P   48 

C    'p   5n 

DIMENSION  »RFA<3)»  FaCTOR(3>  D   S2 

DIMENSION  ISURF15)"  IPUNCj)  p   S» 
DIMENSION  X(inno).  Y(loUO).  ^(1000).  A(21.S),  XX(21<21)|  XINV(2]i2   P   S6 

Ili5)  p   SB 

DIMENSION  RR(2l),  LIMIT(b),  BHR(21),  ZR<10)»  C(=.0»50)  D   60 

EQUIVALENCE  (X,C)  P   ft2 

RtAL  INCHES  p   64 

INTEGER  SYSIN.SYS0T»SUE1 »SoE2»SDF4  D   66 

DATA  SYSIN,SYS0T,SDEIiSUE2.SpF4/5i.ft.l .2.4/  p   ftH 

INTEGER  ROWS.COI.S  P   70 
COMMON  /INFO/  KNUM. JNUm . L Im  1  T » NCARPS . ARFA « FACTOR , XmaX , x MI N, Ymax , YM   p   72 
1  IN. ROWS. COLS. INrK,NK£S.NNCA«U.KTR, I  SURF . I PIIN.CON. TOUR ,l INES, INCHFS   P   74 

?»/:mIN,7MAX,  ITYPE.NTAPE.KT  p   7ft 

COMMON  /MATRIX/  X.Y»2.A»XX,XINViBB  p   7H 

N=l  (K*1)«(K*2)  )/2  P   B,) 

GO  To  <S,2o,35.45,120,l4b,lSo,lt>5»160)  ,  KODE  P   82 

C    . n   R4 

C    KODE  =1  p   Bft 

C    PRIN1  IOENTITY  MATRICES  P   Rm 

C    n   <,n 

5       WRITE  (SYS0T.19S)  K  p   92 

WRITF  (SYS0T.190)  p   94 

DO  IS  1=1, N  p   9ft 

no  in  j=i,n  0  RR 

RB8(J)=0.0  P  100 

no  in  L«1»N  p  in2 

10         BBB(J)=XX  (I  ,L)AxINV(LtJ.N)*HBB(J)  P  104 

IS     WRITF  (SYSOT.200)  1  .  (BHM ( J)  . J= 1  ,N)  p  10ft 

RETURN  n  10H 

c   ,,  no 

C    KOOE  =2  D  112 

C    PRINT  UPPER  TRIANGLE  OE  XX  MAlRIX  R  JJ4 

C    PRINT  THF  HR  MATRIX  p  \\(, 

f    D  118 

20     WHITE  (SYSOT.2SS)  p  120 

WRITF  (SYS0T.27S)  AREA, FACTOR, K  p  \?? 

WHITE  (SYSoT.lftS)  P  !24 

00  2S  1*1, N  0  12ft 

WPITt  (SYSOT,170)  I, (XX(I.J) ,J=] ,N)  p  i?„ 

2S     CONTINUE  p  no 

WRITF  (SYSOT.17S)  D  132 

DO  30  1  =  1, N  0  j34 

WRITE  (SYSOT.IBO)  I,BB(I)  P  j36 

30    Continue        „  p  J3H 

c      RETURN  0  140 

c  -"=" ;  \:i 

C    PRINT  THF  COMPLETE  XX  MATHIX  0  U6 

35     WRITE  (SYS0T.18S)  2  !*0 

oo  40  r-i.N  °     *° 

WRITE  (SYSOT.170)  I , < XX ( 1 , J) , J= 1 .N)  n  [54 

40     CONTINUE  n   g* 

RETURN  !?  J?" 

0  lftO 
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40. 
41. 
42. 

43. 
44. 
45. 
46. 
47. 
46. 
49. 
50. 
51. 
52. 

53. 
54. 
55. 
56. 
57. 
58. 
59. 
60. 
M. 
62. 
63. 

6*. 
65. 
66. 
67. 
68. 
69. 
70. 
71. 
72. 
73. 
74. 
75. 
76. 
77. 
78. 
79. 

80. 
Ell. 


82. 
83. 
»4. 
85. 
86. 
87. 
88. 
84. 
90. 
91. 
92. 
93. 
94. 


95. 
46. 


97, 


C 

C 

C 

C 

45 

50 


55 

60 


65 


C 
70 


75 
80 


85 


95 
100 


105 
110 


C 
115 


KOOE  =  4 

PRINT  ALL  RESIDUALS  FOR  OBSERVATION  LOCATIONS 

GO  TO  (50.60.50) i  NkES 

WHITE  (SYSOT.240) 

WRITE  (SYSOTi245) 

OO  55  1=1.10 

ZR(I)=0.0 

IF  (NRES.EQ.l)  GO  TO  65 

WRITE  (SDF4)  NC*RDS 

REWIND  SDFl 

IF  (NNCARD.GT.1000)  GO  TO  7U 

KNUMsl 

KTRbNNCARD 

KTSW=1 

GO    TO    75 

KTR=moo 

KTSW=0 

KPT  =  1 

00  110  Kel.KNUM 

READ  (SDFl)  (X(I) tY(I) tZ(I) iI'lfKTR) 
DO  110  1=1. KTR  _ 

IF  <X(I)  .GT.XMAX.OR.XU)  .LT.XMIN)  Qo  TO  110 

IF  (Y<I) .GT.YMAX.OR.Y(I) .LT.TMIN)  GO  TO  110 

XD  =  X< I) 

YO=Y ( I) 

ZD=Z(I> 

DO  90  KO=l  .5 

IF  (LIMIT(KO))  85t90.t)5 
CALL  POLY  (PT.XD.YO.KD) 
ZR(2#KO»] )=PT 
ZR(2«KOI=?D-pT 

CONTINUE 

GO  TO  (95, 100.^5)  .  NRES 

WRITE  (SYSOT.250)  KPTiX(l) iY(I) iZ(I) i (ZR(KK) .KK«1«10) 

IF  (NRES.EQ.l)  60  TO  105 

WRITE  <SPF4>  KPT»X<I>  .Yd)  .Zd>  »  (ZR(KK)  .KK«1.10> 

KPT=KPT«1 
CONTINUE 

IF  (KTSW.EO.l)  GO  TO  lib 
KlRaJNUM 
KTSW=1 
GO  TO  80 

CONTINUE 
RETURN 


KODE  a  S 

PRINT  AND  PUNCH  GRID  VALUES 


C 
C 
C 

c       - 

120    WRITE  (SYS0T.205)  K, COLS, ROW? 

WRITE  (SYSOT.210)  XMAX , XMIN i YMAX , YMIN 

DO  125  JJ"l»ROWS 

WRITE  (SYS0T.215)  (C(JJ.^K) ,KK«l,COLS) 
125    WHITF  (SYS0T.220) 

IF  UPUN(K))  130.140.130 
130    WRITE  (SYSOT.230) 

PUNCH  ?35.  FAC10R.K. AREA, ROWS. COLS. ITYPE 

00  135  J»l,ROWS 

PUNCH  225,  <C(  J.KK)  ,KK  =  1.C0|.S) 
135    CONTINUE 
140    CONTINUE 

RETURN 
C 

c      kodl  =  6 

C    THIS  MATRIX  SINGULAR. PRINT  STATEMENT 
C    — 

145    WHITE  (SYSOT.260)  K 

RETURN 

C    .. ... ... . 

C     KODE  a  7 

C    MO  THFNOS  REQUESTED, PRINT  STAlEMFNT 

C    ... 

150    WRITE  (SYS0T.265) 

RETURN 


D 

162 

D 

164 

D 

166 

n 

168 

□ 

170 

D 

172 

D 

174 

D 

176 

D 

178 

D 

180 

D 

182 

D 

184 

D 

186 

D 

188 

D 

190 

0 

192 

0 

194 

D 

196 

n 

196 

D 

200 

D 

202 

0 

204 

D 

206 

D 

208 

D 

210 

D 

212 

D 

214 

D 

216 

0 

218 

D 

220 

D 

222 

D 

224 

D 

226 

D 

228 

D 

230 

D 

232 

D 

234 

n 

236 

D 

238 

D 

240 

D 

242 

D 

244 

0 

246 

D 

?46 

0 

250 

P 

252 

D 

254 

D 

256 

D 

258 

0 

260 

D 

262 

n 

264 

D 

266 

P 

268 

D 

270 

D 

272 

D 

274 

P 

276 

D 

27H 

P 

?ao 

P 

282 

P 

284 

n 

286 

0 

286 

n 

290 

o 

292 

D 

29* 

P 

246 

P 

298 

[) 

3O0 

P 

302 

P 

304 

P 

306 

P 

308 

P 

310 

P 

312 

P 

314 
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99. 
100. 
101. 


102. 

160 

103. 

104. 

105. 

C 

c 

106. 

165 

107. 

170 

108, 

175 

109, 

180 

110. 

185 

111' 


112. 

113. 
11*. 

115. 

116. 
117. 
118. 
119. 

120. 
121. 

122. 


123. 
12*. 


125. 

126. 
127. 


128. 

129. 

nn>FOF**« 


c    

C    KOHL  =  8 

C    PRINT  GCAHS  PAGE  HEADING 

C    — ---• 

155    WRITE  (SYS0T.255) 

WHITE  (SYS0T.275)  AREA . FACTOR ,K 

RETURN 

C 

C    K0DE=9 

C    PRINT  JOB  TITLE  PAGE 

c — ........ 

WRITE  (SYS0T.255) 

WRITE  (SYSOTi270) 

WHITE  (SYStlT.275)  ARE  A.  FACTOR,  K 

RETURN 


FORMAT  <1HO»*7X,35HPRINT.OUT  OF  ALL  MATRICES  REQUESTED, /30X ,58HTHE 
1  UPPFR  TRIANGULAR  PORTION  OF  THE  XY-TRANSPOSE «XY  MATRIX,//, 4H  ROW) 

FORMAT  (1H0,I3,2X,11E10.3,/6X,10E10.3) 

FOHMAT  <1H0»////»52X»25HTHF.  XY-TRANSPOSE ,Z  MATRIX ,/// ,4H  ROW) 

FORMAT  (1H0,I3,2X,E16.4)  ~  . 

FORMAT  (1H0»////«43X,45HTHF  COMPLETE  SYMMETRIC  XY-TRANSPOSE. X Y  MAT 
1RIX,///»4H  ROW) 

FORMAT  (1H0,52HTHIS  MATRIX  SHOULD  (IDEALLY)  BE  AN  IDENTITY  MATRIX. 
1  i/,lH  »55HIF  IT  DEPARTS  GREATLY  FROM  THIS  IDEAL.THE  CORRESPONOING 
2,/»lH  .44HDEGREE  TREND  SURFACE  SHOULD  PE  HELD  SUSPECT,,///) 

FORMAT  HHOiBlH  THE  PROUUCT  OF  THE  XY-TRANSPOSE  XY  MATRIX  AND  ITS 
1INVFRSF  FOR  THE  FIT  OF  DtGHEE  ,12»///4H  ROW) 

FORMAT  (1HO,I3,2X,11F10.7,/6X,10F10.7) 

FORMAT  (1H0O6HGRI0UED  VALUES  OF  SURFACE  OF  DEGREE  ,I1»9H.  GRIO  IS 
1»I3»12H  COLUMNS  PY  ,I3,bH  HOWS.) 

FORMAT  (1H0,55HMAXIMUM  ANO  MINIMUM  COORDINATE  VALUES  READ  IN  ArE  X 
1MAX=,F10.3,6H  XMIN=»F10.3,6H  YMAX=  .  F  1 0 . 3 • 6H  YMTN= ,F 10 .3//) 

FORMAT  (1H  .12F10.3) 

FORMAT  (1H0) 

FOHMAT  (8F10.3) 

FORMAT  (1H0,67H*«»  PUNCH  OPTION  REQUESTED  **•  THIS  MATRIX  BEING  P 
IUnCHFO  ON  CARDS.) 

FORMAT  (3A10.3A10.8H  DEGREE  ,  1 1 ,  1X.2I 3, 12) 

FORMAT  (lH0,l8X,f«2HCOORDINAlES, ORIGINAL  AND  COMPUTED  VALUES, AND  RE 
1SIOUALS  FOR  ALL  SURFACES  REUUESTEO. ,//) 

FORMAT  (121H0  ID      X  Y  Z 

1Z-QUAD.   RES-QD   Z-CUBIC   RES-C   Z-QUAR. 
2T/1X) 

FORMAT  (1H  ,I6,2F8,2,F9,3,5(F9.3,F9.3) ) 

FORMAT  (1H1,55X,29H ——--..  —  . ....  .  ,/,56X  ,29HI       T 

1HE  GCARS  SYSTEM      I,/»56X»29Hl  PURDUE  UNIVERSITY  VERSION  I,/,56X 

2,29H— — — — —  ,///,40X,62HSFARCH  ACTIVITY  -  DAT 

3A  ANALYSIS  PHASE  -  TREND  SURFACE  ANALYSIS,//) 

FOHMAT  U6H  THF  COEFFICIENT  MATRIX  FOR  THE  FIT  OF  DEGREE  ,I1,70H  I 
IS  MACHINE  SINGULAR.  THE  PROGRAM  WILL  TRY  THE  OTHER  FITS  ASKED  FOR 
2.  ) 

FORMAT  U04HTHF  LIMIT  PARAMETER  IS  ZERO  FOR  THIS  GROUP.  THIS  MEANS 
1  NO  FIT  WAS  ASKED  FOR,  GOING  ON  TO  NEXT  GROUP,     ) 

FORMAT  <1H0»*3X.30HTHE  EQUATIONS  HAVE  THE  FORM  -- »//,45X,50HEST,Z 
1"  A1*A2(X)*A3(Y)*A4(X»«2)*A5(X) (Y).A6(Y»«2) ♦,/,53X,*7HA7(X«»3) *A8( 
2X»«2) (Y)*A9(X) <Y»»2) ♦A1Q(Y«»3) ♦     »/»53X«4  7HAn (X*»4) *A12 <X»»3)  < Y 
3)*A13(X»«2) (Y*«2)»  ,/,53X,47HAl4(X>  ( Y»»3) ♦ Aj5 ( Y<>*4 ) ♦ Al6 ( X*» 

45)*A17(X»«4)  <Y>*  i/»53X»47hAj8(X*«3) (Y»»j) ♦A19(X«*2>  (V*#3) *A20(X) ( 
5Y«*4>«   ,/,53X,10HA21 (Y*«5) .) 

FORMAT  (lHo,42X,22hIN  THIS  JOB  REQUEST  — ,//,4-?X ,27HTHE  AREA  BEING 
1  ANALYZED  IS  ,3A10«'»4lX«29HTHE  FACTOR  BEING  ANALYZED  IS  ,3A10./»* 
21X.49HTHE  HIGHEST  DEGREE  OF  TREND  SURFACE  REQUfsTED  IS  .12) 

END 


190 


195 

200 
205 

210 

215 
220 
225 
230 

235 
240 

245 


250 
255 


260 

265 
270 


275 


Z-LINEAR   RES-L 
RES-OR   Z-QUInT    RES-Q 


.116 

31S 

320 

322 

324 

3?6 

3?8 

330 

332 

334 

336 

338 

340 

342 

344 

346 

348 

3S0 

352 

354 

356 

358 

360 

362 

364 

D  366 

D  368 

D  370 

D  372 

D  374 

D  376 

I)  378 

D  380 

D  382 

D  384 

D  386 

D  388 

D  390 

D  392 

D  394 

D  396 

398 

400 

402 

404 

406 

408 

410 

412 

414 

416 

418 

420 

422 

424 

4?6 

428 

D  430 

D  432 

D  434 

D  436 

D  438 

D  440 

D  442 

D  444. 
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v)     FORSTAT    CROSS 

SUBROUTINE  RITtR 


-REFERENCE    LISTINGS    FOR    SUBROUTINE    RITER 


ENT  NUMBER 

DEFINITION 

5 

15. 

14. 

10 

21. 

18. 

15 

22. 

17. 

2n 

24. 

14. 

25 

29. 

27. 

30 

33. 

31. 

35 

35. 

14. 

♦  0 

38. 

36. 

4S 

40. 

14. 

SO 

41. 

40. 

55 

44. 

43. 

60 

45. 

40. 

65 

47. 

45. 

70 

53. 

48. 

75 

55. 

52. 

80 

57. 

79. 

85 

66. 

65. 

90 

69. 

64. 

95 

71. 

70. 

100 

72. 

70. 

ins 

74. 

72. 

110 

75. 

56. 

11"5 

80. 

76. 

120 

B2. 

14. 

125 

B6. 

84. 

130 

88. 

87. 

135 

92. 

90. 

1*0 

93. 

87. 

1*5 

95. 

14. 

150 

97. 

14. 

155 

99. 

14. 

160 

102. 

14. 

165 

106. 

FORMAT 

26. 

170 

107. 

FORMAT 

28. 

175 

108. 

FORMAT 

30. 

180 

109, 

FORMAT 

32. 

185 

110, 

FORMAT 

35. 

190 

111. 

FORMAT 

16. 

195 

112. 

FORMAT 

15. 

200 

113. 

FORMAT 

22. 

2ns 

114. 

FORMAT 

82. 

210 

115, 

FORMAT 

83. 

215 

116. 

FORMAT 

85. 

220 

117. 

FORMAT 

86. 

225 

118. 

FORMAT 

91. 

230 

119, 

FORMAT 

88. 

235 

120, 

FORMAT 

89. 

240 

121. 

FORMAT 

41. 

245 

122. 

FORMAT 

42. 

?bO 

123. 

FORMAT 

71. 

255 

124. 

FORMAT 

24. 

260 

125. 

FORMAT 

95. 

2*5 

126. 

FORMAT 

97. 

270 

127. 

FORMAT 

103. 

275 

128. 

FORMAT 

25. 

REFERENCES 


20, 


65. 


58.  59, 


60. 


37. 


99.       102. 


100.       104. 


55 

STATEMENT  NUMBERS 

DEFINED 

4 

MAXIMUM  NUMHER  OF 

RFFERENCES 

65 

TOTAL    NUMHER  OF 

REFERENCES 

1.2 

AVERAGE  NUMHER  OF 

RFFFRENCES 

•••FOF»»« 
•  Oi-L»FnF» 
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W) 


2. 

3. 

4. 
5. 
6. 

7. 
6. 
9. 


10. 
11. 
12. 
13. 

14. 

15. 
16. 
17. 
IB. 
19. 
20. 
21. 
22. 
23. 
24. 
25. 
26. 
27. 
28. 
29. 
30. 


31. 
3?. 
3J. 

34. 


LISTING    OF    SUBROUTINE    INVER 

SUBROUTINE  INVER  (k.ni)EGRE.KSOBI  B  2 

C    ......... — . ....  —  ..-.._....--...... —  .......   8  4 

C                           SUBROUTINE  INVER  B  6 

C    ............ .... — --. — ........... — .........................   fl  p 

C    PURPOSE"  B  10 

C          TO  COMPUTE  THE  INVERSES  OF  ALL  THE  XY-TRAMSPOSE.   XY  MATRICES     8  12 

C    FOR  ALL  THE  SURFACES  BEING  FITTED,  B  14 

C    ......... . .............................   R.  U 

c              the  gcahs  system  b  is 

c               Purdue  university  b  20 

C  B  22 

C                                                     A.  K,  TUHNERi  8  24 

C                                                 CIVIL  ENGINEERING,  8  2ft 

C                                                 JUNE  196*.  B  28 

C    ... — . — --   8  30 

C    USAGE--  B  32 

C               CALL  INVER  (XX, XINV, K,NDEGRE, KSOB)  B  3* 

C          WHERE  -  XX      a  MATRIX  FOR  WHICH  INVERSE  IS  REQUIRED  B  36 

C                  XINV    a  THE  CALCULATED  INVERSE  B  38 

C                   K        a  SIZE  OF  MATRIX  (NO,  OF  COLS,  ANU  ROWS)  8  40 

C                 NDEGRE  =  DEGREE  OF  SURFACE  8  42 

C                  KSOB    a  IF  XX  MATRIX  IS  SINGULAR  KSOB  SET  fQUAL  TO  ONE   B  44 

C                            AND  CONTROL  RETURNED  TO  TREND,  8  46 

C         NOTE-   THIS  SUBROUTINE  UTILIZES  THE  LABFLLED  COMMON.  B  48 

C                    /MATRIX/  B  50 

C    ..... ....................................................... — ...  H  53 

C    DESCRIPTION  OF  CONTROL  CARDS-  B  54 

C         NONE  REQUIRED  8  56 

C ...  - ............................. .......... — ..... — .__   h  58 

C    REQUIRED  SUBROUTINES-  8  60 

C         NONE  B  6? 

C    OPTIONAL  SUBROUTINES-  8  64 

C         NONE  B  66 

C    . ... ..... ........ ....................   b  68 

DIMENSION  XX(21,21),  XINV(21,21,5) ,  AM(?1,22>  8  70 

DIMENSION  N(21)  8  72 

OIMENSION  DUMK2643)  8  74 

COMMON  /MATRIX/  AM , DUM1 , XX , XINV  B  76 

00  5  I»l iK  B  78 

DO  5  Jel,K  8  80 

5      AM(I,J)aXX(ItJ)  8  82 

KK»K*1  B  84 

C    ............ . . . g  86 

C      INITIALIZE  THE  ROW  PERMUTATION  CELLS  TO  ZERO.  B  88 

C    ..... . R  qo 

DO  10  Iel.K  B  92 

10     N(I)«0  B  94 

DO  85  JaJ.K  B  9ft 

I"l  B  98 

JJ"J-1  B  100 

15        DO  20  L»1»JJ  B  102 

IF  (N(L)-I)  20.25.20  B  104 

20        CONTINUE  8  106 

IF  (AM(I.l))  35.25.35  B  108 

25        I"I*1  B  110 

IF  (K-I)  30.15,15  B  112 

30        KSOBal  B  114 

RETURN  B  116 

35        AMA»ABS(AM(I.l) )  B  118 

II»I  B  120 

*0        II»II*1  B  122 

IF  (K-I I )  60.45.45  B  124 

45        DO  50  IZ»ltJ  8  126 

IF  (N(IZ)-II)  50. 40. SO  B  128 

50        CONTINUE  B  130 

AMBaABSIAMIIIll) )  B  133 

C    ... .... ..  B  134 

C      TEST  FOR  MAXIMAL  SIZE  PIVOTAL  ELEMENT,  B  136 

C B  138 

IF  (AMA-AM8)  55.40.40  8  140 

55         AMAaAMB  B  142 

1*11  8  144 

GO  TO  40  8  146 

r  8  148 
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35. 

60 

C 

C 

c 

M(J)»I 

ADJOIN  THE  CORRECT  IDENTITY  COLUMN  VECTOR. 

36. 

DO  65  IZ»1.K 

37. 

65 

AM(IZ,KK)aO.O 

36. 

AM(I,KK)=1.0 

39. 

C 
C 

c 

DAaAH(Iil) 

NORMALIZE  THE  PIVOTAL  ROw. 

40. 

DO  70  IZs2,KK 

41. 

70 

AM(I,IZ-l)»AM<ItI2)/DA 

42. 

DO  85  IZ«1.K 

43. 

IF  (U-I)  75.85.75 

44. 

75 
C 

c 
c 

DBbAMUZiI) 

ZERO  OUT  NONXPIVOTAL  COLUMN  ELEMENTS. 

45. 

DO  80  JZ=2.KK 

46. 

80 

AM(I7.JZ-1)«AM(IZ.JZ)-AM(I.JZ-1)»DB 

47. 

85 

CONTINUE 

48. 

DO  100  laliK 

49. 

00  90  J»1«K 

50. 

IF  (N(J)-I)  90.95.90 

51. 

90 

CONTINUE 

52. 

95 

DO  100  L"liK 

53. 

LL«N(L) 

54. 

100 

XINV(L.I.NDEGRE)sAM(LL.J) 

55. 

RETURN 

56. 

END 

8  150 
B  152 
B  154 
B  156 
B  158 
8  160 
B  162 
8  164 
8  166 
B  168 
8  170 
8  172 
8  174 
R  176 
8  178 
8  180 
8  182 
8  184 
B  186 
B  188 
B  190 
8  192 
B  194 
8  196 
B  198 
8  ?00 
8  ?02 
B  204 
B  206 
8  208 
8  210- 


x)     FORSTAT     CROSS-REFERENCE    LISTING    FOR    SUBROUTINE    INVER 


SUmROUTINE  INVER 
STATEMENT  NUMBER    DEFINITION 


REFERENCES 


5 
10 
15 
20 
25 
30 
35 
40 
45 

50 
55 
60 
65 
70 
75 
80 

as 

90 

95 

100 


11. 
15. 
17. 
19. 
21. 
23. 

25. 

27. 

?9. 
32, 

35. 
37. 
4  1  . 
44. 
46. 
47. 
51  . 
5?. 
54. 


6. 

7. 

10. 

20. 

15. 

16. 

16. 

1H. 

20. 

lb. 

28. 

31. 

26. 

27. 

28. 

31. 

26. 

36. 

40. 

43. 

45. 

12. 

4?, 

49. 

50. 

50, 

48. 

5?. 

20 

STATFMENT  NUMBERS 

DEFINED 

3 

MAXIMUM  NUMBER  OF 

REFERENCES 

30 

TOTAL    NUMBER  OF 

REFERENCES 

1.5 

AVERAGE  NUMBFR  OF 

REFERENCES 
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y)       LISTING    OF    SUBROUTINE     POLY 


9. 
10. 
11. 

12. 


13. 
1*. 

IS. 
lb. 
17. 


lb. 


SUHROUTINt  POLY  (PT,XD,Yn,K) 

SUBROUTINE  POLY 


PURPOSE- 
TO  EVALUATE   ANY  DESIRED  TREND  SURFACE  EQUATION  AT  ANY  HEOUESTE 
X  ANO  Y  COORDINATES  IN  OROEH  TO  CALCULATE  THE  APPROPRIATE  ESTIMATES 
OF  Z. 


THE  GCARS  SYSTEM 
PURDUE  UNIVERSITY 


A,  K,  TURNER, 
CIVIL  ENGINEERING, 
JUNE  196*. 


USAGE- 
CALL  POLY(PT,XDtYD,K) 
rfHERF  -  PT  a  THE  OESlRED  ESTIMATE  OF  Z 
XD  a  X-COORDlNATE 
YD  a  Y-CUORDlNATE 

K   «  DEGREE  OF  EQUATION  BEING  EVALUATED 
NOTE-   THIS  SUHHOUTINE  UTILIZES  LABELLED  COMMON. 
/MATRIX/ 


DESCRIPTION  OF  CONTROL  CARDS- 
NONE  REQUIRED 


REQUIRED  SUBROUTINES" 

NONE 
OPTIONAL  SUBROUTINES- 

NONE 


DIMENSION  A(21,5>,  DUMMY(3000) 
COMMON  /MATRIX/  DUMMY, A 
GO  TO  (  10,15,20,25,30) ,  K 

FIRST  DEGhEE  COEFFICIENTS 

PT»A(1,K) ♦XO»A(2iK) *YD»A(3,K) 
GO  TO  35 


SECOND  DEGREE  COEFFICIENTS 


PT»A(1,?)*XD»A(2,2) ♦YD»A(3,2)*XD«XD«A(4,2> *XD»YD«A ( 5,2 ) *YD» YD* A (6, 
12) 
GO  TO  35 


third  degree  coefficients 


XD2«XD«XD 
YD2=YD»YD 
PT»A(1,3)*X0«A(2,3) »YD<»A(3,3)*X02*A(4,3) ♦XO»YO*A (5,3) ♦YOZ»A (6,3) *X 

1D2»XD»A(7,3)*XD2*YD»A(B,3) ♦XD«Y02»A (9, 3) ♦ Y02«YD«A ( 1 0. 3) 
GO  TO  35 


FOURTH  DEGREE  COEFICIENTS 


XD2«XD»XD 

XD3«XD2»XD 

YD2oYD»YD 

YD3aYD2»YU  - 

PTbA(1,4)*XD*A<2»4) ♦  YD«A(3,4)»X02»AU,4) »Xn»YD*A ( 5 ,4 ) ♦ YD2«A ( 6 , 4 > »X 
lD3»A(7,4)»X02«YD»A(8,4)*XD«YD2*A(9,4)*Yn3»A(10i4)*XU2««02»A(ll,4)» 
?XD3«YD«A (12,4) ♦XD2»YD2«A(13,4)*XD»YD3»A(14,4)«Y02*YD2»A(15,4) 

GO  TO  3S 


c 

? 

c 

4 

c 

6 

c 

8 

C 

10 

c 

I? 

c 

14 

c 

16 

c 

1H 

c 

20 

c 

22 

c 

24 

c 

26 

c 

28 

c 

30 

c 

3? 

c 

34 

c 

36 

c 

38 

c 

40 

c 

4? 

c 

44 

c 

46 

c 

48 

c 

50 

c 

5? 

c 

54 

c 

56 

c 

58 

c 

60 

c 

6? 

c 

64 

c 

66 

c 

6H 

c 

70 

c 

72 

c 

74 

c 

76 

c 

78 

c 

8  0 

c 

82 

c 

84 

c 

86 

c 

88 

c 

9  0 

c 

*?. 

c 

94 

c 

^6 

C 

98 

c 

100 

c 

1  0? 

c 

1  04 

c 

106 

c 

1  08 

c 

1  10 

c 

1  12 

c 

1  14 

c 

116 

c 

1  18 

c 

120 

C 

12? 

c 

124 

c 

126 

c 

128 

c 

130 

c 

132 

c 

134 

c 

136 

c 

138 

c 

140 

A-83 


19. 

20. 
21. 
22. 
23. 
24. 
25. 

26. 


27. 
2n. 


C ....  C  1*2 

C    FIFTH   DEuHtt  COEFICIfcNTS  C  1*4 

C —  C  j4ft 

30     XD2«XD*XU  C  148 

XD3«XD2«XD  C  150 

XD4oXD3«XD  C  152 

YD2"YD«YD  C  154 

Y03"YD2»YD  C  156 

YD4»Y03»YD  C  158 

PT»A<1,5)»XL>«A(2,5>*Y0«A(3.5)*XD2«A<4,5)*X0OYD»A(5.!5)*yD2»A(6»5>*X  C  160 

lD3»A(7.5)*XDi!*YD*A(d.5)  *XD*Y02»A(9.5)*Yn3<»A(10i5)  C  162 

PT»PT*XD4*A (11  .5) *XD3»YD»A (12,5) ♦X02»YD?»A (13.5) «XD*Y03«A (14.5) ♦ YD  C  164 

14»A<15.5) ♦XD4«XD«A (16.5) »XD4*Y0«A (17.5) ♦XD3*Y02»A( 1«.5) ♦XD2»YU3»A (  C  166 

219.5)»XP«Y04«A(20.5) ♦YD»YD4»A(21,5)  C  168 

35     RETURN  C  170 

END  C  1 72- 


z)    FORSTAT     CROSS-REFERENCE    LISTING    FOR  SUBROUTINE     POLY 


SUHROUTINE  poly 
STATEMENT  NUMBER    DEFINITION 


REFERENCES 


5 
10 
15 
20 
25 
30 
35 


4. 

5. 

7. 

<*. 
13. 
19. 
27. 


4. 
4. 
4. 
4. 
4. 
6. 


12. 


18. 


7 

STATEMENT  NUMBERS 

DEFINED 

4 

MAXIMUM  NUMBER  OF 

REFERENCES 

9 

TOTAL    NUMBER  OF 

REFERENCES 

1.3 

AVERAGE  NUM8FR  OF 

REFERENCES 
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2, 

3. 
4.. 

5, 

6. 

7. 

8. 

».. 
Id. 
1). 
12.. 
13. 
14, 
IS). 
16. 
I  /. 
If. 
19, 
20. 
21. 
22, 


cc)  LISTING  OF  PROGRAM  RESMAP 

"*"  PROGRAM  HESMAP(INHU'T»0UTwUT»TAPE3«P.UNCH«TAPE5"iNPUT»tAPE6»0UTPUT)  A  ? 

C — — — —  4  4 

C  _                               PROGRAM  RESMAP  A  ft 

C   — " "- A  8 

C   PURPOSE  -  a  in 

C ID  PRODUCE  RESIDUAL.  MAPi  FROM  TREND. SURFACE  GRID  VALUES.  A  12 

C  A  14 

C    A  \ft 

c  _            the  scars  system  a  is 

C                      PURDUE  UNIVERSITY  A  20 

C  A  22 

C                                                        A.  KEITH  TURNER,  A  2'* 

C                                                     CIVIL  ENGINEERING,  A  2ft 

C                                                     JUNE  IVfty.  A  ?8 

q   _ A  30 

C  A3? 

C    DESCRIPTION  OF  CONTROL  CARDS  -  A  34 

C  .___  A  36 

C    1 )  TITLE  CARD  A  38 

C          A  GENERAL  TITLE  IN  COLUMNS  1-80.  A  40 

C .........  4  42 

C    2)  OPTION  CARD  A  44 

C         NSURK  -  NUMBER  OF  TRENDS  FOR  WHICH  RFSIDUALS  ARE  TO  BE  A  4ft 

.C CALCULATED  .  (  COL  2  )  A  48 

C          NTAPF  -  IF  N1ARE=1,  RESIDUALS  WILL  RF  STORED  ON  1APF3.   (COL  4)  A  SO 

C          IRON   -  IF  IPUN=1,  RESIDUALS  PUNCHED  ON  CARDS  (COL  ft).  A  52 

C  A  54 

C    TMF  FOLLOWING  CARDS  ARE  REQUIRED  FOR  THE  ORIGINAL  DATA  ARRAY, AND  FOR  A  5ft 

C    EACH  TREND  SURFACE  UATa  ARRAY.   THE  ORIGINAL  DATA  ARRAY  MUST  6E  THE  A  58 

C    FJRST  ARRAY,  THE  TREND  SURFACE  ARRAYS  CAN  THEN  FOLLOW  In  ANY  ORDER,  A  60 

C  Aft? 

C  A  64 

C   A)  format  card  a  ftft 

C         FORM  (  COLS   1-80)  -FORMAT  FOR  DATA  ARRAY ( NORMALLY  (HF10.3)).  A  68 

C  A  70 

C    8)  HEAOFR  CARD  A  72 

C          FACTOR  (COLS   1 -40 ) -ALPHANUMERIC  IDENTIFICATION  UF  FACTOR.  A  74 

C         AREA    (COLS  41-701-AI PHANUmERIC  IDENTIFICATION  OF  AREA.  A  7ft 

C          NCOL     (COLS  71-73)-NIIMEFR  OF  COLUMNS  In  THE  ARRAY.  A  78 

C         NROW    (COLS  7<»-76)-NilMhER  OF  HOWS  IN  THE  ARRAY.  A  80 

C          HYPE   (   COLS  76  )-IF  DATA  IS  TOPOGRAPHIC  RESIDUALS , I  I YPE  =  1  .  A  8? 

C                               IF  DATA  MEASURES  COSTS  •  ITYRE=?  .  A  84 

C                                 IF  DATA  MEASURES  BEnEF I TS » I T YPEa3  ,  A  8ft 

C  A  88 

c   o  data  cards  a  go 

C  A  <J? 

C a  94 

C    OPTIONAL  TAPF  -  A  9ft 

C         TAPE3  may  BE  REQUIRED  .  A  48 

C  A  1  00 

C   a  jog 

C    REQUIRED  SUBROUTINES  -  A  104 

C          -  NONE  A  10ft 

C    OPTIONAL  SUBROUTINES  --  A  108 

C         -  NONE  A  110 

C  All? 

C A  114 

DIMENSION  F0Rm(8),  FACTORU),  AREAI3),  RES  (  3 )  ^OR  (  50  ,  5")  ,  Z  (  50  ,  SO)  ,  A  lift 

1  T ITL  (8)  A  1  18 

INTEGER    SYSIN,SYSOT»SOF3,M.iRM  A  121 

DATA    SYSIN.SYS01 , SDF3/5 , ft • 3/ , RES ( 1 ) / 1  OH    RESIDUALS/  A  12? 

5               RMO     (SYS1N,70)     TJTL  A  124 

IF      (FOF.SYSIN)     60,10  A  12ft 

1(1             READ     (SYSIM,65)     NSURF ,  M  APF  .  I  Hi  IN  A  128 

hFAO     (SYSlN.70)     FURM  A  130 

WRITE    (SYSDT.75)     HTL  A  132 

READ     (STSlN,90l     FACIOR,AREA,NCOL,NROw, ITYPF  A  134 

WRITE     (SYSOT18O)  A  l 3ft 

DO    15    I=l,NROW  A  138 

READ     (SYSIN.FOHM)      (ZOR(T.J) »J=l«NCOL>  A  140 

15           WRI1F    (SYS0Ti85)     (ZOR <I» J) , Jal.NCOU  A  14? 

DO    55    K  =  l, NSURF  A  144 

READ  (SYSIN.70)  FORM  A  14ft 

READ  (SYS1N,90)  F ACTOR , ARE A , NCOL • NHOW , I T YPE  A  148 

."RITE  ISY.S0T»9S)  FACTOR  A  150 

no  2(i  I  =  i  ,nrow  A  l  5? 

RfAD  (SYSIN.FOHM)   (Z(IiJ) tJBltNCOLl  A  154 

<?U  .     -   "RITE  (SYSOT.ba)   1 1.  (  I .  J)  ,  J= )  ,NCOL)  A  156 

"•RTfi-'  (SYSOT.lUO)  A  1  5« 
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24. 
25. 
26. 
27. 
26. 
29. 
30. 
Jl. 
32. 
33. 
}4. 
3b. 
36. 
37. 

38. 

J9. 

4(1. 

41. 
42. 

43, 
44. 
45. 

46. 
47. 


00  30   I«l»NHOW 

U0  2b    J«1»NC0L 
2b  Z<ItJ)=ZORU«J>-Z(I»J> 

30        WHITF  (SYSOT.Hb)   U ( I . J) . J=l .NCOU 

IF  (nIAPE.NE.1)  60  TO  40 

rfRITE  (SDF3)  FwCIOH.RFS.NCOL.fvHUW,  ITYPE 

00  3b    I=ltNROW 
Jb        WHITE  (SDF3)  (2(IiJ)»JenNC0L) 
40         IF  ( [PUN.NE. 1 )  60  10  50 

PUNCH  90?  FACTOH.RtStNCOL.NHOW.ITYPE 

DO  4b  I=1»NR0W 
4b        PUNCH  FORM.  (2  (  I.  J)  ,J=1  ,NCO[.) 
50         CONTINUE 
bb     CONTINUE 

60  TO  5 
C 

60     STOP 
C 

6b     FORMAT  (31?) 
70     FORMAT  ( HA  1 0) 
7b     FORMAT  (1H1.8A10) 

BO     FORMAT  (1H0»27H  THE  ORIGINAL  BRIDGED  DATA  ) 
ab     FOHMAT  (  1H0.12F10.3./ (]X, 1 ?  F 1 0 . 3 ) ) 
90      FORMAT  (4A1 Oi3A10i2I3t 12) 
9b     FORMAT  (lhli4A10) 

100    FORMAT  (lHlf22H  THE  liRIDOEl)  RtSIOUAUS) 
FNO 
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1  60 

A 

16? 

A 

164 

A 

1  66 

a 

)68 

A 

1  7'i 

A 

1  7? 

A 

1  74 

A 

i  rt, 

A 
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